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ABSTRACT 


A study is made of the guiding properties of a long electrostatic 
quadrupole that has been uniformly twisted about its axis. Analytic 
solutions have been developed for particle trajectories and they show 
that stable orbits exist if the structure parameters are properly 
chosen. A confinement criterion, which sets the acceptance limits of 
the twisted structure, has been derived using the analytic solutions. 

The accuracy of the analytic solutions and the validity of the confine- 
ment criterion are examined by comparing the analytically computed 
trajectories with those computed by numerical integration. The acceptance 
limits of the twisted quadrupole structure are compared to those of the 
classical electrostatic quadrupole channel. Twisted electrostatic 
quadrupole structures of finite lengths are found to be capable of 
identically reproducing the injection conditions of a system of particles. 
Such lenses have been fabricated and their imaging properties have been 
tested using a beam of contact charged microparticles. The experimental 
results are in good agreement with the theoretical analysis. A study 

has also been made of the effect of bending the axis of the twisted 


structure, and it was found that stable orbits can exist. 
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CHAPTER 1 


1.1 Background 


During the past few years an increased interest has developed 
in the acceleration of very heavy charged particles. Serious efforts 
are now being made to ionize and accelerate elements as heavy as 


lee) The interest in obtaining energetic heavy ions is largely 


argon 
directed towards the emerging field of nucleosynthesis. However, the 
interaction of such particles with atomic, molecular and biological 
systems promises to yield much vital information in such areas as 

the study of radiation damage in solids and the investigation of 

cellular disorders. Also recently, attention has been focused on charged 
micron-sized solid and liquid particles. When accelerated to high 
velocities the applications for such particles range from sprayed ink 

(6) 46 mass expulsion systems for 
(8) 


sized particles moving at velocities of the order of 10° ki lometers/sec 


jets for extremely high speed printing 
Spacecraft oi It has even been suggested that micron- 
could trigger thermonuclear reactions in a controllable manner. A 
further use of these macroscopic particles, which is of particular 
interest at the University of Alberta, is for the simulation of 
micrometeoroids with velocities up to several tens of kilometers per 


second. Such artificial micrometeoroids can be used to study impact 


- 5 a 


= 


ities 
a Oe 


_ . 


all 


a | 
"fen “sn! * dae Lal 


> 


- 7 26. uve Sif 4 Pea i 5; 4 he pela af, one axtay) at abes priiod wom 9 : 
| é- 

eon’ pl) 2t-2ne! yer 3! se yone pel ararde Pity jtesraial ont A ty. agp 
_ airs 4 13¥ae0! 2 tpattacenaisui table V) qitpre ae wid ote war sala 


7 Tavtpolots baa “taluoeloem ,oiwnds AdiW zelotgasg Rope TH net foe 1810 


7 7 
: —_ 
pagotqech es teaustai horesyec! mp 2am “a 3269 a0) ante 


etaeFis\suotse $=. tal alt vag: bayintia was vray %0 notsersiezs sft 


7 


| : te peas Hout af nol saoridini fotiv mou Olary G2} zoninawd Bmes2y = 
WO aotsaphi maya? on. one | toe “? apemnh aotyeraar VO bude a3 

: xy tne wo Wareoot ew! cad motiaet*o. ys tee orth .erabyoe te set i 7 
7 

: 7 fufi-et batewsiaone wets celaizoeg Giuphl die bilog beste-now ‘a 

7 mt bsyerge wr) apnsy estaltver nove vot geettaabloqge eda 2eds hse a 

“Qe? agmgeve aol 2! vors zene 02! Masten “gq Pete Ape cfeamayg <9 we? eval J 
| —— sani Wein ned love cot a1 .\“notatnaore 3 Detsesege 


ad 


_ _ paphevstonat aT 70 vabig ost to esttfoaley te entvem ealatevag b 


sia . a? crane fdel fov80a> 5 ai enotioge Wstaumerens veppend wes 
: peng. 
(sali vag 
: _ a ore 
: wiiaiel \e? et pprvaalA 2 vItenovinl! ote 25 seoregAt 
- ? 
- : ai her Jwvevee 02 ap 283 taotav riz any oom 


on nes ebtovmatenostin felshtit 


i“ 


— Aahitw p29(atiasq oldoseoTseN 5 aenit Ye 92u* 


7 —_ 


phenomena and the resulting hazards to men and equipment outside the 
earth's atmosphere. For example, high velocity aluminum particles, of 
approximately 4 micron diameter, have been used to investigate the 
abrasion of optical surfaces exposed to micrometeoroid pomeaudnean 2. 
Radio-frequency linear accelerators have been suggested as the 
most suitable machines for producing very energetic heavy particles! !9) . 
In such accelerators an r.f. source is used to establish an electric 


field across a linear array of accelerating gaps, as shown schematically 


im figure (1.1). 
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Figure (1.1) Accelerating Gaps and Drift Tubes of a 


Heavy Ion Linear Accelerator 


As particles pass along the axis of the structure they experience an 


essentially field free region within each drift tube. However, the 
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dimensions of the drift tubes and of the accelerating gaps are so 
adjusted that the particles reach each successive gap at a time when they 
- will be accelerated by the electric field. The field across a typical 


accelerating gap is illustrated in figure (1.2). 


a electric field lines 
taal 


fd 


Figure (1.2) Field Distribution Across a Typical Accelerating 


Gap 


Note that in the first half of the gap the radial component of the 
electric field will exert a force on particles that is towards the axis 
of the accelerator. In the second half of the gap the particles will 
experience a defocusing force away from the axis. Since the accelerated 
particles will spend less time in the second half of the gap, one might 
imagine that the net effect would be a focusing force action tending to 
keep the particles near the axis of the structure. Unfortunately, to 
ensure longitudinal stability of the motion, particles must pass through 
the gap when the amplitude of the r.f. field is iAcrescing This 
increase in the field is sufficient that the defocusing forces are 


increased to the point where the net effect of the accelerating gap is 
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4 
to deflect particles away from the axis of the machine. Space charge 
forces, stray fields, and misalignments have the effect of further 
_defocusing the particles. In a long accelerator it becomes essential to 
compensate for these defocusing actions as well as to make provisions 
for handling particles that enter the structure with non-zero transverse 
momenta. 

In electron and proton accelerators this compensation is almost 
invariably provided either by axial magnetic fields or by magnetic 
quadrupole lenses. Such focusing devices work well as long as the charge- 
to-mass ratios of the particles are large. However, as the charge-to-mass 
ratio, i decreases, the required magnetic fields can become excessive. 
Even the lightest of heavy ions have values of charge-to-mass ratios 
that are several orders of magnitude smaller than that of an electron. 
Singly charged argon for example has a charge-to-mass ratio of only 
Bese 10° coulombs/kilogram. The values of y for typical solid or 
liquid microparticles are much smaller still and may range from as smal] 
as 0.001 coulombs/kilogram to an upper limit of approximately 500,000 


(12-13), The use of magnetic fields to guide and contain 


coulombs/kilogram 
beams of such heavy particles is quite impractical. Consider, for example, 
the problem of using an axial magnetic field to confine a particle with 

a transverse velocity of u = 1000 meter/sec. The particle has a charge- 
to-mass ratio of 1000 coulombs/kilogram and is to be confined within a 


cylindrical aperture of radius, R = 2 cm. The required magnetic field is 


easily calculated from 
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In this case B = 500 kilogauss. At the present state of technology, 


magnetic fields of this magnitude are almost impossible to maintain, 


nee Moreover, the ohmic losses associated 


‘especially over large apertures 
with the large magnet currents are prohibitively high. The power 
requirements can be reduced using superconducting magnets but only at 
the expense of costly cryogenic installations. On the other hand, electric 
field focusing devices consume almost no power and the required fields 
can be established in reasonably sized structures at modest voltages. In 
the preceding example, an electric field of only 500 volts/cm is needed 
to produce a force equal to that exerted on the particle by the magnetic 
field. 

The purpose of the present thesis is to present the analysis, 
design and testing of a novel electrostatic system for guiding and focusing 
very heavy charged particles. It is intended that the structure be 


compatible with the micrometeoroid simulator under construction at the 


University of Alberta. 


ince sone Existing sb lectric Field Focusing Devices 


In conjunction with the development of heavy ion linear accelerators 
two approaches have been proposed that use electric fields to provide the 


(15) | The first approach is to simply modify 


necessary focusing forces 
the radio frequency fields in the accelerating gaps so that the net effect 
of the gaps will be a focusing action, Such modification of the r.f. 

fields can be achieved through fitting the successive circular drift tubes 


with foils, grids or fingers or by using rectangular drift tubes which 
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6 
are rotated 90° with respect to each other. In the second approach the 
focusing action is provided by external electrostatic focusing devices 
~which are placed within the drift tubes. These devices almost always 


consist of some form of electrostatic quadrupole. 


Te2 nee On lSeandeur ds 

When the entrance to the right hand drift tube of figure (1.2) is 
closed by a metal foil thin enough for the beam to penetrate without 
damage, electric field lines in the accelerating gap are as shown in 


Figure (1.3)'19), 


electric field lines 


Figure (1.3) Field Distribution Across a Gap Between Drift Tubes 


with a Thin Foil on the Entrance of the Right Tube 


It can be seen that the defocusing radial component in the second half of 
the gap has disappeared and the whole of the gap is therefore convergent. 
In practice it has been found that the foils cause unacceptable scattering 
of the beam and cannot withstand even occasional breakdown between the 
drift tubes. Moreover, the ability of the foils to dissipate thermal 


energy is very low and the beam density is therefore extremely limited. 
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It was then found that fitting the entrance of the right hand side 
drift tube in figure (1.2) with a metal grid, which closes only a portion 
of the aperture, modifies the r.f. electric field in a manner similar to 
that of the foil. Several grid configurations have been suggested and 


used eieceseruliy ae 


1.2.2 Rectangular Drift Tubes, Rotated Grids and Finger Structures 


To obtain focusing in the radial direction while maintaining the 
longitudinal stability of the motion, linear accelerators with drift tubes 
of rectangular cross-sections have been Piposed ve) 17) The successive 
rectangular drift tubes are rotated 90° with respect to each other, as 


shown in figure (1.4a). 


Figure (1.4a) Rectangular Drift Tubes 
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Figure (1.4b) Electric Field Distribution Across the Gap 
Between Rectangular Drift Tubes in the X-Z and 


the Y-Z Transverse Planes 


As illustrated in figure (1.4b), the shape of the lines of force of the 
r.f. field in a gap is such that there is focusing in one transverse 
direction and defocusing in the other transverse direction. This action 
is reversed from one gap to the next. Thus, as the particle travels 
along the accelerator it will alternately experience focusing and 
defocusing forces, in both the X-Z and the Y-Z planes of figures (1.4a) 


and (1.4b). With proper choice of dimensions and accelerating voltages, 
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this alternating force action can lead to net focusing in both transverse 
directions. 

A similar alternating force can be obtained when using circular 
drift tubes whose entrances are fitted with two parallel wire grids. 
From one drift tube to the next the direction of the grids is altered by 


90°, as shown in figure (1,5) 629) | 


Figure (1.5) Circular Drift Tubes with Two Parallel Wire Grids 


Another structure that can provide an alternating force can be 
constructed by fitting the entrances and the exits of the circular 


drift tubes with fingers as shown in figure (1.6) (21) 


The difference 

in potential between the successive drift tubes,and hence between the 
fingers will exert focusing forces on a particle traversing the gap in 
One transverse direction and defocusing forces in the other transverse 
direction. Across the next gap the action is reversed, hence the particle 
will experience, in either one of the two transverse directions, 


alternating focusing and defocusing forces. Proper choice of the 


dimensions and the accelerating voltages again leads to a net focusing 
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Figure (1.6) Drift Tubes with "Fingers" 


action in both transverse directions. 

The alternating force action provided by the structures described 
above is known as alternating gradient focusing. Such action can also be 
provided using external focusing devices as will be explained in detail 
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1.2.3 Weak Electrostatic Lenses 

Weak electrostatic lenses consist of two or more cylindrical 
electrodes which are at different constant potentials. The focusing 
action is provided by the weak radial components of the d.c. electric 
field in the gap or gaps between the electrodes together with the small 


accelerating or decelerating action of these gaps. One of the weak 
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Figure (1.7) A Unipotential Electrostatic Lens 


electrostatic lenses used in linear accelerators is the unipotential 


lens which may be in the form shown in figure (1.7) 642), The beam 
entering the lens is allowed to spread out and then is bent toward the 
lens axis, with the focal spot near the end of the lens. To provide 
such action, the electrode potentiais should be comparable to the beam 
accelerating potential. Such weak lenses often are used for shaping the 


input beam for heavy ion secelentore 


Although this type of lens does not have widespread applications 
for guiding heavy particles it is widely used for shaping and focusing 


electron beams in electronic devices such as the cathode-ray tube or the 


electron microscope. 
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1.2.4 Electrostatic Quadrupole Lenses 


The basic principle of alternating gradient focusing, mentioned 


(24) as it applied to magnetic 


in Section 1.2.2, was first introduced 
quadrupole focusing. The principle was later extended to electrostatic 
quadrupoles. Since then, quadrupole structures have become an almost 
universal feature of beam transport systems. Since the principle of 
operation of the classical quadrupole structure is quite similar to that 
of the lens system described in this thesis, its operation will be 
reviewed in some detail. 


The fields of a quadrupole section are provided by four symmetrical 


electrodes, as shown in figure (1.8). 


Figure (1.8) Symmetrical Hyperbolic Electrodes 


The operation of an electrostatic quadrupole is based on having the 


potential distribution within the aperture in the form 
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which gives the fields as 


= = 2kX 5 


and Ee aacky 3 Cle} 


Such a potential distribution is obtained when hyperbolically shaped 
V 
—+ where V 


electrodes extending to infinity are used, in which case k = 5 
(25) 


0) 


wn Dw 


is the electrode focusing voltage and a is the aperture radiu 
Extensive investigations of the fields from electrodes of shapes other 

than hyperbolic have shown that electrodes which are parts of cylindrical 
cylinders give much the same field variation as that of equations (1.2), 
over a limited portion of the aperture. More complex electrode shapes 

also have been investigated and have been found to give suitable fields 

over limited apertures). 

Equations (1.2) show that each quadrupole section will provide 
focusing forces in one transverse direction and defocusing forces in the 
other transverse direction. To provide focusing action in both directions, 
doublet, triplet and multiplet combinations are used. These classical 
combinations have a series of sections placed end to end with the potentials 
reversed on alternate sections. Thus, a particle travelling along any of 
these combinations will be acted upon, in both transverse directions, by 
alternate focusing and defocusing forces. Through proper choice of the 


parameters the net action on the particle in both transverse directions 


can be focusing. 
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Figure (1.9) Symmetrical Quadrupole Triplet 


Figure (1.9) shows a symmetrical quadrupole triplet which consists 
of a straight quadrupole section of length 2/2 followed by a drift space 
of length d, then a second straight quadrupole section of length 2 on 
which the potentials have been reversed with respect to the first section. 
The reversal of the potential is equivalent to having rotated the second 
section through 90° while leaving its electrode potentials unchanged. The 


second quadrupole is followed by another drift space of length d, and a 


ms} — 
et 2s. 2 “~ 


om 


: hae pe clols <« ee — cai toed s— 


nd 


- ; 


sefqint staquideul testesemme (€,1) awit 


a 
a ws Nenen ste patatas a) sau ine ins [a dermmee e auane te, f) ang 
. oat, suckte® As Age! to wortwse sfoquvhaup tiple # fe 
: £950 Jlniind ® ‘gh beoliot S15 tens: TORI Ste Ah A 
7 : . snei to 


: J =~ 
mi erst %o. mal tise Sfoqueticun titpletts: Bae a sare 2) Sgee 
J. a a WV alt of vosazes yi hazevey nautl ayad 20eh203 oy ate OTM 
7 : _ a ot a a - - : . ne at — 43 20a , 
a 6f3 bit von wi anateyione 2h f on Se al 
7, Vay i) 1a) es = 
pilin Sel snecon aioives!s 31 pwivest atid “OC HpORM mE 
_— a vial | ; a ; 
Aiewat ForacnGe rtisb redsono ya banca af i ay c 
; 


a 
= : 
a 
7 


- 


15 
third quadrupole of length sewhich is rotated a further 90° with respect 
to the second quadrupole section. The first quadrupole section is 
- oriented such that its positive electrodes are aligned symmetrically 
about the X-axis of a fixed (X,Y,Z) coordinate sytem. Thus, the forces 
exerted on a positively charged particle at Z = O are focusing forces in 
the X-direction and defocusing forces in the Y-direction. Therefore, the 
X - Z plane can be referred to as the initially focusing plane of the 
triplet and the Y - Z plane can be referred to as the initially 
defocusing plane of the triplet. 

If the axial momentum of a particle travelling along the triplet 
is assumed constant, the projections of the particle motion on the X - Z 
and the Y - Z planes are uncoupled. Thereforesthese projections can be 
independently analyzed. In analyzing the motion it must be noticed that 
the basis of operation of the quadrupole combinations as guiding and 
focusing structures is the linear increase of the forces, focusing or 
defocusing, with the increase of the distance from the Z axis. 

Consider the motion of a particle along a triplet with properly 
chosen parameters. A particle injected in the X - Z plane will experience 
focusing forces during its travel along the first quadrupole section. 
These focusing forces will decrease the X directed transverse momentum 
of the particle and will tend to deflect it towards the Z axis. Due to 
the action of these focusing forces, the particle, after coasting across 
the first drift space, will enter the middle quadrupole section at a 
point closer to the Z axis than if the focusing forces were absent. Hence, 
the defocusing forces that the particle experiences as it travels along 


the second quadrupole section are smaller than they might have been if 
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the focusing forces of the first section were absent. The defocusing 
forces of the middle quadrupole section will deflect the particle away 
from the Z axis so that the particle, after traversing the second drift 
Space, will experience large focusing forces as it travels along the 
third quadrupole section. Thus, the particle will experience alternating 
focusing and defocusing forces as it passes along the triplet. However, 
if the triplet parameters are properly chosen, the particle finds itself 
farther displaced from the Z axis while experiencing focusing forces 

than when experiencing defocusing forces. Since the fields and hence the 
forces increase linearly with the distance from the Z axis, the effect of 
the alternating forces will be a net focusing action on the particle. 

A similar action occurs when a particle is injected in the Y - Z 
plane. The particle will first experience defocusing forces that wil] 
deflect it away from the Z axis so as the particle travels along the 
middle quadrupole section it will experience large focusing forces. These 
large focusing forces will tend to deflect the particle towards the Z 
axis so that it experiences small defocusing forces as it travels along 
the third quadrupole section. It has been found that the choice of 
parameters that produces net focusing action in the X - Z plane will also 
produce net focusing action in the Y - Z plane (27) , 

Most beam transport and focusing systems in accelerators are chains 
of triplets similar to that in figure (1.9). A chain of symmetrical 


triplets, referred to as a classical quadrupole channel, is shown in 


figure (1.10): 
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Figure (1.10) A Classical Quadrupole Channel 


An infinitely long classical channel will consist of a straight quadrupole 
section of length 2/2 followed by a succession of drift spaces of length, 
d, and quadrupole sections of length, 2. The successive quadrupole 
sections will be rotated 90° with respect to each other. The mathematical 
treatment of the particle motion along such a classical channel will be 


outlined in Chapter 6. 
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Recently, studies have been carried out for incrementally rotated 


magnetic quadrupole channels for which the angle of rotation between the 


180° 
7 where 


successive straight quadrupole sections was equal to 


tee 2h 165 2, 28-29), These studies showed that the focusing 


strength of the structure is improved as the angle of rotation between 
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the successive straight quadrupole sections is decreased. Preliminary 
investigations of continuously twisted magnetic quadrupole structures 20-3!) 
‘have shown that such structures provide better focusing action than either 
the classical or the incrementally rotated magnetic quadrupole structures. 

The purpose of the present study, as mentioned in Section 1.1, is 
to develop a guiding and focusing system for heavy charged microparticles. 
It is intended that the structure be compatible with the micrometeoroid 
simulator under construction at the University of Alberta. The focusing 
and guiding system must be an electrostatic structure if it is to be 
capable of guiding the artificial micrometeoroids whose charge-to-mass 
ratios are very low. In view of the results obtained by previous 


(30-31) during their investigations of the guiding properties of 


authors 
uniformly twisted magnetic quadrupole structures, it was expected that 

a uniformly twisted electrostatic quadrupole structure would provide good 
focusing and guiding of heavy particle beams. Therefore, it was decided 
to investigate, both theoretically and experimentally, the guiding and 
focusing properties of a uniformly twisted electrostatic quadrupole 
structure designed to handle beams of charged microparticles with very 

low and variable charge-to-mass ratios. 

The twisted quadrupole structure consists of four hyperbolically 
shaped electrodes which have been uniformly and slowly twisted about the 
structure axis as shown in figure (1.11). The radius of the circular 
aperture of the twisted structure is much smaller than the periodic length 


of twist of the electrodes. Each two nonadjacent electrodes are electrically 


connected. Two nonadjacent electrodes are at a positive potential, + Vo: 
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Figure (1.11) The Uniformly Twisted Electrostatic Quadrupole 


Structure 


while the other two are at a negative potential, - Vo: Thus, the force 
exerted on a positive particle on the axis of the first electrode pair 
is a focusing force while the force exerted on a particle an the axis of 
the second electrode pair is a defocusing force. 


In classical or incrementally rotated quadrupole structures the 


change in the direction of the focusing and defocusing forces occurs 
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suddenly as the particle leaves one straight quadrupole section to enter 
the next. In the twisted structure as the particle travels along the 
structure, the particle experiences forces that change gradually from 
focusing to defocusing and so on. However, as in the case of the 
classical triplet, the particle while experiencing the focusing forces 
will be farther displaced from the structure axis than while experiencing 
defocusing forces. Thus, the net action of these gradually alternating 
focusing and defocusing forces will be a focusing action and the twisted 
structure will be capable of guiding the particles. 

Beside its present application as a guiding and focusing system 
for heavy charged microparticles, the twisted electrostatic quadrupole 
structure could also be used to guide electrostatically sprayed liquids 
used in non-impact electrostatic ea as well as for guiding 
exhaust streams of colloidal propulsion peices: 

The work of this thesis is described in the chapters which follow. 
The theoretical analysis of charged particle motion along the slowly and 
uniformly twisted electrostatic structure is developed in Chapter 2. 
Then, Chapter 3 presents the derivation of the acceptance limits of an 
infinitely long structure and these limits are discussed for some special 
cases of injection. To examine the accuracy of the analytic solution 
developed in Chapter 2, particle trajectories computed using the 
analytical solution are compared, in Chapter 4, with the trajectories 
computed using numerical integration of the original equations of motion. 
The effect of earth's gravity on the motion of the heavy microparticles 
is studied in Chapter 5. The guiding properties of the twisted electro- 


static quadrupole structure are compared,in Chapter 6,with those of a 
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classical electrostatic quadrupole channel. 

A theoretical study of the imaging properties of a twisted 
electrostatic quadrupole structure of finite length is presented in 
Chapter 7, showing that the structure acts as a novel identical imaging 
device, provided the structure parameters are chosen properly. This 
imaging device is eaeeie of identically imaging the injection 
coordinates and momenta of a beam of particles and the imaging occurs 
with very small aberration. The imaging properties of a structure 
constructed from two equal but oppositely twisted electrostatic 
quadrupole sections are also studied in Chapter 7, showing this structure 
provides identical imaging with less aberration than a single twisted 
section. The theoretical studies are followed by experimental 
investigations of the imaging and guiding properties of the twisted 
channel. A high density beam of charged 3 micron carbonyl iron micro- 
particles was used in the experimental investigations. The descriptions 
and the results of the experimental investigations are in Chapter 8. 

Finally, Chapter 9 presents the analysis of the motion of a 
charged particle as it travels along a twisted electrostatic quadrupole 
structure whose axis is bent in a circle. This analysis is relevant to 
beam bending devices and to storage rings of heavy charged particle 


beams. 
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CHAPTER 2 
PARTICLE DYNAMICS IN THE UNIFORMLY TWISTED ELECTROSTATIC QUADRUPOLE 
STRUCTURE 


The structure consists of a long electrostatic quadrupole whose 
four hyperbolic shaped electrodes have been uniformly twisted about the 
structure axis at the rate of 8 radians per meter. The objective of this 
chapter is to determine the analytic solutions for the motion of charged 
particles travelling along the structure so that the guiding properties 
of the structure can be easily established. Only non-relativistic solutions 
are considered since the main objective is to provide a guiding system for 
microparticles, which simulate micrometeoroids whose velocities range from 


32) 


a fraction of a kilometer per second up to 72 kilometers per fered: : 


2.1 Coordinate Systems and Potential Distribution 


Figure (2.1) shows the four hyperbolic shaped electrodes of the 
structure. The distance a is the radius of the quadrupole aperture in 
meters and Is is the magnitude of the d.c. focusing voltage applied to each 
of the electrodes in volts. The polarities of the d.c. focusing voltages 
are shown in figure (2.1). In the analysis of the charged particle motion 
along the structure two rectangular coordinate systems are used, namely 
the (X,Y,Z) and the (x,y,z) coordinate systems. 


The quadrupole structure axis is taken as the Z-axis of the fixed 
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rectangular coordinate system (X,Y,Z). The axes (X,Y) coincide with the 
electrode axes at the input plane of the structure, Z = 0. The (x,y,z) 
coordinate system is a rotating coordinate system whose axes (x,y) 
rotate, as z increases, at the rate of 8 radians per meter so that the 
axes (x,y) always coincide with the electrode axes at any distance z. 
The Z and z axes of the two coordinate systems coincide. Figure (2.1) 
shows the orientation of the quadrupole electrodes and the axes (x,y) 
with respect to the fixed axes (X,Y) at a distance Z from the input plane 
of the quadrupole structure. 

It is shown in Appendix A that if the electrodes are slowly twisted 
so that the pitch length L = oT is much greater than the aperture radius a, 

2 


that is (ga)~ << 1, then for hyperbolic shaped electrodes and for points 


within the aperture, the potential distribution is given by 


V = k[tx2 - v2}cos 2827 + {2XY}sin 26Z] (2.1a) 
or 
V = k{x? - y?} (2.1b) 
V 
where k = ou : 
a 


Equation (2.1b) shows that at any cross section of the twisted 
structure the potential at points within the aperture in terms of the 
rotating coordinates (x,y,z) is simply that of a conventional quadrupole 


and is unaffected by the twisting action provided that 
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2.2 The Hamiltonian and the Equations of Motion 


Particle dynamics in the uniformly twisted straight electrostatic 
quadrupole can be conveniently discussed in terms of the Hamiltonian 
of particle motion, H. In the present case, H is a constant of the 
motion since the potential and hence the fields do not depend explicitly 
on fine oo In terms of the fixed coordinate system (X,Y,Z), H is 
simply the sum of the kinetic energy of the particle, T, and its potential 
( 33) 


energy » GV, where V is the electric potential given by either equation 


(2.1a) or equation (2.1b). Thus the Hamiltonian is 


a 
i} 


T Fray 


or, 


2 
if 


mae i 2 2 
H = oy [Py +P, + PL te X igz) (2.8) 
where q is the charge of the particle in coulombs and M is the mass of 
the particle in kilograms. Py, Py and Po are the mechanical momenta 


canonically conjugate to the (X,Y,Z) coordinates and are given by 


Py = MX ; 
Pe tgite : 
Ps = MZ : (2.4) 


where the dot notation represents the total time derivative. 
Substituting for the potential V, using equation (2.1la), H is 


given by 
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+ gk[{x* - Y*}cos 2e2 + {2XY} sin 28Z] (2.5) 


Using the canonical equations of Bamaton: os the equations of motion 


in the (X,Y,Z) coordinate system are found to be: 


X= - 2ok [Xecos 262 4 ¥esinezeZ| : 


Woke oak [X sin 282 - Y cos 28Z] 
eee £OK eLt2Xy} cos 267 ux> = Y-tsin ez) © (2.6) 


These are three coupled nonlinear second order differential equations for 


which no direct analytical solutions seem to exist. 


At this stage the problem can be simplified by expressing H in terms 


of the rotating coordinates (x,y,z) and the corresponding conjugate momenta 


(p,. Pye p,). The coordinates of the fixed and the rotating coordinate 


systems are interrelated through the following equations 


X = xX COS BZ - y Sin BZ A 
Vee xerSatienz 4 .V) COS BZ ; 
Bt=07 , (287) 


Differentiating with respect to time gives: 
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P 
X= = £X CoSy8z -'¥rsin Bz)’ = Bzix sin Bz + Y cos sz} : 
: Py 
Meee ON Ee yy COsmEz st Bz{x cos Bz - y sin 8z} ’ 
P 
oe a (2.8) 


Using equations (2.8), the kinetic energy of the particle is given by 


Spit: 1232) shale, tha? 
T 7 [P, + Py + Ps ] 
T= Foxe + y? + 282{xy - yx} + ge 2" {x2 + y*} + sal (2.9) 
Since the potential energy, qV = qk (x2 - y*), is velocity independent, the 
momenta (ps Py p,) adjoint to the rotating coordinates aan 
) 2 : 
pore af M[X - 6zy] : 
eee 
Py i, ay MLy + B2zx] r) 
. : : Bree 2\: . 
p, = a M[z + 8(xy - yx) + Bo(xo + y%)z] = Mz + B(x Py - P) 
(2.10) 


Using equations (2.9) and (2.10), the kinetic energy of the particle in 


terms of the rotating coordinates (x,y,z) and their adjoint momenta 


(P,» Py» P,) is given by 
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‘Hence, in the rotating coordinate system the Hamiltonian is given by, 
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1 2 2 
H = oq [p. + p wy hod 


y + {p, - B(x p 


vf 


+ gk(x® - y*) : pasa) 


Examination of (2.11) shows that H is independent of z. Thus z is a 
cyclic coordinate and Ps is a constant of the motion given by 


Peal taeke: (he D re ey Bela (2.12) 


y 


Using equations (2.7), (2.8) and (2.12) it can be shown that the momenta 


(py Pye p,)s in terms of the fixed coordinates (X,Y,Z) and their adjoint 


mechanical momenta (Py, Ps Po). are given by 


Bes? Py cos 62.4 Py Sin BZ , 


1S) 
iT] 


jes Py SInesZ + Py cos gZ . 


Pet ee(X Py = YP (2,13) 
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Using the Hamiltonian given by equation (2.11), the equations of motion 


in the rotating coordinate system are given as, 
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Equations (2.14) are four first order nonlinear coupled equations in 


Xs Ys Py and p. and they can be used to obtain two second order nonlinear 


y 
coupled differential equations in x and y alone. 
The nonlinear nature of equations (2.14) is removed when the 


transverse momenta conjugate to the rotating coordinates are such that 


<4 1% (2.15) 


It will be shown later, in Chapter 3,that these inequalities are satisfied 
2 


for all particles successfully guided by the structure, provided (ga)~ << 1. 


Under the conditions of equation (2.15), the Hamiltonian can be reduced to 
Z 
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pidsp, (epevery pdur ak(x® = y 
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This Hamiltonian can be used to generate the following linear, but coupled, 


equations of motion 
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ar uw [8P, py - 2Mqk x], 

py = q [-8p, p, + 2Mqky) 

x= q [p, + 8p, ¥] 

y=ylpy-ep,x) (2.17) 


There are two approaches for solving these four linear coupled 
first order differential equations. The first approach is to resort once 
more to Hamiltonian mechanics to generate, through contact transformation, 
a new coordinate system in which the equations of motion are uncoupled. 
The second approach is to put equations (2.17) in a matrix form and then 
to solve the resulting linear first order vector differential equation. 
Both of these approaches are used and described in the two sections which 


follow. 


2.3 Solution of the Equations of Motion by Contact Transformation 


The difficulty encountered in solving equations (2.17) is due to the 
fact that these equations are coupled. However, uncoupled equations of 


motion can be obtained through contact transformation to a new set of 


variables (x,y) and(p, Py). The appropriate generating function ‘34? 


is found to be 
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when the factor y is chosen properly. 
Since G, is time independent, the value of the Hamiltonian remains 


‘unchanged 24) , The canonical variables of the rotating coordinate system 


(Xs Ys Pys Py) and those of the new coordinate system (x, J, ae Py) are 
34) 


related as fol tote: 
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The Hamiltonian of equation (2.16), when expressed in terms of 
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The coupling between the variables (x,y) is represented by the last 

term in equation (2.20) and in order to avoid such coupling the coefficient 


of {x p. - dau should vanish. Thus,the proper choice of y is 
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This Hamiltonian represents uncoupled motion in both the x and y 


directions and the equations of motion are 
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We la 19 
Mom y 4 
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Py =e2gke roa: (2.26) 


Hence,the uncoupled differential equations in x, y are 


y+ is-1jy=0 (2.27) 
Equation (2.25) shows that the sign of s is independent of the sign 


of 8. Taking the above into consideration,it is clear from equations (2.24) 
and (2.27) that the motion of the particle in the (x, y, Pye Py) coordinate 

system is independent of the sign of 8. On the other hand, equation (2.25) 

Shows that the sign of s is dependent on the sign of the charge of the 


particle, q. For a particle with a negative charge.-q.the factor s will 


be related to that of a particle with a positive charge:+q;as, 
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and the equations of motion of the negatively charged particle, obtained 


by substituting Say into equations (2.27), are 


Comparison of Hg and the last two equations with equations (2.24) and 
(2.27) shows that the change of the particle charge from +q to -q 
corresponds to an interchange of x into y and y into x. 

Equations (2.27) represent stable harmonic motion with respect to 


xX and y provided the factor s, which will be referred to as the stability 


factor, satisfies the following condition: 
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The solutions for x and y are direct solutions for equations (2.27). The 


solutions for Py and Py 


xX and y with respect to time, then substituting the resulting expressions 


were derived by differentiating the solutions for 


for x and y into the first two equations of (2.26). In equations (2.29), 
n = 0 for 8 positive, corresponding to a counter-clockwise twist of the 
electrodes, and n = 1 for 8 negative, corresponding to a clockwise twist 


of the electrodes. The other parameters in equations (2.29) are 
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The values of the constants Cy» Cy» oy and Oo depend upon the conditions 
of the particle at injection. Notice that in equations (2.29), (2.30) 


and in all equations to follow a positive square root will be assumed. 


7 i alee 
(ez.$) ‘ ¢ ve oS ga 


aft ASS. $) anoticups ro? 2nottuloe toaith sms has ¢ OF enot ute m 


vot miitiion atti anlistinate?t?o vd tevi tsb 50" x bas J "oY neltalaes 


anotacorns wridi nee a gnc ie st ye aged ,gmtt a? Joeqes Ate ut x 


2 
= : 
r J (0S-4) i bh ow! JeatY ONT ommt % ne ov 

108.5) enofteupe ot .(05.4) to anof LAupe © —_ 


“ait ta sebwd szivigols-teriaoo 6 of tir onbageNTde « svi sieoa 4 1or 0 = f 


: jetw! getwsonts : of} onthnogesitio. .secrppeh 4 407 [s/n bra aac - 
a dye 103.5) anoiseups o! ew esne1sq VSNRy aA! setiortoste: a3 oo 
a 7 
: _ 
SOT lide! iil i 
7 
_ _ t af ~ 4,3 = et 


stanoo arts satel — g® ine rf vicodly sain ola 
re 408.8) ent us née 2é6ns saison — 


vat (H y00" 9Vbure 60fs$20q & wonton ae 


36 
Thus, from equations (2.23) and (2.29) the motion of a charged 
particle in terms of the rotating coordinates (x,y) and their adjoint 


_momenta (Py p.) is given by 


A) 


X = C, cos F + 


(2°31) 


provided that the inequalities expressed by (2.15) are satisfied. These 
equations give the rotating transverse coordinates and their adjoint 
momenta as functions of the time t, where t = 0 at the instant of injection 
of the particle into the structure at Z = 0. 

For the particle trajectories to be completely determined it is 
necessary to find the axial distance Z as a function of the time t and 
the injection conditions of the particle. Equation (2.12) gives the axial 


velocity Z as 
SN (2-32) 


The constant of the motion, Pp,» can be expressed in terms of the injection 


conditions of the particle by using equation (2.13) as, 
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where (Aral) are the coordinates of the particle at injection, and where 
(Py : Py is ) are the components of the momentum of the particle at 

) 0) fe) 
injection. U, is the axial velocity at injection. As shown in Appendix 
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B,equations (2.31), (2.32) and (2.33) give Z as a function of the injection 


conditions and the time t as 
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Integrating equation (2.34) with respect to time gives the axial 


distance traversed by the particle as a function of the time t as, 
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Equations (2.31) together with equation (2.35) completely specify 
the particle trajectory in the rotating (x,y,z) coordinate system, with 
ee 0 for 8 positive and n = 1 for 6 negative. The equations specifying 
the trajectories in the fixed (X,Y,Z) coordinate system are obtained from 
equations (2.31) by using the transform equations (2.7) 

The constants C,> Cos 0 and Gy can be evaluated in terms of the 
injection conditions of the particle. Setting t = 0 in equations (2.7) 
and (2.13), the relationships between the rotating and the fixed coordinate 


systems are 
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Doe Er + B{X, Py -Y) a (2.36) 


Setting t = 0 in equations (2.31) and using the relations (2.36), the 


constants Cys Cys ot and Oo are found to be 
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2.4 Solution of the Equations of Motion by Matrix Techniques 


The second approach to finding an analytical solution for the 
particle trajectories utilizes a matrix formulation of the equations of 
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Equation (2.38) is a first order linear vector differential equation of 


the form, 
AC toN= [Tlntt) (2.40) 
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For the case where the eigenvalues of [T] are distinct, the solution of 


equation (2.40) is known to pe (39) 
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where his Ags Ag and hg are the four eigenvalues of the transform matrix 
[T] and ny: Nos ng and ng are the corresponding four eigenvectors. The 
values of the constants Ay» Aas A, and Ay depend on the components of 
n(t) at t = 0, that is they depend on the injection values of x, y, Py 
and p.. 


y 
The four eigenvalues of [T] are the roots of the equation 


|T - AI] = 0 where I is the 4x4 unit matrix. 


Substituting for [T], the above equation becomes 
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The solution of this equation gives the four eigenvalues as 
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where f, and f. are the same as defined by equations (2.30). In other 


words , 
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where 


Solving equations (2.44) for each of the eigenvalues gives the four 


eigenvectors as, 
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Substituting equations (2.43) and (2.45) into (2.42) gives the solution 


of equation (2.38) as 
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Stability of the solution implies real values for i and f,. eS 
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This stability criterion is the same as that given in (2.28). 
To evaluate the constants Ay» Ay» A3 and Ay consider equations (2.36) 
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Setting t = 0 in equations (2.46) and substituting from equation (2.47) 
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Solving these four algebraic equations, the constants Ay» Ao» A. and Ay 


in terms of the injection conditions are: 
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vail PP, 1 BP2 BP 
pore Seca) Xe 2 eek y Lees 
Bp Bp Bp 
by 7 all eee Pez 
Sa = 2 Bqkw Px ~ J ZFS > Agee Eye (2.48) 


As shown in Appendix C by substituting equations (2.48) into equations (2.46) 
and by rearranging the terms, the equations describing the trajectories 
of the particle in the rotating coordinate system are obtained as functions 


of the injection conditions in the matrix form: 


x Xo - 
y Np i 
= 1 O(t)4 =. Webel ; 
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Py Py. Py (2.49) 


where the matrix Q(t) is given by 


. Loaded ry 
> Lae ah 7 Bi 
mi: ‘eynt {8*, 2) zwof supe oie 


| eatintoatons wit paldboesh caotrgipy Bil’, eates ae i innirenen wh 
anataoaut dx pantatdo ssh caters Sten thts gn Ieroy GA Ai sty tasen via t to 


. 
~_ 


rey? vbstaw 20d pf srotstines natiaoln “— 


1 i270) = 


qd aeyip-2t (30 abagom ald vst 7 


47 


(0S°2) 
i | 
-SA bWZ/ -SA 
4°4 S09 ! Pie cen ree 
rs Z | “* ybyz 
! LS dg 71PWe/ 
i 
2A eS Se oe ee Gk ee ee a See 
L+S/ — 4DWZ/ ' | 
== a} 3S0) | 0 
AyeuLS de ] 
l 
a ey eae ger ee meme ae ee cere Ae ee ce MR Ge cet We | ec I eee re Ra 
| ] 
1 
ybw2/ L+s/ bw 
c F ] c L 1 bw2 2 
ee (7°4 sod - y°4 soo) Ae a SOD 
ee alyuis Fy dg 
JOG, 2a ee RS a > se 
] 
vty soo - 1°4 $09) AOW2 : iy wwe" ak me Nee Le PWe/ 
“dg | 42) urs we ele hk aon is: ode 


: 0 

| 

| 

on L+S/ 

| ULS > 
) gene sore 

| 

| 

| 

} 

| L+s/  ybw2/ 
qty ULS “dg 

! 

jer ie i eH Ke KH Ke 
! 

i] 

; aby soo 


or we 
_ - 7 — 
= 


nie te ee ae rae 


a — = ee eee 


a 


gh 4 a ie 


e i 
ee - -— — — — = = 


ayvare . 
——— at. yes y d 
Te, | Fey 


—s —_-——— — —_ — — <_ — ~ _ - = _ —, - — — -_ — = — —_ — ~~ 


48 


Appendix C shows that the solution given by equations (2.49) and 
(2.50) is obtainable directly from equations (2.31) and (2.37) through the 
expansion of the trignometric terms and the evaluation of the constants 
Cy» Cos 4 and Ot in terms of the injection conditions. One feature of 
the matrix approach is that it circumvents the need for determining a 
Suitable contact transformation. 

The first form of the solution, given by equations (2.31) and (2.37) 
is very useful for investigating the guiding properties of a long uniformly 
twisted electrostatic quadrupole beam transport system. Equations (2.31) 
can be easily used to obtain an expression for the radial displacement 
of the particle and by imposing the condition that the maximum radial 
displacement must be less than or equal to the aperture radius, it is 
possible to derive the acceptance limits of the structure. 

The matrix form of the solution is useful for studying the imaging 
properties of the twisted structure. By inspecting the different elements 
of the transfer matrix, it is possible to investigate the parametric 
changes required to obtain specific transformation properties for a 
structure of finite length. For example, in Chapter 7 it is shown that it 


is possible for Q(t) to be a unit matrix. 
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CHAPTER 3 


GUIDING PROPERTIES 


Se liees Lally, 


The analytical expressions for the particle motion, either 
equations (2.31) or equations (2.49) and (2.50), show that the harmonic 
motion of the particle will be bounded and stable if the two frequencies 


f, and fo are real, that. is 
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These two conditions show that the motion of the particle is bounded and 
Stable if s satisfies the condition, 


fe Z 


p> 
Is| = loggy-l > 1-0 (3.1) 


Substitution for P, from equation (2.33) gives this stability condition 
in terms of the structure parameters and the particle injection conditions 
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‘For particles which prior to injection have been axially accelerated from 
rest through a potential V, the axial momentum is given by Ps = V2qMV. 


0 
For such particles equation (3.2) takes the form 


2,V Xo 0 One Ove 
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3.2 Particle Confinement 


So far it has been shown that particle motion along the uniformly 
twisted quadrupole structure will be stable and oscillatory provided that 
the stability condition as well as the inequalities (2.2) and (2.15) are 
satisfied. However, these conditions are not sufficient to ensure that 
a particle will remain bounded within the aperture of the structure. For 
a particle to be successfully guided, its maximum transverse displacement 
from the axis of the structure must never exceed the aperture radius a. 


In terms of the rotating coordinates 


Substitution using equations (2.31) gives 


ZEAL Z nies : 
Cae C cos F, + 2(-1) I ° 65 sin F. cos F, 


S ees 2 2 
+ ——— C- sin Fo t C5 cos Fy 


oan a Fe a ee 
yd wpyig 2! maser Salne Nd y¥ TE wok hiv 
| wva¥ fp siey (88) nolan jon dn i 


ny q =: - 
(ee) 0 M-Sat = 


_ a " 


a 


: i. a 
giereattnd 9 Grater nc! fom ef pfevey Fer ngerte dina ant 2F WY Oe 
HAT BablyonG. yrote! Foe bon, Bidet? 50 (few soutsonté efomerivap bateteh 
eva (21.5) ton (5:8) eottifeupar! ai? oa Liew Berngherendy cat Bist ot _ 

jaid Seva of doatattiar Yon ow au0hi thngy arent euiviet bottetiee 
‘Gel arias art to svar eye ag: nrigtw babwuod intent 1him i 
rmensitap iB Aerwteanany cant oon 237 haerae eT betas oe ba) ath 

BRS areheags sit Laer > raver dead syyeauite: SAt To hme ony aan 


~ ae: cotanibenags sent aero 

Se3 igh 7 
ee) 2 (ee 
7) _ 


4u2 
ae 


Piteid 0 Ay AE 
ot Sen B+ «3 Fate a . 


- : 


2avip’ (ff :3)-artoideupe gate 4 


“ 


n S : S Pe ee 
- 2(-1) Veet (1% sin F, cos F., ses ea C sin F, 


= Ci[cos F, + Perea Th Fi] 
2rs - | 2 Aw 
ira C,, 3 cos Fo + sin F] 
n Suck |. 
Taek) C1C,, “Es, [sin Fo cos F, oa ean F, cos F.]. 


However, each of the three terms between brackets can only have a magnitude 


equal to or less than unity. Thus, 


or 


S 
Bee Very 7! eo) 


Hence, to ensure successful guiding of the particle by the twisted 
structure, Ic, | and |C.,| should be such that the maximum radial displacement 


satisfies the confinement criterion 


ie ES 
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This inequality sets the limits for the maximum values that |C,| and |C, | 


can have as 
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In Chapter 2, the solutions for the particle motion were based on 
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the assumption that (—)° and (sa) are much less than unity. Consideration 
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Substitution of the maximum values of C, and C. from equations (3.7) into 


equation (3.8) yields: 
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Thus, for any particle whose stable motion has a maximum radial displacement 


less than the aperture radius, 
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Since (ga)? << 1.0 and |s| > 1.0 for stability considerations, equation 
(3.9) shows that the inequality (2.15) is satisfied as long as the 
‘inequalities (3.1) and (3.6) are satisfied. 

Substituting the values of C, | and |C,|. given in equations (2.37), 
into equation (3.6) gives the confinement criterion in terms of the 


structure parameters and the particle injection conditions. Thus, 
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When the injection conditions of a particle are such that this 
confinement criterion is satisfied it will be successfully guided by the 
structure without being intercepted by any of the four twisted electrodes. 
A general discussion of these injection limitations is difficult and only 


some special cases of injection are considered. 
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3.3.1 Parallel Beam 


When a charged particle is injected parallel to the structure axis, 
the Z-axis, its transverse momenta Py and Py are zero and its initial 
0) ) 
transverse displacements are given by X, and Y,: Equations (2.37) under 


these conditions show that, 
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SOT AIAG es Coleen (3.11) 


Consideration of equation (3.10) shows that particles will be 
contained by the uniformly twisted structure when xy and us are within 


the parallelogram - shaped region depicted by figure (3.1) (plotted for 


s = 2.0 and (ga)? = 0.1). This region is bounded by the straight lines 


2 Ee “o, VS evils fy Staal (3412) 


Equation (3.12) shows that the maximum allowable value of Xo? 
when Ie = 0, is equal to the aperture radius irrespective of the value of 
S. Equation (3.12) also shows the maximum acceptable value of Y, When 


: Ss - ] : , 
X. = 0 is equal to ae thus y increases from zero to a maximum 
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value equal to the aperture radius as S increases from 1.0 to infinity. 


It is clear from equation (3.3) that for parallel injection 


2,V 
[s| = (8a) iy 
For a given structure and fixed (ga), increasing s for a particular 
particle corresponds to an increase of its accelerating potential at 
injection or to a decrease of the quadrupole focusing voltage, ‘oe 
Study of figure (3.1) shows that the twisted structure exhibits 

a feature common to classical quadrupole structures. The maximum 
allowable X-displacement for parallel injection in the initially focusing 

X-Z plane is larger than the maximum allowable Y-displacement for 


parallel injection in the initially defocusing Y-Z plane. A particle 
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may be injected in the X-Z plane with an initial displacement as large 
as the aperture radius. At first, the particle will be acted upon by 
‘large electric fields that deflect it towards the Z axis. As it progresses 
along the structure the forces acting on the particle will change and 
tend to deflect it away from the axis. This defocusing action will be 
followed by another focusing action as the forces change again and tend 
to deflect the particle back towards the axis. Thus, the particle will 
experience alternating focusing and defocusing forces as it passes along 
the structure. However, each time the forces are focusing the particle 
is farther away from the axis than it is during the time when it sees 
defocusing forces. Since the fields and hence the forces increase 
linearly with distance from the axis,the effect of those alternating forces 
will be a net focusing action on the particles. 

On the other hand, particles cannot be injected in the Y-Z plane 


with ve close to or equal to the aperture radius. A particle injected 


with VA initially larger than the limit 2 : a will encounter large 


defocusing forces that will deflect it away from the Z-axis so that its 


transverse displacement will exceed the aperture radius and it will be 
See 
S 


lost. For values of ve equal to or less than a, the particle will 
be deflected away from the axis initially but the subsequent focusing 
action will be sufficient to turn the particle back towards the axis 
before it can be lost. As before the alternating focusing and defocusing 
forces will result in a net focusing action. 

Unlike in classical quadrupole structures, the maximum acceptable 


value of Xo for parallel injection along the twisted structure is dependent 


on the initial Nis displacement of the particle, as shown from equation 
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(3.12) and figure (3.1). This dependence is due to the coupling between 


the X and Y components of the particle motion along the twisted structure. 


3.3.2). Point Source 


For particles emerging from a point source at the input plane of 
the structure the initial transverse displacements, Xo and Yo? are zero 


angry. 
a ue 
x. and AP are 
limited to values within the momenta space bounded by the curve 


at Z = 0 and the injection transverse momenta are given by P 


In this case the confinement criterion shows that P 
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where by equation (2.33) Pea P, . Equation (3.13) is plotted (see 
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Appendix D) for several values of s in figure (3.2). The axes have been 
normalized with respect to Ps and the curves are plotted assuming 

O 
(ga)¢ = 0.1. It is shown in Appendix D that the maximum allowable 


transverse momenta in the X and Y directions are: 
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lea P he values of s approaching unity, the value of [Py | | approaches 
oO max 

—z—— while [Py | approaches zero. As s increases the value of 
) Oo max 
[Py | decreases until it finally approaches zero as s approaches 

O max 
infinity. On the other hand as s increases the value of |P, | increases 

oO max 
until s = 1.3295 is reached at which point [Py | reaches a maximum 
o max 

|.© > For Values of s larger than 1'.32955°(P," | decrease: 
Z X 
0 O max 
as Ss increases until [Py | approaches zero as s approaches infinity. 

Oo max 


value of 0.226 |ga P 


Similar to classical quadrupole structures the allowable values of 

[Py | are always greater than the values of |Py | for all values of 
Oo max Oo max 

s- Comparing the alternating focusing and defocusing forces near the 
point of injection provides a qualitative explanation of this difference. 
A particle injected on the axis with Y-directed transverse momentum 
initially experiences a small defocusing force and is accelerated radially 
away from the Z-axis. By the time the forces on the particle change from 
defocusing to focusing, the particle will have been displaced a considerable 
distance from the Z-axis. Here, the focusing forces acting on the particle 
will be stronger than the defocusing forces previously encountered. Thus 
the net force in this initial region of the structure will be focusing. 
On the other hand, a particle injected with X-directed transverse momentum 
first experiences a small focusing force followed by a larger defocusing 
force as it moves away from the Z-axis. In this case the net force in the 
region near the point of injection is defocusing. In other words,a particle 
injected on the axis with X-directed momentum is initially defocused while 
a particle injected with Y-directed momentum is initially focused. Thus, 
it is to be expected that the allowable values of Py are larger than the 
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A comparison of the allowable transverse momenta for point source 
injection in the case of the twisted and the classical quadrupole structures 
will be presented in Chapter 6. 

Considering the cases of parallel injection and point source 
injection for the twisted structure, it is found that large values of s 
are desirable for parallel injection while they result in a reduction of 
the momentum acceptance for point source injection. A compromise value 


of s = 2 seems appropriate for both cases. 


3.3.3 Injection on the X-Axis with X-Directed Momentum 


A particle injected on the X-axis with X-directed transverse 


momentum has Y. = * = 0. For such a particle p, = ae [S| = (Ba) ol 
and equation (3.10) gives the confinement criterion for this case as, 
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Taking into consideration that |s| > 1.0, the confinement criterion finally 
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becomes: p 


X 
Nase! - a\/s“ - s}¢ 9 
is Xo 
oy eee AG | (3.16) 


2 On = 
{a(s - 1) Shy tafs? - 1} 


pees for such a particle to be successfully guided by the structure, 


XS and a should be represented by a point within the area of intersection 


of the two-hyperbolae represented by equation (3.16), as shown in figure 
(3.3). The axes in figure (3.3) have been normalized, Xo with respect to 


a and Py with respect to Ps . Three cases were plotted in figure (3.3) 
6) 
for values of s = 1.5, 2.0 and 3.0. For all values of s the intersections 
X 
of the hyperbolae occur on the GS axis at = = ite semihis CeSuli 1s in 


agreement with the results obtained for the special SSE of parallel 
X 

injection since figure (3.3) clearly shows that for p= 0, together with 
Z 


Y, = Py = 0, the maximum allowable IX is equal to av Figure (3.3) also 
0 
Shows that as s increases, the allowable transverse momentum decreases. 


Note that all values of s considered are larger than s = 1.3295 at which 


|P. | reaches its maximum value. The values shown in figure (3.3) for 
PO max 
am » at X. = 0, are the same as those predicted by equation (3.14) 


o Max ; 
for the values of s considered. 


3.3.4 Injection on the Y-Axis with Y-Directed Momentum 


A particle injected right on the Y-axis with Y-directed transverse 
momentum has ie = 0 and Re = 0. As was the case for the other special 
cases, p, = oh, and s = (5a) "Ip. The confinement criterion for this 
special case of injection, obtained by substituting in equation (3.10), is 
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This equation shows that such a particle is successfully guided 
by the structure if re and (eo) are represented by a point within the 
area of the intersection of the two parabolas represented by equation (3.17). 
Several such parabolas are plotted in figure (3.4). The axes in figure 
(3.4) have been normalized, ee with respect to a and Py with respect to 
the Three cases were plotted in figure (3.4) for bhiee Gh (Ser ld 
2.0 and 3.0. For the three values of s the intersection of the parabolas 
occurs on the ats but the point of intersection depends on the value 


of s. As expected from the special case of parallel injection, the 


maximum allowable values of Ya when Py = 0 are equal to ye : and 
fo) 
decrease as s decreases. Also as expected from the special case of point 


source injection, the acceptable values of Py for a particle injected at 
0 
a point up with Xo = P, = 0 decrease as S increases. 


Xo 


3.3.5 Point Source at a Distance d From Input Plane 


Consider a point source a distance d from the input plane of the 
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twisted structure. A particle emerging from this source will have its 


transverse coordinates at the input plane, Z = 0, related to its momenta 


’ by 
oe 
Xo eat > 
0 
ay 
i Sipe tae . (25033 
Ly 


Thus (X, ey - ihe Ea = 0 and es i for such particles. Substituting 


equations (3.18) into equation (3.10) gives the confinement criterion for 


this special case as 


p Py p, @ 
fea Ow SOe sy Sis Foe 
a P Ba Po 2,¢ 2 
@) ) Ly 
Py p p, 2 

d 2 : 

ty ie 2 - Be + 1) oye <1.0 (3.19) 
) Ly B de 


Equation (3.19) gives the values of the allowable divergence that 
a beam, emitted from the point source distant d from the input plane, can 


have so that all particles will be guided. The acceptance area in the 


X vd 
transverse momentum es By hase plane is bounded by an ellipsoidal curve 
O O 2 
WER eri re lancciuse of che icrosss tenis: PAB Es 4 Phewee pena and 
Xo Le No Yo NG Yo 
é y é that will appear in the expression of equation (3.19) if both 
0 fo) 


sides of the inequality are squared. The degree of tilt of the acceptance 
curve is dependent on the value of < considered as well as the structure 


parameters. 
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For ¢ - 0.0, equation (3.19) will reduce to equation (3.13) 
obtained for point source injection at the injection plane, Z = 0.0, and 
in this case the acceptance curve will not be tilted. 

Figure (3.5) shows the acceptance limits in the transverse normalized 
momentum phase plane. The axes are normalized with respect to Ps . The 
acceptance curves are for s = 1.25 and s = 2.0 for a AES ee 


twisted structure, 8 positive, having (ga)? = 0.1. The point source is 


at a distance d = 5a from the injection plane. 
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CHAPTER 4 
AN EVALUATION OF THE ANALYTIC SOLUTION 
4.1 Introduction 


The analytic solution developed in Chapter 2 is based on the 


p p 
assumptions that (ga)° << 1 and (X)4, (vyé << 1. It was shown in 


p p 
Chapter 3 that these assumptions aca justified for any particle whose 
maximum radial displacement from the structure axis is less than the 
aperture radius. Therefore, the approximate analytic solution is not 
restricted to small displacements from the Z axis. Further, it was 
proved in Chapter 3 that a particle will be successfully guided if its 
transverse coordinates and its momenta at injection satisfy the acceptance 
or confinement criterion given by equation (3.10). The first objective of 
this chapter is to establish the accuracy of the approximate analytic 
solution and in particular to test the validity of equation (3.10) in 
predicting the acceptance limits of the twisted structure. This objective 
is achieved through comparison of the trajectories computed using the 
analytic solution with the actual trajectories of the particles. The 
latter trajectories are computed through numerical integration of equations 
(2.6) which are the original equations of motion. These equations of 
motion as well as the approximate analytic solution are based on the 


approximate expression of the potential distribution given by equations 
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(2.1a) and (2.1b). The second objective of this chapter is to present 

a mode theory that is based on the general examination of the nature of 

- the trajectories and which describes the projections of the particle 
trajectories on the fixed X-Z and Y-Z planes in terms of four spatial 
modes. This mode theory explains the apparent amplitude modulation in 
some of the particle trajectories as well as the presence of small ripple 
components superimposed on longer wavelength fundamental components. 
Using the expressions provided by this mode theory,for X and Y parametric 
in Z, it is possible to compute and plot average approximate particle 
trajectories instead of using the long expressions of equations (2.31) 
and, (2335)).. 

A number of trajectories were computed and plotted using the analytic 
expressions of equations (2.31) and (2.35). The trajectories were plotted 
for a large number of injection coordinates and momenta satisfying equation 
(3.10). These trajectories will be referred to as the analytic trajectories. 
The first program listing in Appendix E shows the steps of computing these 
analytic trajectories. A second set of trajectories was then computed 
through the numerical integration of equations (2.6), the original non- 
linear coupled equations of motion, using the Runge-Kutta method of 
numerical ntegRet ibn ee These trajectories were plotted for the 
same injection conditions considered for the analytic trajectories. The 
steps of the computations of these trajectories also are shown in the first 
program listing of Appendix E. They will be referred to as the Runge-Kutta 
trajectories and will be taken to represent the actual trajectories. 

The numerical solution provided by the Runge-Kutta method is a step- 


by-step solution. The accuracy of such a step-by-step solution of 
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differential equations is often difficult to se The Runge- 
Kutta method has no error checking equations. However, the truncation 


9 where h is the time interval used for 


error is near the order of h 
computing X, Y, Z. Improvement on the accuracy may be achieved by using 
smaller time intervals. However, a decrease of the interval size, h, 
adds to the amount of computing time and increases the possible round-off- 
error. For choosing a suitable value of h so that both the truncation 
error and the round-off-error will be negligible, a set of Runge-Kutta 
trajectories were computed and their projections on the X-Z and the Y-Z 
planes were plotted using different values of h for the same particle. 
The value of h was varied from sau up to aut where L is the 
periodic length of the twisted structure and Us is°the initial axial 
velocity of the particle at injection. The eeanlitd cons were performed for 
a structure whose length is 24 meters, that is 800 L approximately. It 
was found that there is excellent agreement between the trajectories, 
computed and plotted for this long structure, corresponding to values of h 
from SOOT UP to a The consistency of the trajectories for this 
wide range 8f variation of h indicates that the value of h should be 
chosen between soo UT and pe It was decided to use a time interval 
h equal to sep in’ the computations of the Runge-Kutta trajectories. 
The same value of°h is used in computing points of the analytic trajector- 
ies. 

For the purpose of assessing the accuracy of the approximate 
analytic solution the analytic and Runge-Kutta trajectories are compared 
for a structure having an aperture radius, a = 1.5 cm and a periodic 


aaa 


length L = 29.68 cm which corresponds to (ga)” = c = 0.1, hence 
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13 
Satisfying equation (2.2). As shown in Chapter 3, a stability factor 
s = 2.0 will allow good use of the aperture for parallel injection and 
‘it also yields a relatively symmetric transverse momenta acceptance 
region for point source injection. It was also shown that for both 
parallel and point source injection s = (sa)“|]. Therefore, s can be 
set equal to 2.0, for parallel and point Bt ce Ae isis by setting 
the quadrupole voltage, oe equal to 1.5 KV for particles preaccelerated 
from rest through a potential V = 30 KV. It is evident that for 
particles with general injection conditions the value of s must be 


calculated for each particle using 
2V le eee 
Sha Ga) ain] Wise (6a). aoe preset (4.1) 
6) 


The discussion and the listing of the computer program used to compute 
and plot the analytic and Runge-Kutta trajectories are given in Appendix 
as 

Since only electrostatic fields are involved and all the particles 
are presumed to have been preaccelerated from rest through the same 
potential prior to injection, the trajectory of a particle along the 
structure is independent of the particle's charge to mass ratio. For 
computation of the trajectories the particles were assumed to have a 
charge to mass ratio equal to 30 coulomb/kilogram, which is typical of 
particles used in micrometeoroid simulation studies. In order to compute 
the analytic trajectories the values of the injection conditions Xo? 5? 


y and Ps were used to evaluate the constants of motion s, 
fe) fo) 0) 


C1» Cos Oy» A and P, using equations (4.1), (2.33) and (2.37). . Then 
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by increasing the time t in steps, each equal to h seconds, and using 
equations (2.31), (2.34) and (2.35), the values of (x, y, Pra Py) and 
‘the corresponding values of the axial velocity Z and the axial distance 
Z were evaluated for each value of t. The analytic trajectories were 
then plotted in terms of the fixed (X, Y, Z) coordinates using equations 
(2.7). On the other hand, the Runge-Kutta trajectories computed using 
equations (2.6) are already in terms of the fixed, (X, Y, Z) coordinates. 
The computations of both sets of trajectories were carried out 
using an IBM 360 digital computer and the projections of the trajectories 
on the X-Z plane, the Y-Z plane and the X-Y plane as well as Z versus Z 
were generated directly by the computer using a Calcomp Plotter model 


770/663. The various trajectories are presented in figures (4.1) to 


(4.19). 


4.2 Comparison of the Analytic and the Runge-Kutta Trajectories 


The trajectories discussed below have been computed and plotted for 
a large number of injection conditions. Some of these injection conditions 


were general, that means particles having nonzero values of Xoo Year 


X 


0° 
0 


and P The restof the injection conditions considered were for the 


y° 
0 
special cases of parallel and point source injection at Z = 0. 


Figures (4.1) to (4.7) show the X-Z, Y-Z and X-Y projections of the 
trajectories for cases where the injection conditions are well within the 
limits set by equation (3.10). Figures (4.1) and (4.2) pre for the 
general injection conditions (X, = 0.1 cm, y; = 072 CMs Dae 0.02 and 


Py Px 
a " a Og. Ope Ove 
ae 0.02) and (X, = 0.3 cm, Y, = 0.2 cm, P 0.01 and P 0.01) 
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respectively. Figures (4.3), (4.4) and (4.5) are for parallel injection 
at the points (X, = 0.45 cm, le = 0 ), (X, = (0.337 cm, Mes = 0.079 cm) 


and (Xo = 0.225 il, ee 0.159 cm) respectively. For point source 


injection figures (4.6) and (457) show the trajectories for paneicles with 
X y 

injection transverse momenta (—— = 0.016, a= 0.011) and (5 = 0; 007. 

Pp Z 

Y 0 O 


= 0.022) respectively. For see of the particles, there is nearly 
exact agreement between the analytic and the numerical trajectories. For 
all seven cases, figures (4.1) to (4.7) show that the X-Z and the Y-Z 
projections of the trajectories feature a dominant oscillation superimposed 
upon which there is a small ripple. The amplitudes and positions of the 
dominant oscillations and the ripples for each of these seven examples are 
found to be exactly the same in both the analytical and the Runge-Kutta 
trajectories. Moreover, the (X-Y) projections of the trajectories in 
figures (4.1) to (4.7) show that the maximum displacement is always less 
than the quadrupole aperture, a = 1.5 cm. 

Figures (4.8) to (4.11) depict trajectories for particles with 
injection conditions chosen such that, as expected from studying equation 
(3.10), their maximum radial displacements are equal to or greater than 
half the aperture radius. Figures (4.8) and (4.9) are for parallel 
injection at (X, = 0.675 cm, ue = 0.159 cm) and (X, =rQt. J ie = 0.53 cm) 
respectively. Figures (4510) and (4.1) are for point source injection 


: Xo S ome Cue 
with transverse momenta ‘pigs ONE pane 0.051) and(5— = 0.047, 
Z 


P 
Y O O 0 
> = 0.012) respectively. All four cases show substantial agreement 


fe) 
between the approximate analytic and the Runge-Kutta trajectories. The 


agreement is nearly exact up to a length equal to approximately 6L, or 


- ie - oe! . 
a “woliostnt | faltevay 107 ste (2 Boren (nio} oh ‘Se ”" said baits 
(ms. 20.0°= 1 fms, 06.0 = A) tO gai io | 

syne 2atog. 1? .¢lav broaqesY ine “— ahs en, 
ashe — 70% 25!narcele A eft wars ie bre silat 
« 1048: 79 we (Creo r (370.0 ae ti 
aie wient , 2a! eae aii? to ie 10% . yfewltapagern. 190.0 «= 
«St .caradtssints Te: own ont bop 2NGTane ae Aare ridin 398 
coy oe ee Red aol aig owt C88) jon TP) eatin? , oso never Ti 
Gaznqaivsque hotdelttow snestt® s ata? zeisosssiew ant to = 
aie Vo-enbtsheog bat eotweliqm wl? = .aiqgm ffana & ct. ovat? totdw n 
Eye eslquane naves stad! Fe dyes Ty tiga ad Ons enptisi {idee snstiba 
Sp4t-sendlt ane ban toothy ene sh? 1700 #) anes aft ulsaexe ad of bruo¥ 
a) egtiineiet? sit to. tnoriosiowg (4-0) ant ,vavusrall. . 2s) 0329 s 


» 


erat augdly. Zt ansmeasiqet! naujxen gtr Iaheee t).ET 02 (1.0) eowe 


- 


7 
. 


pet 
a . 


2 


« 


i) 2 | = & pa vilitegs ‘SToquibaup. odd ne iv 

Aste aatatacey qi exdonsosiars ToS UTE) of (8,8) ever | 

netfiupe gatvbute aes? patoecee ci . Jens (Cue meenes shotsihnes nolsostar 
Alay *iastp-9 6) Tae maieelgeth (shhey dante vio? {08 h 
forfeveg 36% ow (P8) ona (8.5) steph” curber ras i vi 


stat 
noreaapit 054 Jaton tet ave (fi A) fine a eon alanis r yee 
90 = = (130.0 = At.0 = | stemmon senevenest Aah 

| aon - q 

Sheaves (sisnedcowe wore. eae suet £ ia sts 


7 7 


: . oft egtretasta9 KituNeagius Sy tire ich sd 


7 “” oi sclala OF Teupa MpAal os OF: Qi 


lea CE,0 = 1. ) = awe bine ‘im Ge). O = |) fie tai = oh) oe) no 


8] 


0 
"(wo 6SL°0 = A WO G/9°0 = 


‘WBOT LA TBNY 


B1LLNY-39NN¥ 


°y) USYM UOLJIACUT LBL {eied wos (L*y) aunBly ul sy 


o°O 


WOT LANE Tet 


“31 
BLLAYM-SONNY TS"! 


(8°~) aunbL4 


ae 
Tr 


Soot 
BLLAY-SONNG T3t- 


To 


«ti 
0 
ru eter @F 
a 
io 
* tl 
Meeloy 
antzs 
aint. tsi 
Tewa 
q 
at 
“Ut } 
up 
ull a? ah 


82 


SST! J WN 


“s =" BLLM-3ONNE 


‘“ 


X) UaYM UOLZDaECUT LaLLeueg “OF (Lt) BUNBLy uL sy (6°) aunbL4 


se nunca WIT LA TONE 


is 


UBLLNY-3SONNY 79 °t- BLINY-39NNY 


+ 


b x i = 
* — hugs ol P e ee | 7 
7 v4 y 
a 
| 
Pf 
7 
7 ae = A) ee an) am, 
a cee “pene” ees (ies | 
- 
or 
: — 


a 


cs 0.9) oa 


a : no : 


Ame 20 < 3 0 > ki nets nottoatal fellow onl (Phat vi 


_ 


83 


BLLINYA-3ONNY 


0 
i 
= =-d) U9YM UOLZDSCUT BdUNOS YULOd MOF (L*H) BUNBLY uL sy (OL’p) aunbLy 
X 
d 


(W) Z 
o'e Sat 0'O 
——}—+$ + ++ 4 
WOT LATYNY 


le 


BLLAY-39NNY 73"! 


ae 
fre 


san 


‘9 
—t 


(W) Z 
sh O"e Sit 0°O 
$+ + —-$ + "+ 


TOT LA TBNG TS?" t- 


HLLAY-39NNY 


a 
a 


* 
~*~ ~ j eae 
a. f a P| 
; A 
~ — | aes = 
7 j 
— 1 
7 Ta 4 
@ 
a 7] \ 
ae 
= * 
P| 
oe 


a, 
es 


)0 


u 
c+ 
Lhe CS) Us eyreq a 6 O'S) 
A 
7 


® 
b 


(0) O 
Z 
AZPOt0p="— > T2p020 = —4) yaym uolqoafuy aounos qulog 4os (L*p) adnBLy UL sy (LL‘p) eunbL4 


WWST LA ING WIT LATING 73" ‘WOT LA NY 


= 


SLLAY-3ZIONNG BLINY-SONNY TS '- BLINY-SONNY TS" '- 


i 


6 


ae 


a ES — 


od 


* rom : ® = : y ' , 2 - 7 
[S08 “a 10,0 © gia ners sets arms recoo--e® 0). Piaget ab ah Rie 


o oe oe : —s : 
a a ' a 7 - 
+e _ 7 7 a 
- = 7 
> 


85 


about 1.8 meters. At this point a difference appears between the analytical 
and the Runge-Kutta trajectories in the form of an apparent phase shift 
between the trajectories. This small apparent phase shift arises because 
the numerical trajectories cross the Z axis before the corresponding 
analytic trajectories. A study of the previous four cases shows that the 
phase shift increases as the injection conditions approach the limiting 
values. This is especially evident in the positions of the ripple 
components. Note that the X-Y projections of the analytic and the Runge- 
Kutta trajectories, shown in figures (4.8) to (4.11) show that, as 
predicted by equation (3.10), the maximum radial displacements of the 
particles are less than the apterture radius. 

Figures (4.12) and (4.17) represent particles for which the 
injection conditions are exactly equal to the limits set by equation 
(3.10). Figures (4.12), (4.13), and (4.14) show the trajectory projections 
for parallel injection at the limiting acceptance points (xX, =a =e omen: 5 
i FO) (xX, = 1.125 cm, ve = 0.265 cm) and (xX, = 0.187 cm, ia = 0.928 
cm). Figures (4.15), (4.16) and ale) are for sok source eo 


aa Medeor Xo _ ce om 
with the limiting transverse momenta (ome 0.0599, fo = Or ire Oe0525 
Z 


10) 0 be 


“2 
as) 
aw) 


X if 
2 =) 0.038) and == 0, ~~ = 0.079). Even in these limiting cases the 


Zo Lo Lo 


essential features of the motions are quite accurately predicted by the 


“3 


analytic solution. However, examination of the trajectory projections 
again shows that the actual (Runge-Kutta) trajectory projections cross 
the Z axis somewhat before the corresponding analytic trajectory projections. 
In spite of the differences between the analytic and the actual trajectories 


for these limiting injection conditions, it should be noted that the X-Y 
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projections of the trajectories show that the maximum radial displacement 
is equal to or very slightly less than the aperture radius. In other 
words equation (3.10) has quite successfully predicted the acceptance 
limits of the twisted structure for these cases. 

The study of all the above trajectories, both analytic and Runge- 
Kutta, shows that for all injection conditions accepted by the twisted 
structure the projections of the trajectories on the fixed X-Z plane and 
the fixed Y-Z plane feature a dominant oscillation superimposed upon which 
there is a shorter wavelength ripple. The study also shows that the 
relative amplitudes and positions of the ripple component with respect to 
the dominant oscillation depend on the particle's injection conditions. 

For particles injected with small transverse coordinates (X_, Y_) 


0) O 


and/or small transverse momenta Pye Py) that is for particles whose 
injection conditions are well within the limits set by equation (3.10), 
the comparison shows excellent agreement between the analytic and the 
Runge-Kutta trajectories. For such injection conditions, the amplitudes 
and positions of either the dominant oscillations or the ripples are the 
Same in the analytic and the Runge-Kutta trajectories. Such excellent 
agreement indicates the degree of accuracy of the analytic solution and 
shows that the approximate assumptions of equation (2.15) are well 
justified for injection conditions well within the limits set by equation 
(310) « 

The comparison of the analytic and the Runge-Kutta trajectories 
for particles with injection conditions larger than those discussed above, 
but still within the limits set by equation (3.10), shows that there is 


still good agreement between the analytic and the Runge-Kutta trajectories. 
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Any differences between the two solutions are not evident for the first 
few periodic twists. Then a slight apparent phase shift starts to appear 
between the projections, on the X-Z and the Y-Z planes, of the analytic 
and the Runge-Kutta trajectories. The magnitude of this apparent phase 
shift increases as the injection conditions of the particle are increased 
towards the limiting values set by equation (3.10). This apparent phase 
shift arises because each time the Runge-Kutta trajectory projections 
cross the Z axis they do so earlier than the corresponding analytic 
trajectory projections. It would appear that in the actual structure the 
forces on the particle cause it to be returned towards the axis sooner 
than indicated by the analytic solution. In other words, the focusing 
action on the particles displaced from the Z-axis is actually somewhat 
stronger than is predicted by the analytic solution. 

For injection conditions equal to the limits set by equation (3.10), 
there is still agreement in the general nature of the analytic and the 
Runge-Kutta trajectories. However, the apparent phase shift is larger 
than before indicating that the analytic solution fails to account fully 
for the strong focusing action of the actual structure. 

It is important to note that examination of all the X-Y projections 
of the analytic and Runge-Kutta trajectories indicates that the maximum 
displacements of particles whose injection conditions satisfy equation 
(3.10) are equal to or less than the aperture radius, a. This result 
strongly supports the validity of using equation (3.10) in predicting the 


allowable injection conditions of the particles. 
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4.3 The Axial Velocity Z 


Equation (2.34) gives the axial velocity of the particle at any 
position along the twisted structure as a function of its axial velocity 
at injection, sce its injection transverse coordinates and momenta, 
and the time, t, elapsed after injection. Equation (2.34) shows that Z 
varies about ie as the particle moves along the structure. It is also 
Shown in Appendix B that, for particles whose maximum radial displacements 
from the Z axis are equal to or less than the aperture radius a, the 


maximum deviation of Z from U, is given by 


O 
2 2 
, [1 = tea ; ; [1 a 1 
Hp) y= ee, 


2 


0.95 i 1.05 for s = 2 and (ga)~ = 0.1 
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which indicates that the axial velocity of any particle successfully 
guided by the structure is always approximately equal to the particle's 
axial velocity at Pie ae at all points along the twisted 
Suructure. 

For the large number of injection conditions considered in the 
comparison of the actual and analytical trajectories, the axial velocity 
Z was plotted versus the axial distance Z. In all cases the deviation 
between Z and U, was found to increase as the injection conditions of 


0 
the particle approached the acceptance limits set by equation (3.10). 
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It was also observed that for all actual numerical and analytic 
trajectories investigated the deviation was less than 3%, as shown in 
figures (4.18) a, b, c, d, e and f, for the particles with limiting 
injection conditions whose trajectories are shown in figures (4.12), 
(4.13), (4.14), (4.15), (4.16) and (4.17) respectively. Thus, without 
Significantly affecting the accuracy of the solution, it may be assumed 


that, for particles confined by the twisted structure 
(22) 


Figure (4.19) shows the trajectory projections for a particle with 
the injection conditions, Xo = 0.9 cm, Le = 0,.4cm and os = ae = 0. 
The (X-Z) and the (Y-Z) trajectory projections at the top of figure (4.19) 
were computed using equations (2.31) and the axial distance Z was computed 


as 
t 

Z = { Z dt as given by equation (2.35). 
0 


The X-Z and the Y-Z trajectory projections at the bottom of figure 
(4.19) were computed using equations (2.31) and the approximate value of 


the axial distance Z was computed as 


Comparison of the two trajectories verifies that the accuracy of the 
analytic solution is not significantly affected when using the above 


approximation for Z instead of the value given by equation (2.35). 
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4.4. Spatial Modes 


The trajectories investigated exhibit three dominant features. 
Each consists of a long wavelength motion upon which is superimposed a 
small ripple of much shorter wavelength. Further, some of the 
trajectories display a slow amplitude modulation. All of these 
characteristics can be explained by expressing the trajectories in terms 
of cosine - like and sine - like modes of four fundamental spatial 
frequencies. It is shown in Appendix F that the trajectories given by 
equations (2.31) can be transformed to the fixed (X,Y,Z) coordinate frame 
by a simple rotation of coordinates, through an angle gZ, using equations 


(2.7). Thus, im terms of X, Y, Z and t the particle trajectories are 


given by: 
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| 
+5 Ci{1 - [Sql 3sin(t,t + 6Z + a4) 


Sap tic 
Ss - 1! 


] 
+ x Cot] -[ tcos (fat + eZ + On): i) 


Using equation (4.2), the expressions for fit and fot for particles 


confined by the structure become: 


ee eee aay alegre Cy, 
] ] U5 M PS S 
6) O 
OSS1 Tit toe. leqk 1/2 M Secale. 
fot = f, 0, Z nae Pp SB relia! (4.4) 
O (6) 


Substituting equations (4.4) into equations (4.3), the trajectories become 


linear combinations of modes at the following spatial frequencies: 
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Consider the trajectory shown in figure (4.12). For this case 


C, GG WI Relea ie Cy = 0.0, a1 = 0.0 and Oo = 0.0. Hence from equations 


ey (0.4) and (4.5) % 


. 
Pe 


(ys 4 te ¢ ay Vintet tr oo 
a uy 5 


ge + S42 ion Ep , Sa 
tests 7) ae 


salbieea 407 9,7 bio 4,4 voT 2notactiqns ott (80) nekseugs igntav a 


rere 
: gq hte, 25 A, Mey + fils * ee 
- bt iat 
. feb) eV tS) a4 Mer = ee, = ttn P| cer Hef ' 
e aoe llene 7 


smeped estrvotzeters edt .(E-)) eentanupe odat (PLAY anotseupe has oe 
a c 
-sofsioupe?t Ishlege eniwalict of? 26 eebom Yo enottentdmoo 7 ey 


m4 

a a 
wre sta # S\Sa-= ay) = yo ( 
AME ort a Reng * 
ah *Mieah 2 S\(Sa 4 ay) = ge _ 
: . 
| At 

of SMe) a. Sh aT = ye 


_ (STR) oruptt mt woe \iodnatens 


101 


XK = 5 atl + [Sy]! /*icos wz 
+ 5 atl - oS fi 7]'!*}cos WZ : 
Y= - atl + (Sy)! 3sin wiz 
- 5 atl 2 fe sa] /43sin WZ . (4.6) 


Clearly, W1 < W and the X-Z and Y-Z projections of the trajectory consist 


of a dominant oscillation at Wy with corresponding wavelength, 


,l/2 


Me (ee CF=n 
ee me es 


upon which is superimposed a ripple at W3 > with corresponding wavelength, 


,l/2 


Ff gl 
i ee 


r 
In figure (4.12), s = 2, thus ae 4.45L and 4, = 0.45L. Similarly if 
injection conditions were such as to excite only the Wo and Wy modes, as 
in figure (4.9), the trajectory would consist of a dominant oscillation 
at wo with dy = 3.42L (s = 2) and a ripple at wg with hy = 0.586L (s = 2). 
In general, as in the other cases investigated in Section 4.2, all four 
modes are excited and the interpretation in terms of the spatial frequencies 
becomes more difficult. However, it is found that a general trajectory 
will feature a dominant oscillation with an apparent wavelength between 


dy and ho (approximately four times the periodic length of the twisted 
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structure for s = 2). Superimposed upon this oscillation will appear a 
small ripple with an apparent wavelength between h3 and da (approximately 
equal to half the periodic length for s = 2). 

When all modes are present a slow beat will be evident between the 
two large amplitude modes (wy » Wo). As a result the dominant oscillation 
will appear to be amplitude modulated. An excellent example of this 
modulating effect is seen in figure (4.11). In this case equations (2.37) 
show that C, = 0.561 cm, C, = 0.448 cm, a 


= 258° and a, = 169.6°. Thus, 


] 2 ] 2 
from (4.3), (4.4) and (4.5) the expressions for X and Y become 


(wy - 0) (wy + wy) 
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0.03 cOS(woZ + 10.37) 40. 049 cOS(w,Z = 10.37") Lom). 


The major portion of the motion consists of an oscillation at a frequency 
lying midway between Wy and Wo which is amplitude modulated at a rate 


equal to Hu ~ wy). The corresponding modulating wavelength is given by: 
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CHAPTER 5 


EFFECT OF EARTH'S GRAVITY ON PARTICLE TRAJECTORIES 


As mentioned previously the main objective of the analysis of the 
uniformly twisted electrostatic quadrupole is to provide a structure 
capable of guiding heavy charged particles, in particular charged micron- 
Sized particles. Since the charge to mass ratios of such particles are 
very low, the earth's gravitational force could noticeably affect the 
motion of these particles along the twisted guiding channel. The 
objective of this chapter is to determine to what extent gravity modifies 
the particle trajectories and to find values of the focusing quadrupole 


potential, vs beyond which the gravitational effect can be neglected. 


5.1 Particle Motion 


Taking the horizontal axis of the twisted structure, the Z-axis, 
as the zero gravitational potential level, the gravitational potential, 
Vg at any point is 


V_ = MGY = MG(x sin gz + y Sin Bz) (51) 


where the (X,Y,Z) and the (x,y,z) coordinate systems are the same as 
defined in Chapter 2, where M is the mass of the particle and where G is 


the gravitational acceleration. 
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The earth's gravitational forces are accounted for by adding the 
potential Ue to the Hamiltonian of motion given by equation (2.16). The 


new Hamiltonian in terms of the rotating coordinates is, 
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+igk(x” - y~) + MG[x sin. 6z + y cos gz] (5.2) 


where Py» P and p, are as defined by equations (2.10) and it is still 


y 


assumed that ae and ee << 1. The equations of motion generated 
Z z 
from (5.2) are linear, inhomogeneous and coupled. By resorting to 


contact transformation using the generating function 
e seeltes = 
Fee ee hy = Dame Peete (553) 
a new coordinate system (x, y; Z; Pye Pye p,) is generated in which the 


equations of motion are uncoupled provided the constant p is selected 


properly. The new coordinates are related to the old ones as, 
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The following assumptions are now introduced: 
Teel eiseassumed ‘that P, = Pz. This follows from (5.4) AF ‘the 


first order quantity, B(x p. - y Py)» is neglected. 


: p 
2. It was shown in Chapter 4 that the axial velocity Z = ‘a 
is essentially constant and equal to its value at injection U5 2 Palithals 


fe) 
within the same degree of approximation to assume that 


Dee Peep SMU (5.7) 
Taking (5.7) into consideration, equation (5.6) shows that the 
term (x Py - y p,) is the only coupling term between the coordinates 
Xe Vs P, and Pye and this term vanishes if p is selected as, 


p = Ps = Py (5.8) 


Substituting (5.7) and (5.8) into (5.6), the Hamiltonian reduces to, 
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= = Oo- = 
MGLy cos 6z + ZaMk Py sin 6z] (5.9) 
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= 2k at ; 

xX + = (s, + 1)x [2s,, + 1]G sin B a tes 

Gaeneds.(os aii y eee emilarcaen an ee (5.10) 
M ) y 0 Ly : 


and two similar equations for Py and Py where So 1S a stability factor. 
It is evident from equations (5.10) that the motion is stable if 


Iso > 1 where 


a Sor | ace (5.11) 
In equations (5.10) the approximation z = Z z U, t is introduced. It 
is evident from equations (5.4), (5.7) and the assumption that ais 
essentially constant, that this approximation is justified. 

The solutions of the inhomogeneous differential equations (5.10) 
consist of the homogeneous solutions, exactly analogous to the solutions 
given by equation (2.29) for the particle trajectories in the absence of 


gravitational forces, and the particular solutions given by 
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When transformed to the fixed (X,Y,Z) coordinate system, the particular 
solutions resulting from the gravitational forces are, 
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Figure (5.1) The Mean Orbit of Particle Motion when the Effect 


of Gravitational Forces is Taken into Consideration. 


The projection of the particular solutions of equations (5.13) onto 

the fixed X-Y plane is plotted in figure (5.1). This figure indicates 
that the mean orbit of the particle trajectories, which was the 
horizontal Z axis when gravitational forces were neglected, has been 
modified to a helical orbit. The axis of this helical orbit is parallel 
to the Z axis at a distance rae the negative side of the Y-axis 

and the radius of this helical orbit is mK As expected, equations 
(5.13) and figure (5.1) show that the effect of the gravitational forces 


becomes noticeable as H decreases. 
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For the effect of gravitational forces to be neglected, the 
distance between the axis of the helical orbit and the Z axis should be 


much less than the aperture radius. This condition implies choosing the 


focusing potential to satisfy 


MG 
4S 0 2qk ES Gl 3 
that is 
pe (5.14) 
(HV 
(@) 


As stated before the value of i for micron sized charged particles 
may be as low as 0.01 coulomb/kilogram. For G = 9.8] m/sec? and a = 0.015 
meter equation (5.14) becomes 
me 
.294 7 So 
ee ees (5.15) 
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For s,= 2.0 the inequality (5.15) implies that 


Mi ee 58.6 Volts for i = 0.01 coulomb/kilogram. 


The parameters of the lens structure which is considered in 
Chapters 7 and 8 are 
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Equation (5.14) shows that for these parameters gravity becomes negligible 


provided that the charge to mass ratios of the particles satisfy 
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CHAPTER 6 


A COMPARISON OF THE PROPERTIES OF THE UNIFORMLY TWISTED ELECTROSTATIC 
QUADRUPOLE STRUCTURE AND THE CLASSICAL ELECTROSTATIC QUADRUPOLE CHANNEL 


6.1 Introduction 


The theoretical and experimental investigations of the guiding 
properties of the uniformly twisted electrostatic quadrupole show that 
it can successfully guide charged particles provided that equations 
(3.3) and (3.10) are satisfied. It is of interest to compare the 
guiding properties of the twisted structure with those of a classical 
electrostatic quadrupole channel. The latter will be assumed to consist 
of symmetrical quadrupole triplets. Each triplet consists of a straight 
quadrupole section of length = followed by a drift space of lengthd, 
then a second straight quadrupole section of length 2 which is rotated 
90° with respect to the first quadrupole section. The second quadrupole 
section is followed by another drift space of length d, and a third 
straight quadrupole section of length 2/2 which is rotated a further 90° 
with respect to the second quadrupole section. The infinitely long 
classical channel constructed using such symmetrical triplets, will 
consist of a straight quadrupole section of length > followed by a 


succession of drift spaces of length, d, and quadrupole sections of length, 


2. The successive straight quadrupole sections will be rotated 90° with 
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respect to each other. 

Both the twisted and the classical structures to be compared will 
be assumed to have the same aperture radius, a, and the same electric 
field gradient. The classical channel has two symmetrical triplets in 
the same axial length as one complete twist of the twisted structure, as 


Shown schematically in figure (6.1). Thus, 


L = 4(2 + d) ore) 


where, d = length of the drift space between the quadrupole sections of 
the classical channel. 
2 = length of a complete quadrupole section of the classical 
channel. 
L = periodic length of the twisted structure. 


The fixed (X,Y,Z) coordinate system used in the previous chapters 
is used in the comparison. The classical and the twisted structures are 
oriented such that initially at Z = 0, the positive electrodes are aligned 
symmetrically about the X-axis. Thus the forces exerted on a positively 
charged particle at Z = 0, by either the classical or the twisted 
structures, are focusing forces in the X direction and defocusing forces 
in the Y direction. Therefore, the X-Z plane will be referred to as the 
initially focusing plane while the Y-Z plane will be referred to as the 


initially defocusing plane. 


thw noragnes se) aeizounee Ceotentt> a2 Sa 

,, iPatnae aake edd tine 6 zutber erudreqe ame) 
: ab Gietania: Teirssmays ows 20 Toners 

ga pbebbounte batetwt st: to tater atelgeay” 


e awit 4 (f 3) 


(7) 


: 
So anhlgsse Gloquibau od? nuanred soege JTW ant To eaone - Foals 
fannada favtvesto ats Ja 
featz2nfo ads Yo Wwortese si cquibsup eialguas 6 Yo nspaal= 2 _ 
onan = 
-yutouere bartind ons im naqnat athorveg = J _ 
arsiqens 2u0tve s ms i! breu mereye eran iinnaon {5,7.2) bextt et i 
eye eowiputse betziw? or ona jasleegts gut agg twqnes oy nt beau 2 


a 


- Ree | 
ait 4 = = de yh tort doue be2ns iio 


benpits ova eoborzasfe wy is rece 
; ee 
 Yisviatzoq.s san baivees payvot sili Sud? calmed ey ayode git : 
tadetat sit 0 livreesis ada vadtTo ila 16 liste é — 

sear? paleuadtsh brs nottoenth f oni ni tensa gibasod? ave ;: i 

: ; : , ie 

od ap ad eerste ot [TH smote + st erate aici my ane! 
sft 28 of Deven xi (ftw anely 5-° a7 a) (iw aaafg gaPawao? tina 

: - a | sinks pitauaateb | 4 #-3 hi 


| i 29 


—_— 


i. 7 
- a 
a > © 
- ‘aos 


. 7 
© 


ei) 


6.2 Acceptance Limits 


6.2.1 The Classical Structure 


If the axial momentum of a particle guided by the classical channel 
is assumed constant and equal to its value at injection, he the 
projections of the particle motion on the X-Z and the Y-Z planes are 
uncoupled. Therefore these projections can be independently analyzed. 
in, thes X%X-2 plane: the.analysis of the motion of a particle along the ith 
triplet gives the transverse coordinate x. and the normalized transverse 
momentum x of the particle at the exit plane of the triplet in terms of 
the transverse coordinate Keo and the normalized transverse momentum Xo 


of the particle at the entry plane of the triplet as 


Xi Xio Xio 
ee) aR tn (6.2) 
X Mal fe 
where, 
xr = = (atid ee 
Z dZ/dt mn Ps 


[T] = transformation matrix of the whole triplet in the X-Z plane, 


LT, J = transformation matrix of the first half of the triplet in 
the X-Z mane 
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[T.] = transformation matrix of the second half of the triplet in 
the X-Z plane(3”) 
ye de Sa 8 
fore SIO Weiieas ealal coshotes Lewsanh 
z y 2 
Bere ty yk ee ye 
y Sin 5 COs + 0 ] y sinh 5 cosh 5 
(6.3b) 


V 
r-ivVy? (6.4) 


where a, V and ve are the same as defined in Chapter 2. 
It is readily shown by performing the indicated matrix multiplication 


that the transformation matrix [T] can be put in the form, 


cos wu — sin u 
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“Yy sin u cos u 


where, 
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Az 
{sin y% cosh y2 - sinh y2} + yd{sin y2 sinh y2 - (1 - cos y2)cosh y2} 


a pee 
{sin y& cosh y2 + sinh y2} + yd{sin y% sinh ye + (1 + cos y%)cosh y2} 


Die 
- 1s (1 - cos y&)sinh yz Ve 
* al a 0 ena al een 
+ 1s (1 + cos y&2)sinh y2 (6.5) 


Then provided the stability condition 


=< COS 
or 
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is satisfied it is apparent that in the X-Z plane the triplet exerts an 
overall focusing force on the particle. The projection of the particle's 
motion on the Y-Z or initially defocusing plane can be described in an 
analogous manner. Indeed, provided equation (6.6) is satisfied an 
overall focusing action also will result in the Y-Z plane. 

A long chain of such triplets can be used to guide charged 
particles provided the components of the initial displacement and 
momentum along the X and Y axes are suitably restricted. In order that 
particles be confined within the structure it is necessary that the 
maximum displacements in the X and in the Y directions be less than, or 


at most equal to, the aperture radius of the structure. Thus, since the 
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<a and |Y| a: 


motion in the X-Z and Y-Z planes is uncoupled, |X| ae 


max 
Smith,L. et a, (38) have shown that these conditions will be 


fulfilled when the following conditions are satisfied: 


X xX! 
042 0 42 
eee (6.7) 
and 
042 Oo 42 
= ] 
gis sawas.* (6.8) 
where, 
Ky = initial X coordinate of the particle at injection, 
Py 
X! = = Py and Ps are the injection X-directed transverse 
Z O 0 
° momentum and the injection axial momentum respectively, 
ie = jnitial Y coordinate of the particle at injection, 
Py 
ue = oe Py is the injection Y directed transverse momentum, 
fa 6) 
y? = 22 
a 
11 


The parameters ayy and Aono are the two diagonal elements of the trans- 


formation matrix [T,] of the first half of the triplet in the X-Z plane, 
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Thus, 


9 COs t* cosh i+ sin <* sinh = + yd cos X= sinh ne 
y& yh yu 
2 Z 


cos 1 cosh 5 - sin 
(6.9) 


In the limiting cases equations (6.7) and (6.8) describe the elliptic 
boundaries of the acceptance areas of the classical channel in the X - X' 
and the Y-Y' phase planes. For a particle to be successfully guided by 
the classical channel its injection conditions in both the X-X' and the 
Y-Y' phase planes should be within these elliptical boundaries. A 
particle with injection conditions outside these elliptical boundaries 
will have a maximum transverse displacement exceeding the aperture 
available and the particle will be lost. The areas bounded by the 
acceptance ellipses in the X-X' and the Y-Y' phase planes are equal and 


given by 
A Zep, Ey ts (6.10) 
a ae ay ; 


For the special case of parallel injection along the classical 
structure, the maximum acceptable transverse injection conditions are 
independent of each other and their values are obtained from equations 


(6.7) and (6.8), by setting Je ey and Y' = 0, as 


|X | = 4a ’ he 2 


o'max Coa) 


Therefore, the acceptance area of the classical structures in the X-Y 
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transverse coordinate plane is a rectangle whose area is 
Ay _y 5 = » where » is given by equation (6.11). (6.12) 


Similarly, for the special case of point source injection at Z = 0, 
the maximum acceptable values of the normalized injection transverse 
momenta are independent of each other and are obtained from equations 


(6.7) and (6.8), by setting IX, | = 0 and Yo! = 0, as 
[XO BREW ye Ys LYS = wyya (6513) 


Therefore, the acceptance area of the classical structure in the X'-Y' 


normalized transverse momenta plane is a rectangle whose area is 


C Ys 


Rey eal gyi (6.14) 


6.2.2. The Twisted Structure 


For the twisted structure values of ae and Vase cannot be specified 
separately as was the case for the classical quadrupole channel. In order 
that the particles be confined it was seen in Chapter 3 that the radial 
displacement must be at most equal to the aperture radius, a. Thus 
Brae oe For particles confined within this circular aperture, equation 
(3.10) gives the acceptance limits of the twisted structure and shows how 


the maximum acceptable value of any of the injection conditions of a 


particle depends on the values of the other injection conditions of this 
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particle. It is difficult to compute the acceptance limits for particles 
with general injection conditions and the comparison between the acceptance 
limits will be restricted to the special cases of parallel injection, 
point source injection, injection in the X-Z plane with X directed 
transverse momentum and injection in the Y-Z plane with Y directed 
transverse momentum. For these special cases the acceptance limits of the 
classical and the twisted structures in the X-Y transverse coordinate plane, 
in the X'-Y' transverse normalized momenta plane, in the X-X' phase plane 
and in the Y-Y' phase plane will be compared. The twisted structure will 
be assumed to have (ga)? = 0.1 in the computations which will follow. 

The stability condition for the particle motion along the twisted 
structure, given by equation (3.3), can be written in terms of the 


normalized transverse momenta as, 


1 wv 
aheul le?) 
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ls| = 67a iy tt oe eee 

For the special cases of injection to be considered in the comparison the 
term(X oY 6 - Volos = 0. Also it was shown in Appendix B that for all 
particles successfully guided by the structure, B(X Yo - ey Ie 


Therefore, for all acceptable injection conditions the stability condition 


can approximately be written as 
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Is|_ = 2 [e-)° =] (6.15a) 


(6.15b) 


It is interesting to note that for particles guided by the classical 


channel equation (6.6) leads to the condition 
0< w eye (6.15c) 


when there is no driftspace. 
For particles injected in the X-Z plane with X directed transverse 
momentum, the acceptance limits are given by equation (3.14) which in terms 
of xy and XO becomes 
fo e LE lall a eed al ¢ ko 2 


ea fae) (6.16) 
(Soe) as. - | 


eae 


Equation (6.16) represents the intersections of two hyperbolae and gives the 
maximum allowable values of [Xo and [Xol. for particles injected with 
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| Ae | eae SYS |: 6 = 8 when eae 
Ws is Hl © shi oa 


The acceptance area bounded by the hyperbolae described by equation (6.16) 


I> 
I 


a 
T 
O 


: ay. X 1/2 
4ga [d5vs i) Penis = 1),/S 241 Tee 91 dX 
; E 


+ ivsVs ae 1) R= s+ St hnaé Gare 


where &£n denotes the natural logarithm. 
For the special case of injection in the Y-Z plane with Y directed 
transverse momentum, the acceptance limits are given by equation (3.15) which 


in terms of ies and ee becomes 


Acreear hi < | (6.18) 


representing the intersection of two parabolas. Equation (6.18) gives the 
maximum allowable values of a and Hele for a particle injected with 
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The acceptance area bounded by the parabolas described by equation (5.18) 
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For parallel injected particles the acceptable injection points 


(X03 i Should be within the parallelogram shaped area bounded by the 


straight lines, given by the limit of equation (3.12), 
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and the acceptance area of the twisted structure in the X-Y transverse 


coordinate plane is 
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Wiel aS (6.21) 


For point source injection, equation (3.13) can be rewritten in 
terms of Xo and ss to give the acceptance limits of the twisted structure 


as 


“|[s¢ x24 s(s = 1)yi7] | 
< |Ba| (6.22) 


and the acceptance area of the twisted structure in the X' - Y' phase plane 


is 
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calculated through numerical integration of equation (6.23). 


6.3 Acceptance Comparison 


The first objective of this section is to compare the acceptance 
limits of classical structures having f= 0. ,70.125-and 0.25 with those 
of the twisted structure. Computation of the acceptance areas in the 
X - X' and the Y - Y' phase planes is performed over the stable range of 
Vy and hence the stable range of L for each of these four guiding 
Structures. As mentioned before, the length of each of the symmetrical 
triplets constituting either one of the three classical channels is equal 


to half the periodic length of twist of the twisted structure. Thus for any 


value of “ y2 and Lt are related by 


= are lees S) (6.24) 
For the classical structure with no drift space, c= Oc Roa = 

and from equation (6.15c) 
qe aes (6.25) 


4 


For the classical structures with c= 0.125 and 0.25 the stable nange oF 
y& was computed using equation (6.6), then equation (6.24) was used to 
compute the corresponding range of Las 


0.0 < HE < 1.89 for © = 0.125 (6.26) 


s2nsagasze. oid sigamo> af 2f aortsaee oti 4: i" ido 92" + 3 ie - 
Sed AWAY bos OS, Oo . © patusn aavutbunte Yeotiont> teh 
aid ni 260% soneiqa2e eft Yo EL ail -smsoutte beretw2 oft 7 
qo spray sidase aft ovo harretieq 2f conn om ‘y = ¥ aA? bas 14 aie 
patbtug wo? sesds 30 dono or =} 6 apne’ stdete oft eonad br ra ' 
reotatemime ets 0 toes To flzpns! amy ,wr0ted banorinem 24. 2e5ud3 “— 
Kndpa 2t el onnatio Issiz26to seade arid To seg ‘youtia ont tutttznos sofa 
aypnawrs2 betefw aft to Jetwt 70 Alene! aibotysq ~ Isic : 


cd natetor one 2E bas ay ss 


(es .3) | . é + Fisy = 4 _ 


x = wy 4 te -2 , wOGGe FiIMb on arte wren? Tadrrzai> ony ot - 
(d81.0) norieups mon ba 


ye Wor <2uilT 


; a : 

(28:8) oa. oe 

: : 7 7 

to aot etdete sit a%.0 te OS1.iF = 2 frvtw eaiuitounde ets it 7107 
ot beew dew ABS.5) ootseupe ment a sj nite porn \ gmca 26W 


zn Af wit Spive = nes1102 ahs 


SSE 


asi, d= B yo 


Fee sO? 8'L 


eee 
eal) 


9'L 


“Saueld aSeud ,A-A pue ,X-X BY} UL SPduy adUeZdsd0y BU] 


A! él Ol 3" 


(w2 GL =0) 


9° y Z" 


PAA wg XK 


Daily adunjdad.y 


(2°9) aunbl4 


6) 


OL’ 


SL" 


Ot 


Soe 


Ot” 


Ce" 


i 
i 
i 


128 


and 
o<fF <i. for S= 0.25 (6.27) 
The acceptance areas in the X - X' and the Y - Y' phase planes, 
Ayes here er and AS ie were computed using equations (6.10), (6.17) 


and (6.19). The computations were performed for the twisted structure and 


the classical structures having cs O , 0.125 and 0.25, over the stable 
ranges of fa given by equations (6.15b), (6.25), (6.26) and (6.27) 


respectively. Figure (6.2) shows these acceptance areas drawn versus t. 

It is clear that the increase of ¢ from 0 to 0.25 decreases the acceptance 

areas Aon and af of the’ classical Structures im tne trange 0 =< = 

< 1.54 radians while it increases them over the range of t. higher than 

1.54 radians. It is also clear from figure (6.2) that for values of = 

up to 1.0 the acceptance areas ae and hee of the twisted structure are 

nearly equal to those of the classical structures. For the range 1.0 < 

= < 1.3, the acceptance areas of the classical structure having o = 0 

are larger than the X - X' phase plane acceptance area of the twisted 

structure by less than 6% while they are larger than the Y - Y' acceptance 

areas of the twisted structure by less than 12%. For the range ae les 

the acceptance areas of the classical structure are much larger than those 

of the twisted structure since the latter drop rapidly for B cae |e 
Equations! (6.5) .s(6.7)9_ 06.8) .) (6.9).29(6.16) and (6118) were used to 

compute the acceptance limits in the X - X' and the Y - Y' phase planes for 

the twisted structure and the classical structure having c= 0. These 


acceptance limits are shown in figures (6.3) and (6.4) for three different 


values of quadrupole voltage corresponding tos = 1.5, 2 and 3 = or 
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ie 1.283, 1.11 and 0.907. These figures show clearly that even though 
ae and eee are slightly larger than Mae and Se the acceptable 
maximum values of the normalized transverse momentum at injection, 


xs 


o'max and ne 


mae ore considerably larger for the twisted structure 


than for the classical structure. Note also that for s = 3, |Y for 


alee 


the twisted structure is slightly larger than |Y for the classical 


slifes 


channel. For s = 2 and s = 1.5 the values of |Y for the twisted 


elmax 


structure become smaller than the corresponding values of |Y for ste 


nie 
classical channel. 


Equations (6.7), (6.8) and (6.9) were then used to compute the 
X - X' and the Y - Y' phase plane acceptance limits for the classical 


Structures having c 0.125 and 0.25 for one value of tt - 1.11 which 


corresponds to s = 2. These acceptance limits, as well as those previously 


computed for the twisted structure and the classical structure having c = 0 


at si=%2 as tt = 1.11 radians, are shown in figures (6.5) and (6.6). It is 


clear from these figures that as < increases the acceptable normalized 


and |Y' 


Pere decrease for the classical structure. 


transverse momenta, IXo | 


max 
a Vb ge : ' 
For > = 0.25 and q~ = 1.11 radians, the value of Lee acceptable by the 


classical structure is 0.0499 while that acceptable by the twisted structure 


is 0.0599. Also the value of base 


having ¢ eae Sewer - = 1.11, is 0.0679 while that acceptable by the twisted 


structure 1s 0.079. 


acceptable by the classical structure 


For the special case of parallel injection, equations (6.9), (6.11) 
and (6.20) were used to compute the acceptance limits in the X - Y 
transverse coordinate plane. The computations were performed for the twisted 


structure and for the classical structures having a= Oe 0s l25 and UN2s 
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at s = 2 and the corresponding te = 1.11 radians. The computed acceptance 
limits are shown in figure (6.7) where the axes of the transverse coordinate 
plane are normalized with respect to the aperture radius. The areas in the 
X - Y plane available for parallel injection along the twisted structure and 
the classical structures considered were computed using equations (6.12) 


and (6.21) as 


eal ae, 

Avy = 2.8932 a” forf=0  , 
Ay y = 2.906476 a° for $= 0.125 , 
Ay_y = 2.93595 a for © = 0.25 


It is clear that the acceptance area in the X - Y transverse coordinates 
plane for the twisted structure is much less than the acceptance areas in the 
X - Y plane for classical structures with or without drift spaces. This 
situation is due to the coupling of the motion in the X - Z plane and the 

Y - Z plane of the twisted structure. Recall that in the case of the 


classical structure |X|_.. < a and LY | < a independently while for the 


twisted structure (xé a y*) yet a’. Of course the classical case is 


max 


somewhat idealized since it completely ignores all fringing effects. 
For the special case of point source injection at Z = 0, equations 
(6.9), (6.13) and (6.22) were used to compute the acceptance limits in the 


X' - Y' transverse normalized momentum plane. These acceptance limits are 
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Figure (6.8) Effect of Drift Space on the X'-Y' Acceptance Limits of 
Classical Structures. 
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shown in figure (6.8) for the same quadrupole potential, s = 2 , and the 
Same structures as considered in the case of parallel injection. The 
acceptance areas in the X' - Y' plane were computed using equation (6.14) 


and numerically integrating equation (6.23) using the approximate 


expression 

at j=0 100 O74] 0; 

| 6 lmax 
where ye is the value of uA evaluated from equation (6.22) at XS = 700 


The acceptance areas in the X' - Y' plane were found to be: 


Ayi ys = 0.0149 

Avy _y: = 0.0149 for $= 0 
Ayiys = 0.0145 for 2 = 0.125 
Aveuy: = 040136 for 2 = 0.25 


Based on the degree of accuracy of the numerical integration the 


areas inal and oe are essentially equal for f= 0 while for the 


practical values of ¢ 0.125 and 0.25 Aneve is larger than Avie: 


ol max auy 


Figure (6.8) also shows, as was pointed out earlier, that |X 


Ay are significantly larger for the twisted structure than for the 
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The above comparison shows that the X-Y acceptance of the 
twisted structure is smaller than that of the classical structure with 
or without drift spaces. However, the comparison shows that the twisted 
structure allows more divergence of the beam than that allowed by the 
classical structures, especially those with drift spaces. This larger 
acceptance area of the twisted structure in the X' - Y' phase plane 


makes it very suitable for guiding beams emerging from point sources. 


6.4 The Focusing Strength 


The focusing strengths of the twisted structure and the classical 
Structure having f= O can be compared by tracing particle trajectories 
along both structures. The projections of these trajectories on the 
X-Z and the Y-Z planes were computed and plotted for each particle traced. 
On the X-Z or the Y-Z planes, the projection of a particle trajectory 
along the structure which has a higher focusing strength in the 
respective plane should intersect the Z axis before the intersection of 
the projection of the trajectory of the same particle along the other 
structure. 

Equations (2.31), (2.35) and the transformation equations (2.7) 
were used to compute the X-Z and the Y-Z projections of the trajectories 
along the twisted structure. For the trajectories along the classical 
structure, the X-Z plane projections were computed using the matrix 
representation of equation (6.2) along each of the symmetrical triplets. 
The Y-Z projections along each of the triplets were computed using the 
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where [T,] and [T,] are as defined by equations (6.3a) and (6.3b). 

Particles with various injection conditions were considered in 
the comparison of the focusing strengths of the twisted and the classical 
structure. Some of these conditions were within the acceptance limits 
of one structure only. For all cases considered, it was found that the 
twisted structure has a higher focusing strength in either the X-Z plane 
or in the Y-Z planes than has the classical structure. 

Figure (6.9) shows the Y-Z projections of the trajectories of 
a particle injected at Z = 0 with Xo = ee = Xo = 0.0 and . = 0.07. 

These injection conditions are within the acceptance limits of both 
Structures. It is clear that the twisted structure has a higher 
focusing strength in the Y-Z plane than has: the classical structure. 

For a parallel beam injected at Xo = 5 and YA = 0, figure (6.10) 
shows the X-Z projections of the trajectories along the twisted and the 
classical structures. Note that the trajectory along the classical 
structure returns to X = 0 just after the return of the trajectory along 
the twisted structure to X = 0. Again, the twisted structure has a 
slightly higher focusing strength in the X-Z plane than has the classical 
suructune. 

It is clear from the above examples that the twisted structure 
is more strongly focusing than is the classical structure with no drift 
space. Moreover, because of the absence of drift spaces, the classical 


channel is itself stronger focusing than any practically realizable 
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CLASSIC CHANNEL 


Figure (6.9) The Y-Z Projections of the Trajectory of a Particle Along 
the BN Structure and Along the Classical Structure 
with oe 0, when (X, ay 0p tle = 0). 
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Figure (6.10) The X-Z Projections of the Trajectory of a Particle Along 
the Twisted Structure and Along the Classical Structure 
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classical channel. This stronger focusing action of the twisted channel 
indicates the possibility of shorter focusing structures than are 


realizable with a classical structure. 
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CHAPTER 7 


THE ELECTROSTATIC TWISTED QUADRUPOLE AS A NOVEL IDENTICAL IMAGING LENS 


7.1 Introduction 


The analytic solution for charged particle trajectories along the 
uniformly twisted electrostatic quadrupole was developed in Chapter 2 
using two approaches. These were the contact transformation approach and 
the linear algebra approach. The two different forms of the solution 
were then shown to be equivalent in Appendix C. The first solution, shown 
in equations (2.31) and (2.35) was successfully used for studying the 
guiding properties of a long uniformly twisted electrostatic quadrupole 
channel. The second solution, shown in equations (2.49) and (2.50) and 
describing the particle trajectories in a matrix form in terms of the 
injection conditions is in suitable form to show that a finite length of 
the twisted structure can serve as an imaging device. 

Examination of the matrix Q(t) of equation (2.50) shows that all 
elements are dependent on fit, fot as well as s, BP, and 2qMk. Through 
proper choice of the structure parameters the values of i and f, can be 
adjusted so that at a certain time ty all nondiagonal terms of the matrix 
will vanish while the four diagonal terms will all be + 1.0 or all be 
- 1.0. Hence,at this time t,,and at the corresponding length of the 


structure,the transformation matrix Q(t) will be either +I or -I, where I 
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is the 4x4 identity matrix. 

It will be shown in detail that the twisted lens structure with the 
‘properly chosen parameters will provide identical imaging of the injection 
conditions at Z = 0, the lens injection plane, into a target plane at 
Z = 2, the lens exit plane. Identical imaging means that the exit 
coordinates and momenta in the target plane are equal, with or without 
inversion, to the coordinates and momenta at injection. The target plane 
has its axes rotated an angle g2% with respect to the fixed (X,Y) axes that 
coincide with the electrodes axes at Z = 0 as defined in Chapter 2. 

The study of the imaging of a single particle by the twisted lens 
is followed by a discussion and a numerical study of the aberration that 
takes place when a beam of particles is considered. It is shown that the 
degree of aberration is small and that it can be further reduced by 
following the twisted lens structure by an equal but oppositely twisted 
lens structure. 

This chapter then presents a numerical study of the effect, on the 
imaging properties of the lens, of detuning the quadrupole voltage, about 
the exact imaging value. The study also discusses the characteristics of 
the exit beam cross-section for various values of drift space beyond the 
end of the lens. This last part of Chapter 7 forms the basis for the 


experimental investigations of the lens properties described in Chapter 8. 
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the twisted structure is given by equation (2.49) as, 


x) xe 
y(t) ‘p 
Et ya 
p, (t) , 
Xo 
pit) Py. aa 


where Mail ga Py and Py are the initial transverse displacements and 
a) fe) 


their initial conjugate momenta at t = 0, Z = 0 with respect to the 
rotating coordinate system (x,y,z). x(t), y(t), p(t) and py(t) are the 
transverse displacements and their conjugate momenta at time t, also with 
respect to the rotating coordinate system. For the convenience of the 
reader the matrix Q(t) is reproduced on the following page and renamed 
equation (7.2). 

It is now clear that a particle's coordinates and momenta will 


be imaged in the rotating coordinate system at some distance along the 


structure, if at some time, t = ty» 


Q(t, ) eo reas) 
where I is the 4x4 unit matrix. If Q(t, ) = - I the particle's coordinates 
and momenta will be imaged in an inverted form. In other words x(t) 5 oh 


y(t,) = "os p(t,) = Py and py (ty) = sae The conditions which must be 
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imposed on the elements of the matrix so that it becomes either +I or 


-I are: 
cos fit, = cos fot, le ae A 
sin fit, sin fot, = 0 
or 
fy ty a mr 9 
fot = nt (7.4) 


where m and n are integers, both odd or both even so that the term 
Bp 
Z < . : 
DqMk *©S fot, - cos ft, =/0 =. Substituting for fF, and fo in terms of 


the stability factor s, equations (7.4) become, 


\/ 8 Sick t= mar ‘ 
VS re ar oe 
M S er i haa nt 


It follows that for imaging to occur at t = ty» the stability factor s 


must have the value 
2 2 
oe (7.5) 


and m >on. Also, since the magnitude of s must be greater than unity, 


neither m nor n can be zero. 
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Thus, when the parameters of the twisted quadrupole structure are 
chosen so that s assumes one of the values indicated by equation (7.5), the 
particle will be identically imaged in the rotating coordinate system at 


a time 


1 
| Tact (7a) 


Appendix B shows that the axial velocity Z during the traversal 
of the twisted structure differs only slightly from the axial velocity 


at injection, U, » and that the approximation Z = U, is within the same 
fe) fe) 
degree of approximation used in deriving the analytic solution of the 


particle trajectories. In fact,during extensive computer calculation of 
particle trajectories over a wide range of injection conditions no case 


was found where Z differed from U; by more than 3%. Taking these results 
: 0 
into consideration, Z is assumed to be constant and equal to the axial 


velocity at injection, Us . Thus, the value of 2 may be approximated by 
; 6) 
replacing Z with U, in equation (7.7). Thus, 
fe) : 
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giving, p 


Substituting for s from equation (7.5) gives the length of the imaging 


lens as 
ea el a 0. (7.8) 


At this point it has been shown only that it is possible to image a 
particle in the rotating coordinate frame, where p, (t) and py(t) are not 
true mechanical momenta. Now consider a fixed coordinate system 
(X 


Y,, Z) whose axes X,, Y, are rotated an angle gZ(t) with respect to 


(i ee (Saami 
the fixed coordinate system (X,Y,Z), Z(t) being the distance along the 


structure given by 


t 
Z(t) “i Z(t) dt (7.9) 
O 


The axes Xe ve of such a coordinate system coincide with the rotating 
coordinate axes, x(t) and y(t), at the distance given by equation (7.9). 
The transformation equations relating the two fixed coordinate 


systems (X,, Ya. 2) and (X,¥,Z) ‘are 
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Kes Xphecos BZ(t) - Y, Sinmel (cae, 


ae X, sin BZ(t) + Y, CoS mead) memes 


ae) 
i] 


y Sarr BZ(t) - NG Sin @Z(t) 2 2, 


ao 
i} 


Py sin @Z(t) + P 
‘ "e 


COSsec( ty) os foal 6) 


where Py and Py are the transverse mechanical momenta conjugate to x, 
t t 
and Y,. Comparison of equations (7.10) with equations (2.7) and (2.13) 


SHOWS that x(t), y(t); p(t) and py(t) correspond identically to Xs Ves 


Py ; Py , the latter being coordinates and momenta in the fixed coordinate 
£ ct 
frame that coincides with the twisting axes at an axial distance given by 


equation (7.9). 
Thus from equations (2.7) and (2.13), it is apparent that at Z = 0 


ee Px and Py. are identical to Xo? Ves i and oe Similarly 
equations (2.7), (2.13) and (7.10) show that in a fixed coordinate frame 
(X75 Yr. 


the lens, x(t) = Xy, y(t, ) = Yrs p(t, ) = i and Py (ty) = ie Thus 5 
for a structure whose parameters are chosen properly to satisfy equation 


Z), which coincides with the twisting coordinates at the end of 


(7.5), it follows from equation (7.1) that, 
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15] 
It is evident from equations (2.32) and (2.33) that the axial momentum at 
Z = 2 1s then exactly equal to P5 » the axial momentum at injection. In 
other words the initial fae displacements and the initial momenta 
of the particle are reproduced in the fixed target coordinate frame 
(X 


Yrs Zeraltiesfixed Xz» Yz axes coincide with the x, y axes at Z = 2. 


T? 

Thus, the target reference frame axes are rotated an angle 8% with 

respect to the X, Y axes which coincide with the rotating axes at Z = 0. 
It is now clear that the twisted electrostatic quadrupole lens 

can image a single particle. The usefulness of any lens, however, lies 

in its ability to image a whole system of particles. In this respect the 


twisted quadrupole will generally suffer from aberrations. For instance, 


s is given by equation (4.]) as 
s = (ea)*|y-I[1 + (sa) 52 - S53] 112 
‘ a Ps a P ; 


When a system of particles enters the quadrupole structure as either a 


point source or a parallel beam, then 
KP = YP = 0 j Oyen an) 


and s has the same value for all particles. The value of s can then 

be adjusted, for example by varying the quadrupole voltage Vo? to satisfy 
equation (7.5) for these particles. However, when particles are injected 
under general conditions of transverse displacements and momenta,s may be 
different for the various particles in the beam. Condition (7.5) may not 
be satisfied for every particle. Then,the nondiagonal terms of Q(t) do not 


completely vanish for every particle and aberration will result. From 
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equation (7.12) it is evident that this aberration will be minimal for 
particles injected close to the axis with small transverse momenta. 

Even though a system of particles may have a single value of s 
and hence may not exhibit the aberration described above, a different 
type of aberration may still be present. The latter appears because each 
particle, according to its injection conditions, may have its own variation 
of Z about U, , as described by equation (2.34). Thus, the required length 
of the lens e given by equation (7.7) will differ slightly for each 
particle. This in turn means that for each particle a different error is 
incurred when computing the length of the lens from equation (7.8). As a 
consequence the degree by which the nondiagonal terms of the matrix Q(t, ) 
will differ from zero varies for each particle according to its injection 
conditions, thus giving rise to aberration. 

Computer evaluation of many trajectories has shown this latter type 
of aberration to be negligible when compared with the aberration resulting 
from changes in the stability factor. This phenomenon will be illustrated 


by examples of the trajectories shown in the following section. 
7.3 Results 


Using the analytic solution and the transformation equations (7.10), 
sets of particles with a variety of initial conditions were traced through 
the twisted quadrupole to graphically demonstrate the latter's imaging 
properties. In each case the numerical computations were performed with 
the aid of an IBM 360 computer. A listing as well as the discussion of 


the program used are given in Appendix E. The trajectories were plotted 
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directly by a Calcomp Plotter in the (X,Y,Z) fixed coordinate system. 
Although trajectories of charged particles through such an electrostatic 
Tens are independent of the charge to mass ratios of the particles, for 
the purpose of computation,particles with a charge to mass ratio of 30 
coulomb per kilogram were considered. 

The lens was chosen to have a = 1.5cm and (ga)é = 0.1. Since the 
parameters m and n have to be integers, both odd or both even and 
m> n> Q, the lens was chosen to have m = 3 and n= 1. These values of 
m and n correspond to the shortest possible length of the imaging lens. 
Thus at time th» 

cos fit, = cos fot, = -| : 

and 


ape sin fat, = 


nn 

= 

— 

—h 

Gi 
! 


and Q(t,) = -I for these particles which have 


The calculations were carried out assuming that every particle 
entering the lens structure had been accelerated through the same potential, 


V, and that the angle of divergence of the injected beam was very small, 
such that each particle had an axial injection momentum, Ps = (2qmv) '/2, 
fe) 


The quadrupole potential, V., was adjusted so that the stability factor was 


On 
1.25 for particles injected from a point source or a parallel beam where 


equation (7.13) was satisfied. In this case, from equation (7.12) 
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oa 
(aged 2\V 
s.=..(6a) = (ga) || 
2qMV up 
The accelerating voltage was assumed to be 30 KV which for (ga)? Sal 


required a quadrupole potential of 2.4 KV. 

The required length 2 of the lens was calculated from equation 
(7.8) to be approximately 33.3 cm, assuming either point source or 
parallel beam injection. Hence, imaging should occur in a target frame 
whose axes (X7, Yr) are rotated with respect to the (X,Y) axes by the 
angle B& = 403° = 360° + 43°. 

Figure (7.la) shows the trajectories of a set of particles which 
were injected from a point source at Z = 0.0. Each particle was injected 
with the same value of a but a different value of Cee It is evident 


from figure (7.la) that the point source is reproduced at Z = 2. Figure 


(7.1b) shows a prejection of the trajectories in momentum space onto 


Y 
the fixed Ge a) normalized transverse momentum plane. For sake of 
Ly ~) 
clarity those trajectories for which Py < 0 have been omitted. Notice 


fe) 
that the particles are imaged in inverted form in the target reference 


frame. 

Figure (7.2a) shows the trajectories of a beam which was parallel 
at injection. Each particle was injected with the same value of Xo but 
with different values of ve Notice that the emerging particles are again 
travelling parallel to the Z axis at Z = 2. Figure (7.2b) shows a 
projection of the trajectories in coordinate space onto the fixed (X,Y) 
transverse coordinate plane. Again it may be seen that the injection 


coordinates are imaged in inverted form in the target reference frame. 
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X and Y versus Z for Point Source Injection. 
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Figure (7.2a) 


X and Y Versus Z for Parallel Beam Injection. 


158 


"UOL{IaCUT Weag LAL L[e4eq UOJ SazZeULPY0O) 71X49 puke uoL.AaLU] 


(4z°Z) aunb4 


Mufisatnl mood (eflexs? wi 2otenibood +1xd bee nokfostal (aS.%) avugtt - 


= 


tog 


It should be noted that neither of the two foregoing examples 
indicates the presence of aberration. Since the two examples respectively 
illustrate point source and parallel beam injection, the only source of 
any possible undetected aberration would be a variation in Z which, by 
equation (7.7), would mean that the length of the lens is not correct for 

all particles. 

Figure (7.3a) represents the injection and exit coordinates of 28 
sets of particles. The injection coordinates lie along the periphery of 
the rectangle whose sides are parallel to the X,Y fixed axes. Each cross 
(+) represents five individual particles with the same value of Xo and ie 
but with different transverse injection momenta at Z = 0. These transverse 
momenta are shown in figure (7.3b). They are represented respectively by 
five (+) marks at the center and the four corners of the square whose sides 
are parallel to the pi and pl anes. Figure (7.3b) also shows the exit 
transverse momenta. In general these particles do not have s = 1.25 and 
therefore cannot be imaged without aberration. Figures (7.3a) and (7.3b) 
clearly show the presence of this aberration which is due to the variations 
in both s and Z. 

Figures (7.4a) and (7.4b), which feature the same injection conditions 
as figures (7.3a) and (7.3b), show the effect of joining to the first section 
of twisted quadrupole a second equal but oppositely twisted section. For 
this composite structure the transformation matrix Q(t) in equation (7.1) 

TS (-1)° or +l. Therefore, the injection conditions are now reproduced in 
a target plane without inversion or rotation. 


The computation of the trajectories along the composite structure 


was started by computing the trajectories along the first section, using 
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equations (7.12) then equations (2.31) and (2.35). Then,for each 
individual particle the axial momentum, P_ = MZ at Z = 2.was calculated 
using equation (2.34). The fixed target plane of the first section co- 
incides with the fixed injection plane of the second section. Therefore, 
the coordinates and momenta at the exit of the first section, in the fixed 
target plane which is rotated 403° with respect to the (X,Y) axes, are 

the injection coordinates and momenta to the second twisted section. For 
each individual particle, its input transverse coordinates and momenta as 
well as its input axial momentum to the second section were used to 
calculate the value for the motion of this particle through the second 
section, using equation (7.2). Equations (2.31) and (2.35)were then used to 
calculate the trajectory of the particle through the second section. 

Figures (7.4a) and (7.4b) show that the injection coordinates and 
momenta are now reproduced without inversion in a target plane that co- 
incides with the fixed (X, Y) plane. It is evident that some aberration 
is still present. However, this aberration is greatly reduced when compared 
to the results obtained with the single twisted quadrupole and shown in 


figures (7.3a) and (7.3b). 


7.4 Effect of Detuning s About Its Imaging Value 
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The twisted structure acts as an imaging lens if the structure parameters 


and the particle's injection conditions are chosen to give a value of s 


equal to 
mé + n° 
an ai emma 3 Where m and n are as defined by equation (7.5) 
m -n 
—- owiep= teed 
Any deviation of V, V., (ga) or {—— —— - ~9 9) from the proper values 
0 a PS a PS 


corresponding to s = Soe will lead to Sberration. In practice, an error 

in the length of the imaging lens, an error in the value of the aperture 
radius or any deviation of the fields from those derived from the focusing 
potential of equations (2.1la) and (2.1b) may require adjusting the structure 
parameters to values other than those corresponding to s = S, $0 that imaging 
will occur. Any deviation from this new setting of the parameters again 

will lead to aberration. 

In the present section the length, aperture radius and fields are 
assumed to be exact and the effects of the variation of s, about Sy» on the 
imaging properties of the lens are investigated. A study is also made of 
how the emerging beam parameters are affected, for values of s equal to 
and not equal to So? by introducing a drift space beyond the exit plane of 
the lens. 

A sketch of the arrangement for performing the foregoing studies 
is shown in figure (7.5a). A beam of particles is assumed to originate 
from a point source at a distance Le from the injection plane of the lens 
used in the numerical investigations of section (7.3). Before injection 
the beam cross-section is shaped by using a rectangular aperture E 


centered at Xo = te = Q and having its sides parallel to the fixed X, Y 
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eee eee ! Teel. Freee 
Representation is in the —#— Representation is in the 


(X,Y,Z) reference frame (Xp ¥722Z) reference frame 
1 


Lens 4 particle 
injection exit detector 


plane plane 
f— 6, +} Twisted lens ———++— 8, 


Figure (7.5a) Schematic of the Arrangement Used for Studying the Effect 


of Detuning s. 


Lens 
[Inversion + 
Rotation 43° | 


0 
‘in (X,Y,Z) frame ; in (X,Y75Z) frame 
| 
Figure (7.5b) Lens Equivalent at Imaging Values of s. 


At s = 1.25 (V. = 1200 volts) the effect of the lens is 
inversion and rotation of the coordinates and momenta. 
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axes, with the longer side parallel to the X-axis. The dimensions of the 
aperture E are chosen so that the dimensions of the injected beam 
rectangular cross-section are 1 mm x 2 mm With this aperture the angle 

of divergence of the injected beam will be less than 0.038 radians for 

2. = 6 cm, and the injection axial momentum of all particles will be 
approximately given by a ~ J2qMv. Also, with this small aperture, the 
injection conditions of the particles will be within the limits set by the 
confinement criterion of equation (3.10). 


Figure (7.5a) shows that the coordinates and momenta of the 


particles at the injection plane are related by 


Py 
ee 2 
O 
Py 
Beles Pe (7.14) 
0 


Thus,for all injected particles 


Foe Me -Y P, ]=0 a 
fe) ue 0 Xo 


which together with the approximation Po =\ /2qMV will give s as 
fo) 


|s| = (Ba) 


Thus, for the arrangement of figure (7.5a) and for the lens structure 


considered in Section 7.3 for which Pee 1.25, the exact imaging quadrupole 
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voltage is ee = 1200 volts when the accelerating potential is 15 KV. 

In the numerical investigation, presented at the end of this 
section only particles on the circumference of the rectangular injected 
beam are considered. The coordinates and momenta of these particles at 
the exit plane of the lens are computed using the transformation matrix 
Q(t) given by equation (7.2). The time t in which the particle 


traverses the lens structure is approximated by 


HI 


Q 
ie gers 
With this approximation f,t, = NSE ANE i ee ube) 

V 
g = 33.3 cm as computed in Section 7.3. 


At the exact imaging voltage V5 = 1200 volts, corresponding to 


= 1.25, tne exit coordinates of the particles at the exit plane of 


w 
! 


0 
the lens will be related to the coordinates at the injection plane by 


ee a and Yy = - Y (7.15) 
exe exit 
Therefore, with no drift space beyond the lens exit plane, the exit 
coordinates of the particles will lie on the circumference of a rectangle 
having the same dimensions as the injection rectangle and having its sides 
rotated 43° with respect to tne X and Y axes, that is parallel to the Xr 


and Y, axes of figure (7.2b). Moreover, tne exit momenta will be 
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Z Z | 
P = - Pp = ox Vf SON Y satis) (716) 
ee US OF exit *i 
and thus [X> Py - Y Py ] = 0. Equation (2.32) gives 
exit vers exit ears 


the axial exit momentum for such particles as 


P Pee. (Fe 17) 
Loxit Lo 
Equations (7.16) and (7.17) give the relations between the exit coordinates 
and the exit momenta of a particle, when there is no drift space beyond the 


lens exit plane, in terms of the (X75 Yrs Z) coordinate system as 


Pp 
p ay lexit 
Xt Toxit 4 
exit 
P 
feng 
ih ee wale : (7.18) 
ane exit 7] 


Equations (7.15) and (7.18) show that the coordinates of a particle 


at a distance &, beyond the exit of the lens, as shown in figure (7.5a), 
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i 
ue 
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Thus, at s = is 1.25, thesintroduction of adriet space he 
beyond the exit of the lens will not change the geometrical shape of the 
beam from that at the exit. The length of the sides, however, is 


Q 
increased by the ratio (1 + a. Thus at s = si 1.25, the action of the 
1 


lens itself can be represented by the introduction of a rotation through 
2n + 43° and inversion of the coordinates and transverse momenta as 
presented schematically in figure (7.5b). 

Having considered the properties of the image produced by the 
experimental arrangement of figure (7.5a) at s = 3i 1.25 without and 
with a drift space de beyond the lens exit, a study was made to determine 
the effect of varying s by considering values of Ue other than 1200 volts. 
It should be noted that for particles whose coordinates and momenta are 
related by equations (7.14), the relations between these coordinates and 
momenta and those at the exit of the lens are linear for all values of V,. 
For values of Vo different from 1200 volts the matrix Q(t) at Z = 2 will 
be different from -I but the transformation described by equations (2.49) 
and (2.50) will still be linear and can be written, using equations (2.49) 
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Substituting for the elements of [Q(2)] as O53 i= lG2,244, 3 = l2soe4s 


the above matrix relation can be put in the form 


where the elements ees 


considered. 


combinations of Xo and Ve The coordinates of such particles, whose 
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Equation (7.20) shows X7, Yr, Py and Py to be linear 
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are evaluated from (2.50) at the value of Pe 


coordinates and momenta at the lens exit are given by equation (7.20), at 
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distance de from the exit plane of the lens are, 
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Substituting equations (7.20) into equations (7.21) gives 
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where Al» Ay» A. and Ay are functions of Ley &, and the structure parameters 


at the value of Fe considered. 
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a straight line such that Xo ve are related by 


in = b, Xo 7 bos where bj and b. are constants. (7223) 


Equations (7.22) and (7.23) give the relation between X; and Y, , at 
Q Q 
e e 


a distance de from the exit plane, for a particle whose injection co- 


ordinates are related by equation (7.23) as, 


“Pe MM Pg Pate ms Si 
Ty A, + bj A, yy 2 A, a By A, 


Y (7.24) 

Equations (7.23) and (7.24) show that particles injected with 
their coordinates lying on a straight line will have their coordinates at 
any value of drift space from the lens exit plane lying on a straight line 
defined by equation (7.24). The slope and intercept of the line defined 
by equation (7.24) is clearly dependent on the original orientation of the 
line defined by equation (7.23), as well as on es Le and the structure 
parameters. 

Exit coordinates and momenta of the particles on the circumference 
of the injection rectangle were computed for different values of s by 
considering values of Ve equals tom150, 1180331200, 121d, 1220; 12505 and 
1300 volts. The computations were repeated for values of de equal to 0.0, 
fede) 2.0 and 3.0.cmewhile &. was equal to 6 cm. The computed exit 
coordinates in the target reference frame (Xr Yq Z) are shown in figures 
WieGi)e (7e7a (7.8) and (7.9) for ies 0.0, 1.0, 2.0 and 3.0 cm respectively. 
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VG Se icodevOl te athat 1 Ses a= = 1.25, to be exactly on the circumference 


ao 
of a rectangle whose sides are equal to those of the rectangular cross- 
“section of the beam at the input plane and are exactly parallel to the Xy 
and Yr axes. As i is decreased below 1200 volts, corresponding to s 
increasing to values larger than 1.25, the dimensions of the sides that 
Should be parallel to the Xo axes decrease at a rate much higher than the 
rate of decrease of the other sides that should be parallel to the Yr axis. 
At the same time the exit rectangle is deformed into a parallelogram that 

is rotated in a direction opposite to that of twist of the quadrupole 
electrodes. The angle of rotation of the parallelogram increases and its 
dimensions decrease as V5 decreases. 

On the other hand, as Le is increased above 1200 volts the exit 
rectangle deforms into a parallelogram that is rotated in the same direction 
as the twist of the lens electrodes. At the same time the lengths of the 
Sides increase and the sides that should be parallel to the Xr axis increase 
in length at a rate higher than the rate of increase of the other two sides. 

Figures (7.7) to (7.9) show the exit coordinates of the particles 
injected at Z = 0 on the circumference of the injection rectangle, at values 
of drift spaces de beyond the lens exit plane, of 1 , 2 and3 _ cm.The 
value of hes shown in figure (7.5a), is 6 cm as it was for the case de = 0 
of figure (7.6). Figures (7.7) to (7.9) show that the sizes of the exit 
parallelograms, at all values of We increase as de is increased. The exit 
coordinates at ue = 1200 form a rectangle with its sides parallel to the 
eo edhiany 


iT T 
(7.18). The exit coordinates for values of va less than 1200 volts are 


axes for all values of de considered, as expected from equations 


parallelograms, as for Lo 0 , rotated in a direction opposite to the 
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direction of twist of the lens electrodes. It is observed that the rotation 
of these parallelograms increases as We is decreased and/or the value of 

he is increased. For values of ve larger than 1200 volts the exit 
parallelograms are still rotated in the same direction of twist as that 

of the electrodes for all values of the drift space de beyond the lens exit. 
It is also found that the degree of rotation increases as the length of 

the drift space de is increased and/or the value of Ls is increased. 

The investigations of this section are of great importance since 
they are the basis for some of the experiments of Chapter 8. The 
experimental arrangement used for investigating the imaging properties of 
the twisted lens is exactly that of figure (7.6a), and the results of the 


present section are used to check the experimental results of section (8.4). 
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CHAPTER 8 
EXPERIMENTAL INVESTIGATIONS OF THE IMAGING PROPERTIES OF THE TWISTED 
ELECTROSTATIC LENS 


The theoretical analysis for the motion of charged particles along 
the uniformly twisted electrostatic quadrupole structure shows that the 
trajectories and hence the guiding, focusing and imaging properties of 
the structure are independent of the charge-to-mass ratios of the particles. 
Therefore, either electron beams, ion beams or any other charged particle 
beams can be used in the experimental investigations of the structure 
properties. However, since the present work was aimed at providing a 
successful guiding structure for microparticle beams, a charged micro- 
particle beam was used in the experimental investigations. The microparticle 
beam was generated through contact charging in a microparticle charger 


specially designed to meet the requirements of the experiments. 


8.1 High Vacuum System 


The vacuum system that was used in all experiments is a CHA High 
Vacuum Pumping Station type SS-600-RAM. It consists of a 7 3/4 inch, 
4-stage water cooled oil diffusion pump with liquid nitrogen trap. 
Pressure measurements were made by means of two DV-6M thermocouple gauge 
tubes and one Bayard-Albert hot wire IG-100-N ionization gauge tube. The 


experimental chamber consisted of two pyrex glass crosses, one having 
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symmetrical arms both with inside diameters of six inches, the other 

having arms with inside diameters of four inches and six inches. A 

combination of aluminum and perspex end plates together with rubber 

O-rings was used for sealing the system. The usual operating pressure 
7 


ranged from 1 x 10 torr to 5 x 10° torr. A schematic diagram of the 


high vacuum system is shown in figure 8.1. 


8.2 The Microparticle Charger 


The charger was based on the technique of contact charging of the 


40). 


micron-sized particles: Electric fields are then used to accelerate 


these heavy charged particles. Several microparticle contact chargers 


40 Saat) to produce single particles 


have been designed by other workers | 
with the highest possible charge-to-mass ratios. The charging electrode 

in each of these chargers is usually a very small sphere, of radius less 
than 100 microns, supported by a thin tapered needle. The small size 

of the sphere, for a given sphere potential, assures a very high surface 
electric field and consequently results in high charge-to-mass ratios 

for the particles. On the other hand, the small size of the sphere also 
means that charged particles are emitted at a very low rate. Such chargers 
are best suited for providing single highly charged microparticles for 


(41) 


Studying hypervelocity impacts or for evaluating microparticle 


detectors 42), 

In the present experimental investigations high charge-to-mass ratios 
are not the main objective. As long as gravity effects are negligible and 
provided all particles have been preaccelerated through the same potential, 


particle motion in the structure is independent of charge-to-mass ratio. 


However, since the microparticle beam cross-section is detected visually 
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by intercepting the particles on a greased glass plate, the charger should 
provide a large flux of microparticles for long periods of time, irrespective 
‘of charge-to-mass ratio. 

The present microparticle charger has a conical charging electrode 
with a rounded tip. The large radius of curvature of the tip increases 
the probability of emission of charged particles. Thus, a high density 
microparticle beam is generated although the particles have charge-to-mass 
ratios much lower than in the case of small spherical charging electrodes. 
The charger has a microparticle reservoir from which particles are fed 
mechanically into the charging chamber. The reservoir system permits 
continuous operation of the charger for more than 12 hours. 

The charger is shown schematically in figure 8.2. The charging 
chamber consists of two brass circular walls, A and B, separated by the 
insulating perspex wall, C. The stainless steel cylindrical charging 
electrode with its conical charging tip, having a head angle of 60°, is 
screwed into the central threaded hole of wall A. Two nuts are used to 
set the distance that the charging electrode protrudes into the charging 
chamber. Vertically above the charging electrode is a hole R in wall A. 

P is the microparticle powder reservoir. It consists of a brass body to 
which a vertical glass tube is fixed. P is vertically mounted on the 
brass cylinder T and a vertical hole in T connects the interior of P to 
the hollow interior of T. A threaded auger O passes through the 
horizontal cylindrical hole inside T and protrudes into hole R. The 
auger is mechanically coupled to a teflon insulating shaft which in turn 
is fixed to a thin stainless steel shaft. The latter passes through a 


rotary vacuum feed-through on the left hand side aluminum end plate of 
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Figure (8.2) The Microparticle Charger. 
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the 6-inch pyrex glass cross. The shaft is rotated by a synchronous 
motor at a speed of 10 R.P.M. The charging electrode, wall A, and 

the conducting parts of the feeding mechanism are at a potential of 15 KV 
with respect to a common ground point to which wall B is connected. 

The microparticles in P drop onto the auger 0 through the vertical 
hole in T. As the auger is rotated it transports the particles toward 
the charging chamber between walls A and B. The electric fields that 
exist between the grounded brass peg G, which acts as a field intensifier, 
and wall A induce charges on the particles that have fallen into R. The 
charged particles in R are accelerated towards wall B. Striking either 
the peg G or wall B, these particles exchange charge and are repelled 
towards wall A. This process repeats itself many times and aids in the 
deagglomeration of the microparticles entering the charger. The tapering 
of the interior side of wall A causes the electric fields between A and B 
to have components that force the charged particles, undergoing alternate 
collisions between A and B, towards the center of the charging chamber. 
This action enhances the flux of particles in the vicinity of the charging 
electrode and, as described below, helps to increase the rate of emission 
of charged particles. 

Some of the particles at the center of the charging chamber will 
enter the region of high electric field between the conical tip of the 
charging electrode and the hemispherical depression in wall B. As the 
particles contact the tip of the charging electrode, each particle gains 
a high positive charge and is accelerated toward wall B. Some of these 
particles pass through the aperture while others continue to alternately 


collide with wall B and the charging electrode. The latter particles 
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Figure (8.3) Exploded Front View of the Charging Mechanism. 


Figure (8.4) Exploded Rear View of the Charging Mechanism. 
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may pass through the aperture subsequently or they leave the region of 
high electric field and settle on the interior of the insulating perspex 
icy linden,. GC. 

Since the emitted charged particles leave the tip of the charging 
electrode at various angles, it was found necessary to place a second 
aperture in the path of the particles. This second aperture, D, limits 
the divergence of the emitted microparticle beam and helps to avoid 
contamination of the vacuum system. Exploded views showing the components 


of the microparticle charger are shown in figures (8.3) and (8.4). 


8.2.1 Mlesting ofthe thangen 


The charger was tested to determine the breakdown voltage between 
the charging electrode and wall B. With no particles inside the charging 
chamber, a potential of 30 KV was reached without breakdown. During the 
experimental investigations the microparticle reservoir was filled with 
3 micron carbonyl iron microparticles. When these particles were fed 
into the chamber, breakdown occured at lower voltages. During these 
experiments the value of the breakdown voltage was dependent on the distance 
between the charging electrode tip and the hemispherical depression on 
wall B and also on the rate of feeding of the particles. Through trial and 
error it was decided to use an operating voltage equal to 15 KV in all 
experiments. 

It was then observed that after a period of operation, the auger 
ceased to feed the particles into the charging chamber although the 
particle reservoir P was still full. This failure was attributed to the 


small sizes of the particles that cause them to stick together and to the 
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Figure (8.5) Microparticle Charger and Tapping Mechanism. 


Charge Sensitive Amplifier 


Vout [To C.R.O] 


Faraday Cup 


Figure (8.6) Sketch of the Charge Detection Arrangement. 
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walls of the reservoir. By tapping on the auger cylinder, T, it was found 
that the feeding action could be restarted. During the experiments the 
‘tapping was provided by a teflon-iron plunger actuated by an electrical 
solenoid. Figure (8.5) shows the assembled microparticle charger with the 
tapping mechanism mounted below the auger cylinder, T. The circuit used 
to actuate the solenoid consists of a silicon controlled rectifier which 
is triggered by pulses from a pulse generator. The repetition frequency 
of the pulses was 10 cycles/sec and the width of the pulse was 8.33 msec. 
When triggered, the SCR connects the terminals of the solenoid to the 
secondary of a step-down transformer, whose primary is connected to the 
mains via a manual switch. The latter also turns on the motor that drives 


the auger mechanism. 


8.2.2 Charge Detection 


The objective of this experiment was to ensure that charged particles 
were emitted from the charger. The position of the charging electrode tip 
with respect to wall B, the rate of particle feeding, the rate of actuating 
the tapping plunger, and pitch of the auger thread, were all adjusted 
through trial and error until a reasonable rate of emission was obtained 
without breakdown. 

Since particles are not emitted from the charger as single particles 
but rather as a continuous stream it was not readily possible to carry out 
conventional time of flight avderinents manta determine the charge-to-mass 
ratios of the particles. However, it was possible to detect the charges 
of the particles by using a Faraday cup that intercepted the particles 


emitted from the aperture D. The Faraday cup was followed by a charge 
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sensitive amplifier whose output was fed to a Cathode Ray Oscillograph as 
shown schematically in figure (8.6). 

Figure (8.7) depicts the typical output of the charge sensitive 
amplifier and shows the wide range of the charges on the microparticles. 
This wide variation is due in part to the size distribution of the 
particles. Also, since the charge induced on each particle depends on the 
particle shape, any irregularly shaped particle will accumulate a charge 
different from that on a spherical particle of similar size. Even if all 
particles have the same size and the same spherical shape, two or more of 
these particles may agglomerate to form a larger irregularly shaped 
particle. The latter, when charged, will have a charge on it different 
from those on the single spherical particles. Another factor that may 
cause similar particles to have different charges is the dynamic nature of 
the contact charging process. Little is known about the latter, but it is 
possible that the time constant of this process is much larger than the 
time of contact between the microparticle and the charging electrode. 
Therefore, the microparticle could leave the charging electrode before 
the charging is complete. 

A typical value of the charge on the carbonyl iron particle, as 
calculated from the detected output voltage of the charge sensitive 
amplifier is 


14 


q =.56x10 ~ coulomb 


Assuming spherical particles of carbonyl iron with a diameter of 3 microns 


which equals the average size specified by the manufacturer, the mass of 


7 > | Oo 3 . 


gfe Yo nattedtiseth oxte ond ore 
ont Wo ebnaqel sfottysg Maes ao bsoubat 
opine 4 a¥eTumoae Miw efahiney bagete: 
fe t¥ weed) vowta veltmt2 to efsisieq 1 
io s70n-40 owt weusie Toatvedee omnz oils tne owe . ; 
Haqede ylaslupemt opal 6 mvt oF ‘sdllcheeiba wikeniepeti 
tnavels YF ee epvade ie est TT tw , berets nat settel edt valotaeg 
“sil Jeet toraey yoddanA .eataterag (a2 tasdqe stente atts no" s20Hs mov? 
Yo esuten otmdinh of 2f esprens Inovst th svat oF zafotivag vefimt2 s2u69 
it fed ovettel af% duode wont 2l-of3ts s2esooNgyphtp Bits S9e9n03" ord 
st? oetiy wept Moun zt 2zsoorg zt To tnedenoy eats, ait tend efdteeoq 
esbowsosle enfyteds sit Gna oi stiteqotstm ett noewted tosinoo to smbs 
sheted wbotizsfe ontoinn: siz eves! Divas slottreqorotm att ,erotetadT 
25 ,piotitsg novi fynodtss ot) no spedo xt to sbi sy Testqyt A 
svittense spiets silt to sgetfov tuatuo bafosieb ott mort bets fuates 
8) eh ve ptt ame: 
—— a » ets 
: dmoi(uoo *orxaaeege ,: ” peu) 
Gi ® a “ny 
enortotm € Yor retean?b 6 dttw nowt Iynodtss Yo estattieq: Teokverige potmeek 
Yo\ses afd terusostuvem ony ys ber tioege sxtaysgereve silt spe 


—y 


é 


T91 


Figure (8.7) Photographs of Oscilloscope Trace of the 
Typical Output from the Charge-Sensitive 


Amplifier. 

(a) Output recorded during a single sweep on 
the C.R.0. 

(b) Output stored during several sweeps on the 
Gak..Ok 


Sensitivity: 5 mV/cm, 1 ms/cm. 


Typically, q = 0.56.10°'* coulomb, M = 1.13.10" 
kitogram and ae 0.05 coulomb/kilogram. 
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a single particle is 


M= + r3 
Bai -6,3 
= = x 8000 x (1.5 x 1079) 
= 1.13 x 107!? Kilogram 


Thus the typical value for the charge-to-mass ratio of the particle 


produced by the charger is 


Te 0.05 coulomb/kilogram 
It is evident that for ve = 1200 volts the effect of gravity on these 


particles will be negligible. 


8.3 The Uniformly Twisted Quadrupole Structure 


8.3.1 The Electrode Shapes 


To provide the desired quadrupole field as given by equations (2.1a) 
and (2.1b), the electrodes must be shaped to coincide with a set of 
hyperbolic equipotential surfaces extending to infinity. If electrodes 
of finite size are used the field will be distorted. However, in designing 
the poles or electrodes for straight classical quadrupole sections it has 


(44) 


been found that even ey a) or semi-circular electrodes can give 


acceptable quadrupole fields provided their dimensions are suitably 
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The electrodes for the uniformly twisted electrostatic structure 
are most easily manufactured by producing a cylindrical mould with the 
required electrode shape. As it is no more difficult to machine a mould 
of hyperbolic cross-section than to machine a mould of circular cross- 
section, it was decided to use moulds of hyperbolic cross-section. In 
order to reduce the amount of truncation distortion of the fields within 
the quadrupole structure aperture, the outer surface diameters of the 


electrodes were made greater than 3 times the aperture diameter. 


8.3.2 Structure Parameters 


The aperture radius and the periodic length of the twisted structure 
must satisfy the condition (ga)? << 1. Choosing (ga)? = 0.1 and the 
aperture radius a = 1.5 cm, the periodic length of twist is L = 29.68 cm. 
To test the imaging properties of the twisted electrostatic lens, the 
parameters of the structure were chosen to correspond to the minimum 
possible length of the lens. It was shown in Chapter 7 that this minimum 
length corresponds to m = 3 and n= 1. These values in turn correspond to 
s = 1.25 and give the length of the lens as 2 =) = 33.3 cm. Equation 
(7.12) shows that for these values of s and (ga)¢ and an accelerating 


potential of 15 KV, which is provided by the microparticle charger, the 


magnitude of the quadrupole focusing voltage must be Us = 1.2 KV. 


8.3.3 Fabrication of the Twisted Structure 


A mould is used to fabricate each of the twisted electrodes of the 
structure. The mould, shown in figure (8.8), consists of a cylindrical 


aluminum split sleeve which encases a machined aluminum cylinder. The 
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Figure (8.8) The Mould that was Used to Fabricate the Twisted 
Electrodes. 


Photographs (a), (b), (c) show the mould for the 
clockwise twisted lens at different stages of 
assembly. The electrode is still inside the mould. 
Photograph (d) shows the mould for the counter- 
clockwise twisted electrodes. 
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latter is manufactured by first mounting the aluminum cylinder on a 
milling machine with its axis parallel to the horizontal bed of the 
machine. Then, while the cylinder, through an appropriate gearing 
mechanism, is slowly turning about its axis and while the milling machine 
bed is moving horizontally,a vertical cutting tool mills an initial small 
groove. The vertical position of the cutting tool and the initial 
angular position of the cylinder are then precisely adjusted for the milling 
of a second groove. The process is repeated many times until a mould with 
the desired hyperbolic cross-section is completed. The direction of twist 
of the final electrodes is determined by the direction in which the aluminum 
cylinder is turning during the milling process. The outer surface diameter 
of the machined aluminum cylinder is 9.3 cm, greater than 3 times the 
aperture diameter. The length of the cylinder is approximately 2 cm longer 
than the required length of the imaging lens in order to allow for the 
Shrinkage of the moulding resin as it solidifies and in order that the ends 
of the electrode can be finally machined to the exact required length 
with the end faces exactly perpendicular to the axis of the cylinder. 

Once the mould is machined, the fabrication of the electrodes is 
started by spraying the mould with a standard release agent. The sleeve 
is then placed around the mould. The two liquid components of a general 
purpose rigid low temperature potting resin are then mixed together. The 
mould is filled with the resin and it is left to cure under vacuum. Once 
cured the solid electrode is withdrawn from the mould, the release agent 
is cleaned off and the electrode's ends are cut to the prescribed length. 
The surface of the insulating electrode, shown in figure (8.9a), is then 


sprayed with a conductive silver paint. Four silver painted electrodes 
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Figure (8.9a) A Clockwise Twisted Electrode Before Painting. 


Figure (8.9b) Front and Side Views of the Clockwise, 8 - ve, 
Twisted Lens Before Painting. 
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Figure (8.9c) The Uniformly, Clockwise, Twisted Lens After Silver 
Painting. 
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are then assembled to form the twisted lens. The electrodes are held in 
place by using perspex supporting sleeves. Figures (8.9b) and (8.9c) show 
the assembled lens before and after silver painting. 

Two moulds were machined, one for the counter-clockwise twisted 
electrodes, 8 positive, and the other for the clockwise twisted electrodes, 
B negative. For each of the twisted lenses opposite electrodes are 
connected by a high voltage cable. Two other high voltage cables supply 
the focusing voltages + Vs and - ue to two adjacent electrodes. The 
focusing voltages + VA and - Ve are supplied by a high voltage supply 
built for this purpose. A circuit diagram of this supply is shown in 
Figure (8.10). 

A breakdown test of the twisted lens was performed with one opposite 
pair of electrodes grounded and the other pair at high potential. Breakdown 
between adjacent electrodes occured at a potential of 3.5 KV. There was 


no breakdown at the required operating potential difference of 2V = 2.4 KV. 


8.4 Experiments 


The schematic diagram of figure (8.11) shows the arrangement used 
to investigate the imaging properties of the twisted lens. The micro- 
particle charger and the twisted lens are mounted on a horizontal 
insulating slab inside the vacuum system. The horizontal axes of the 
charger and the lens are carefully aligned to avoid aberration caused by 
misalignment. The twisted structure is oriented such that the electrode 
axes at the injection plane, Z = 0, that is the fixed X and Y axes of the 


theoretical analysis, are in the horizontal and vertical directions. 
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Representation is in the Ee Representation is in the 
(X,Y,Z) reference frame (X75¥ qT?) reference frame 
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Figure (8.11) Sketch and Photograph of the Experimental 
Arrangement Used to Test the Imaging 
Properties of the Twisted Lens Structures. 
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The imaging properties of the twisted lens are examined by shaping 
the cross-section of the injected microparticle beam and then studying 
the beam cross-section as the particles emerge from the lens. At first 
two crossed thin wires were used to shape the injected beam cross-section. 
It was found that these wires caused scattering of some of the particles 
which gave these particles a large transverse momentum, which in turn 
caused the injection conditions of these particles to lie outside the 
acceptance limits of the structure. It was then decided to use a 
rectangular aperture to give the injected beam a rectangular cross-section, 
similar to that discussed in Section (7.4). The detected exit beam cross- 
section can then be compared with the results obtained in Section 7.4. 

The circular disc E, shown in figure (8.11), has a rectangular 
aperture whose sides, equal to ne 1.56 mm and No = 0.85 mm, are parallel 
to the fixed X, Y axes. In order to detect the degree of rotation of 
the emerging beam without ambiguity, the sides of the aperture are chosen 
non-equal. The ratio between their lengths is approximately 1] : 6, the 
longer side being in the fixed X-direction, as shown in figure (8.12). 

To ensure that the tip of the charging electrode can be considered 
aS a point source and hence to make sure that the cross-section of the 
injected beam at the input plane, Z = 0, is a rectangle; the following 
experiments were performed. 

The first experiment was to examine the microparticle distribution 
at the plane of the aperture E to ensure full coverage of this rectangular 
aperture. The experimental arrangement shown in figure (8.13a) was used. 


The rectangular aperture was replaced by a greased glass slide. The 
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Figure (8.12) The Injected Beam Shaping Aperture. 
The rectangular aperture E is mounted in 
front of the lens injection plane to shape 
the injected beam cross-section. 
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Figure (8.13a) Sketch of the Experimental Arrangement Used to 
Test the Injected Beam Cross-Section. 


The greased detector plate is first placed at 
the plane of the aperture E and then at the 
injection plane of the lens beyond the aperture 
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Figure (8.13b) Photograph of Injected Microparticle Beam Cross- 
Section at the Plane of the Shaping Aperture E. 


The detector slide was moved horizontally to 
photograph the detected beam cross-section and 
the aperture E simultaneously. The illuminated 
circular area behind the rectangular aperture is 
the reflection of light from the tapered hole in 
perspex plate supporting the aperture E. 


Figure (8.13c) Photograph of the Injected Beam Cross-Section 
at the Lens Injection Plane. 


The detector slide is still mounted at the lens 
injection plane. 
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greased detector slide was mounted on the disc of the aperture E. 

Figure (8.13b) shows the detected beam cross-section. The detector 

slide was moved horizontally to be able to photograph the detected 
cross-section. Figure (8.13b) also shows the rectangular aperture E 

and shows that the detected beam cross-section is reasonably uniform. 

The diameter of the circular aperture D was 1 mm and the distance 

between it and the tip of the charging electrode was 28 mm. The distance 
between the tip of the charging electrode and the surface of the greased 
glass slide, the position of the aperture E, was 67.5 mm. The diameter 
of the detected beam cross-section was equal to 2.4 mm which is equal to 


] mm. ee . This result indicates that the tip of the charging electrode 


can be considered as a point source. Also the detected beam cross-section 
at the position of the rectangular aperture ensures full coverage of this 
rectangular aperture by the microparticle beam. 

A second experiment was performed to insure that the microparticle 
beam cross-section, at Z = 0, in the presence of the aperture E between 
the charger and the lens as shown in figure (8.11), is a rectangle. A 
greased glass slide was placed at the input plane, Z = 0. With no 
quadrupole voltage and using a charger voltage equal to 15 KV, the 
detected beam cross-section was, as shown in figure (8.13c) in which the 
detector plate is still mounted on the entry of the lens, a rectangle. 
This rectangle had its sides parallel to the electrode axes at Z= 0, 
that is parallel to the sides of the aperture E. The microparticle 
distribution over the detected rectangular cross-section was more or less 


uniform. The sides of the rectangle were equal to 1.73 mm and 0.95 mm, 
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the longer side being in the fixed X- direction. 

The last two experiments show that the tip of the charging 
electrode can be considered as a point source at a distance (d, ce d.) 
from the lens input plane. Therefore the experimental arrangement of 
figure (8.11) is exactly that of figure (7.5a). The injected beam at 
Z = 0 has a rectangular cross-section centered at the center of the 


quadrupole lens aperture. The sides of the injection rectangle are 


parallel to the fixed X and Y axes and are equal to 


1.73 mm x 


SS 
—J 
QO. 

iH] 


0.95 mm : Ce 


where ny and Ny are the dimensions of the aperture E and (d, .d,) are 
as shown in figure (8.11). 

Also by considering the tip of the charging electrode as a point 
source, the transverse momenta of each particle at Z = 0 are related to 


the coordinates of the particle at Z = 0 as, 


Z 
ri d a Xo ; 
0 1 2 
ors 
Py a Fd, ‘o Se 


where P5 is the axial momentum of the particle at injection. The 


0 
dimensions and the position of the aperture E during the experiments were 
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chosen such that the angle of divergence of the injected beam is small. 
Thus, the axial momentum of each particle at injection can be assumed 


equal to its total injection momentum are Hence 


Ps = ee = (2qmy)'/2, where V is the charger high 
0 


voltage = 15 KV. (8.3) 


The analysis of Section 7.4 shows that for an injected beam 
having the rectangular cross-section provided by the aperture E and whose 
particles have their coordinates and momenta related by equations (8.2) 
and (8.3), the exit beam cross-section at the exact imaging values of 
the lens parameters should be a rectangle. This exit rectangle should 
be at the center of the quadrupole aperture and its sides should be 
parallel to the electrodes axes at the exit, the (X7,Y7) axes. The beam 
cross-section should be rotated 27 + 43° with respect to the fixed X, Y 
axes in case of the counter-clockwise twisted lens. The exit rectangular 
cross-section should have the same dimensions as the cross-section of 
the rectangular injected beam at Z = 0. Moreover, the mechanical 


momenta of any of the exit particles in the fixed target plane should be 
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where X- and Yr are the exit coordinates of the emerging particle with 
respect to the (Xr Yrs Z) fixed target coordinate frame. 

According to the relations of (8.4), if a drift space d, is 
introduced beyond the lens exit, the particles should cover an expanded 
rectangular area whose sides are still rotated 43° with respect to the 
X, Y axes. The ratio between the lengths of these sides should be that 
between the sides of the rectangular aperture E, that is npn = eGr 

In the experiments, the microparticle beam cross-section was 
detected on a greased glass slide as the beam emerged from the lens. 
Arrangements were made so that the distance between the lens exit plane 
and the greased face of the glass slide could be adjusted. Also a 
ratchet and vacuum feed-through mechanism were used to move the greased 
glass slide in the horizontal direction parallel to the exit plane so 
that more than one experiment could be conducted without opening the 


experimental vacuum chamber to air. 


8.4.1 Experimental Results 


(a) Single Twisted Electrostatic Quadrupole Lens 


The experiment which is described below comprises two separate 
experiments, one for the counter-clockwise twisted lens, the other for 
the clockwise twisted lens. For all experiments the greased detector slide 
was mounted exactly at the end of each of the single twisted lenses and 
the charger high voltage was constant at 15 KV. For the initial part of 


the experiment the quadrupole focusing voltage was adjusted to the exact 
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theoretical imaging value of 1200 volts. The detected exit beam 
cross-section at this focusing voltage closely approximated a rectangle 
but the ratio between its sides was less than 11:6; the longer sides were 
smaller than expected from equation (8.1). The angle of rotation of 

the sides of this cross-section with respect to the fixed (X,Y) axes 

(the horizontal and vertical axes) was slightly smaller than 27 + 43°. 
The theoretical study, already described in Section 7.4 was then carried 
out and it became apparent that the twisted lens was in effect slightly 
detuned from its actual imaging value. 

The previously mentioned ratchet mechanism was then introduced to 
allow movement of the detector plate in the horizontal direction so that 
several experiments could be carried out without letting the system up 
to air. A series of experiments was then carried out while the focusing 
voltage was varied in steps from 1150 volts up to 1250 volts. It was 
found that a potential of 1210 volts produced the best image for both the 
counter-clockwise and the clockwise twisted lenses. At this value of 153 
the detected beam cross-section is a rectangle whose sides are parallel 
to the electrode axes at the exit (X7¥ ), that is rotated 27 + 43° with 
respect to the fixed (X, Y) axes. The dimensions of the rectangle are 
equal to 1.73 mm and 0.95 mm, exactly equal to those of the injected beam 
at Z = 0. Photos of the detected exit beam cross-section at us = 1210 
volts, with the detector plate mounted at the exit of the lens, are 
shown in figure (8.14) for the counter-clockwise twisted lens and in 
figure (8.15) for the clockwise twisted lens. 

The detected exit beam cross-sections with no drift space beyond 


the exit of the lens and at values of the focusing voltage less than, 
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Figure (8.14) Beam Cross-Section at the Exit Plane of the 
Counter-Clockwise, 8 + ve, Twisted Lens at 
ue = 1210 volts. 


This is a rear view with the detector slide 
still mounted at the exit plane of the lens. 


Figure (8.15) As in Figure (8.14) for the Clockwise, 8 - ve, 
Twisted Lens. 
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Beam Cross-Sections at the Exit Plane of the 
Counter-Clockwise Twisted Lens for Different 
Values of Vo: 


The values of Vo are (from left to right) 1190, 
1200, 1210 and 1220 volts. The cross-sections 
are photographed in the direction of incidence 
of the beam on the detector. 


Scale of reproduction: 1 cm on photograph > .21cm. 


Figure (8.17) 


Beam Cross-Sections at the Exit Plane of the 
Clockwise Twisted Lens for Different Values of 


Ve 
fo) 


The values of Vo are (from left to right) 1180, 
1190; 1200; 1210 and. 1220 volts se thencrass= 
sections are photographed in the direction of 
incidence of the beam on the detector. 


Scale of reproduction: 1 cm on photograph > .22 cm. 


Ors , 


ait Yo enel4 3143 dF 36 engtsase-redyd ane? (31,8) etuptd 
tpaeettid aot enst-bodetwl setwigol}-redaued 
P V ra. 2a) OV 


eff (gdeta of Ts! matt) o6 ov Yo saul off 

any? }gee-ca0ve sar. .2ezfov OSS! bre OTST ,Q0Si 

antah boi? Yo aetsoarlh ald nf bernqaevgos ong 61s 
Aotosseb add no med SAI Yo 7 


eors. « dadypotorq vo m | snolJoubouga) To Sinse 


ant to: sett ttxd ode Je 2e0teast-e20) miod (51.8) owptt 
Yo xeuls¥ semrvett!O Wo? eval betetwT eetwaootl 


"9 
Oat (riety of #teT wot) eis VY Yo 2oufny ont 
<8em oft .zitov OSS! bas OTS ,OOST Ott 


Yo Notdogri ly 4d7 of bedqegolodq as. enotiner 
Aetoaryh oft ao mesd urls Yo sansbdiont 


+i gs. + fasypodorg no nd f snolisvborqe Fo afpo2 


7 - 


raat | 


equal to and greater than 1210 volts are shown in figure (8.16) for the 
counter-clockwise twisted lens and in figure (8.17) for the clockwise 
twisted lens. For values of ve greater than 1210 volts, the detected 
beam cross-sections are nearly rectangular parallelograms whose sides 
are rotated with respect to the X, Y axes through an angle slightly 
larger than 2m + 43°. The angle of rotation as well as the dimensions 
of the parallelograms increase as the value of the focusing voltage 
increases. Also the lengths of the sides that should be parallel to 
the X- axis at ve = 1210 volts, are observed to increase at a higher 
rate than that of the increase of the other two sides. 

As the focusing voltage is decreased to values less than 1210 volts, 
the detected beam cross-sections are nearly rectangular parallelograms 
rotated with respect to the X and Y axes through an angle less than 
27 + 43°. The angle of rotation decreases as the focusing voltage is 
decreased. The parallelograms detected at values of Ne less than 1210 volts 
have smaller areas than the eecencuir cross-section detected at 1210 
volts and it is found that the exit beam cross-section decreases as Ve is 
decreased. It is also observed that the lengths of the sides that should 
be parallel to the X+ axis at 1210 volts decrease at a rate higher than 
that of the decrease of the other two sides. The experimentally observed 
phenomena are in good agreement with the theoretical results of Section 
7.4. However, it is clear from figures (8.16) and (8.17) that the density 
of the particles is higher in the middle of the exit beam cross-section 
than on the boundaries. This exit beam density distribution is different 


from the more or less even microparticle distribution over the rectangular 
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Figure (8.18) 


Beam Cross-Sections at a Distance d3 = 7 mm 
Beyond the Exit Plane of the Clockwise Twisted 
Lens for Different Values of Vee 


The values of Vg are (from left to right) 1150, 1180, 
1200, 1210 and 1220 volts. Direction of photography 
the same as figure (8.17). The bubbles in the 
background are due to some moisture trapped under 

the grease and heated during photography. 


Scale of reproduction: 1 cm on photograph > .25 cm. 


Figure (8.19) 


Beam Cross-Sections at a Distance d, = 17.5 mm 
Beyond the Exit Plane of the Clockwise Twisted 
Lens for Different Values of Vo: 


The values of Vo (from left to right) 1150, 1180, 
1200 and 1210 volts. Direction of photography 
the same as figure (8.17). The overlapping of 
the second and third exit beam cross-sections was 
due to insufficient movement of the ratchet 
mechanism. The effect of trapped moisture is 
more pronounced in this case. 


Scale of Reproduction: 1 cm on photograph + .27 cm. 
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Figure (8.20) 


Beam Cross-Sections at a Distance ds = 33 mm 
Beyond the Exit Plane of the Clockwise Twisted 
Lens for Different Values of ie 


The values of V, are (from left to right) 1150, 
Msco-st200) 1210, 1220 and 1230 volts. Direction 
of photography the same as figure (8.17). 


Scale of Reproduction: 1 cm on photograph > .31 cm. 
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beam cross-section at the injection plane. A discussion of this form 
of aberration will follow the presentation of the remaining experimental 
results. 

While the foregoing series of experiments examined the imaging 
of the particle coordinates, the series of experiments which is described 
below examines the imaging of the particle momenta. Since the two lens 
sections used in the previous series of experiments are identical except 
for the direction of twist of the electrodes, it was decided to carry 
out these experiments for the counter-clockwise twisted lens only. The 
experiments were performed by using different values of drift space d3, 
as shown in figure (8.11), between the exit plane of the lens and the 
greased detector slide. 

Figures (8.18), (8.19) and (8.20) show the exit beam cross-sections 
at drift spaces equal to 7 mm, 17.5 mm and 33.6 mm respectively. At 
Up = 1210 volts, it is clear that the exit beam retains its rectangular 
shaped cross-section and the ratio between the lengths of the rectangle's 
sides is still 11:6. Also the directions of the sides are still parallel 


to the X, and Yr axes and the dimensions of the rectangular cross-sections 


t 
at Ue = 1210 volts in all three cases are exactly as theoretically 
predicted in Section 7.4. In interpreting figures (8.18), (8.19) and (8.20) 


it must be noted that their scales of reproduction are not the same. 


(b) Two Oppositely Twisted Equal Sections 
The oppositely twisted sections that were investigated separately 
were mounted and electrically connected so as to have no drift space 


between then, as shown in figure (8.21). The rectangular aperture E 


ays 
avo? ett a noteeusetb A ,snalq notaoatnt ond 36 
fesvombveqs gatatsmay ott to noltednagerq sr 


ortpaut sit benhnexs atnemtysqKe to soiree antoperet iat 
hatiyoeah at dobitw atnoumtvsqxs to estuee srt .2edentinoos ototng 9 
eval ows og eont2 .ssesmom ofotyreq ‘sre ee os aaiitenh: wor 
Sqenas (sotswabt ote edngmiisqxa *o entree: a | 
vries of beltosh zaw t) , esboyvtaela ert ie said 
eiT .\fno ane! batefwi 22 twioselo-iedmbos easier pene er 
gh sneqe F7Pab Te epufuv prsvtTih pateu uh bs | aq or0N eineuhrogns 7 
ait brs anal 9d Fo snelq dix oft nvawied Ut. 8) matt nf ron To 


— 
_sbtle ‘osostab bovsore 


: 


7 
enottose-22079 weed Jtxo ond wore (09.8) bns (@f.8) . (6.6) cowort i _ 


th .ylevtdonqz2e7 mu 3.€€ bos mm 2.T7 me OF fsups 2s6qe sib te 
wa(upemtony ett 2vtetey msed ttxs sdt tera wola at ti .277ov orsh + v 
z'algnstooy odd to attgrel add poswrsd olfiy edt baw Notdase-220%9 begat 
loffavaq (fht2 eve zebte edt Yo enotiasitb aid OFTA Or Tf Tie 2t aabie 
enotios2-220%9 Yatuynatoet ait %) ernvtenom’o ay bmg tone -Y bas ph afd of: 
vl leoteeroedt 2p ylfoexs sie zeeno aerid (fe of edfov Orst = 4V | V ts 7 

(0$.8) bas (@f.8) . G08) eowort onttevqrvetnt ol 5b,1 notiose Ht aoe : 
ome2 oly Jon ws noltouborgey to esface Fond tarts beton = 


2vorioa2 [sup betetwr leat: owt (a) 


\ lave tagse heseptizavat avow Jel) enotiase ey tesVeeu ar ae 
e96qg2 Si1b on ever oF 25 2 betoenNGd cake 7 9 bas be: ; re 


q ad ysTugnsise1 sit .(/S.8) svupt?. nt 


yy ”* : 
ae. ae : 


Z\o 


Figure (8.21) Photograph of a Lens Structure that Consists of 
Two Adjoining Sections of Equal Length but 
Opposite Twist. 


Photograph shows the mounting base plate and the 
focusing voltage calbes. 


Figure (8.22) Beam Cross-Sections at the Exit Plane of the Lens 
Shown in Figure (8.21) for Different Values of V). 


The values of Vg are (from left to right) 1190, 
1200, 1270912208 and 1230) Volts: "Direct ionvoT 
photography the same as figure (8.17). 


Scale of Reproduction: 1 cm on photograph > .21 cm. 
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previously used was again employed in these experiments. As before the 
Sides of the rectangular aperture were kept parallel to the axes of the 
electrodes of the first twisted Tens at Z = 0, that 7s the sides were 
parallel to the horizontal and the vertical directions. According to 
the theoretical analysis of section (7.3), the exit beam cross-section 
at the exact imaging voltage should be a rectangle with its sides in the 
horizontal and vertical directions. The dimensions of this rectangle 
should be equal to those of the injected beam rectangular cross-section 
at Z = 0, that is 1.73 mm in the horizontal direction and 0.95 mm in the 
vertical direction. The two section lens structure was tested at 
different values of lie and the corresponding beam cross-sections, detected 
right at the end of the structure, are shown in figure (8.22). This 
experiment verified that exact imaging with no rotation with respect to the 
fixed X, Y axes occurs at V7 1210 volts. The beam cross-sections at values 
of be greater than 1210 volts are slightly expanded but their shapes are 
still nearly rectangular parallelograms. These parallelograms are 
Slightly rotated with respect to the X, Y axes. 

The results of the experimental investigations of either the 
single twisted lens or the structure consisting of two oppositely 
twisted sections show that these lenses are capable of identically imaging 
the injection conditions of the particles of the microparticle beam at 
ie = 1210 volts. Also the results shown in figures (8.16) to (8.20) 
and in figure (8.22) are in good agreement with the theoretical analysis 
of Chapter 7. However, these figures show that there is some aberration 


that results in a higher microparticle's density at the middle of the 
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exit beam cross-section than on the boundaries and in some slight 
curvature of these boundaries. 

This aberration is undoubtedly due to several factors, the 
most important of which is the misalignment of the axes of the micro- 
particle injector and the lens structure. In some preliminary 
experiments the microparticle charger and the lens structure were 
separately mounted and aligned optically. The exit beam cross-sections 
were so distorted as to be almost unrecognizable. Subsequently the 
charger and the lens were mechanically mounted to one rigid base plate 
and then aligned optically. With this arrangement it was immediately 
possible to obtain exit beam cross-sections that clearly demonstrated 
the imaging properties of the lens. The effect of misalignment was 
found to be most pronounced in the case where the structure consisted of 
two oppositely twisted sections. In this case proper alignment was 
most difficult to achieve. This sensitivity to misalignment is common 
to all quadrupole et riletures to 

During the experiments the quadrupole focusing voltage was 
kept constant within a few volts by regulating the mains feeding the 
quadrupole voltage supply. However, the mains feeding the charger high 
voltage source was not regulated and some of the observed aberration 
may be due to drift in the value of V during the experiments. 

In Chapter 5 it was found that for V5 = 1200 volts, the 
effect of gravitational forces is negligible 1f it >> 0.0003 coulomb/kilogram. 
During the charge detection experiments described in Section 8.2.2 it was 
observed that the values of the charge-to-mass ratios of the particles 


fluctuate over a wide range. Some of the particles emitted from the 
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charger may well have charge-to-mass ratios low enough that their 
motion will be noticeably affected by the gravitational forces. Such 
particles will not be properly imaged and hence contribute to the 
observed aberration. 

Further, fringing fields at the entry and exit planes of 
the lens structure as well as any tolerances in the length, radius, 
mounting and cutting of the electrodes may be responsible for some of 
the aforementioned aberration. 

The aberrations discussed above are significant only in that 
they affect the identical imaging properties of a lens section. They 
are unimportant in terms of simply guiding and transporting a beam of 
particles. The experimental results clearly indicate that in all 
cases the injected particles are confined within the aperture of the 
structure, inspite of any aberration. In the experiments described 
above,particles have been successfully guided over a maximum distance 
of approximately 67 cm. It is important to note that without the 
twisted structure these same particles would be spread over an aperture 
approximately 19.5 mm in diameter as opposed to the 1.73 mm x 0.95 mm 


rectangle produced by the lens at re = 1210 volts. 
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CHAPTER 9 
THE CIRCULAR UNIFORMLY TWISTED ELECTROSTATIC QUADRUPOLE 
9.1 Introduction 


As a beam transport system is often required to bend the beam 
to direct it to a certain experimental area, it is of interest to analyze 
the effect on the particle motion caused by gradually bending the twisted 
quadrupole structure. Such effect is studied in this chapter by analyzing 
the motion of a charged particle which travels along a circular uniformly 
twisted electrostatic quadrupole structure. Beside being of importance for 
beam bending devices, the above analysis is also relevant to the design of 
storage rings for heavy charged particles. 

Figure (9.1) illustrates the coordinate systems used in the 
analysis that follows. The (X,Y,Z) coordinate system is a fixed coordinate 
system whose origin coincides with the circular axis of the bent structure, 
Which has a radius Ro? at ¢ = 0. The coordinates (X', Y', o) are related 


to the (X,Y,Z) coordinates as follows, 


k=" Cos. ht RC - cos 9) 
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9.2 The Hamiltonian of the Motion 


The kinetic energy of a particle of mass M in terms of the 


(X', Y', #) coordinate system is 


T= OR + 2 + 27] 
see a eaten 1 (9.2) 


If the radius of curvature of the circular axis of the structure, 
Ro? is much larger than the aperture radius, a, and if (ga)? << 1, then the 


potential at any point within the aperture in terms of the coordinates X', 


Y' and ¢ is 
V(X',Y'.6) = KLLX'? - ¥'43c0s 2eR + {2X'Y'}sin 2eR,6]. 


Since V(X',Y',¢) is velocity independent, the momenta conjugate to the 


(X', Y', ¢) coordinates ec 
we Nana re 
eae ae 
page Duck pc 8 
Pv sy" MY 
P, at Mb(R, - xt )* (9.3) 


and the Hamiltonian of motion in the (X', Y', $¢) coordinate system 
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Equation (9.4) can be used to derive the equations of motion in the 

(X', Y', >) coordinate system. However, the equations of motion thus 
derived are found to be nonlinear and coupled. At this stage the problem 
can be simplified by referring the motion to the rotating coordinate system 


(x', y', ¢) where coordinates (x',y') are related to X', Y' by 


X="X" COs BRo¢ =4V." 1S 0M BRO 


eae ex sin BRoo + y' cos BRoe (9.5) 
Differentiating equations (9.5) and substituting into equation (9.2), the 


kinetic energy of the particle in terms of the coordinate system (x', y' ¢) 


is 


+ 2eR A(y! x! KE yt) + BORE §A(xI® + y!*)] (9.6) 


and the momenta conjugate to x', y', ¢ are 
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Using equations (9.3), (9.4), (9.5) and (9.7) it now follows that the 


Hamiltonian in the (x', y', ¢) coordinate system is 
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then to first order the above Hamiltonian becomes 
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(9.10) 


The Hamiltonian of equation (9.10) leads to coupled equations of motion 


and a new coordinate system must be generated in which the motion is 
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uncoupled. It is found that the generating function 
=x po eyo . BP ais 
Cee et Y Rye coup my os (9.11) 
generates through contact transformation the coordinate system (x', y', ¢, 


yrs Pyrs Py) in which the equations of motion are uncoupled, provided that 


the constant p is chosen properly. The new coordinates are related to the 


old coordinates py (34) 
aG - 
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) ek 9 
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p. = LE x : (9.12) 
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Upon substitution of (9.12) into (9.10), the final Hamiltonian is obtained 


as 
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9.3 The Equations of Motion 


The following approximations are now introduced: 


ew Its assumed ithat Py = Po: 


and (9.12) if first order quantities are neglected. Then from equations 


This follows from equations (9.7) 
(9.3), to the same degree of accuracy 


p= MRS G (9.14) 
2. Based on the results for the straight uniformly twisted 
electrostatic quadrupole, where by neglecting first order quantities it 
was shown that the axial velocity Z is essentially constant, it is now 
assumed that the circumferential velocity, Roos is also constant. Thus 


from equation (9.14) 


Py Sage ard. (9.15) 


where Uo is the initial circumferential velocity. 
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The value of the constant p which appears in the generating 


function G is now selected as 
p = M Un = Po (9.16) 


Finally at this stage, the equations of motion are derived from 
equation (9.13). These equations are coupled. However, upon introducing 
equations (9.15) and (9.16), the equations of motion become uncoupled and 


appear as follows: 
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where the stability factor Sy is given by 


Ze 
BP 


So ~ 2Mqk (9.18) 


It now follows immediately from equations (9.17) that 
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y' + oak (5 ee Wid ae Celts 25,5) $ sin gU,t (9.19) 
since, from equations (9.12), (9.14) and (9.15) $ = (U,/R,)t. 

One recognizes immediately that the homogeneous solutions of 
equations (9.19) are exactly analogous to the solutions obtained in 
Chapter 2 for x and y. The complete solutions of equations (9.19) are 
the sum of these homogeneous solutions and the particular solutions, 


where the latter are 


Py (1 + 259) 
x' = - ———_——. cos BUt : 
2Mqk i 


po (1 - 285) 


ae ih BU, 
2Mqk Ro 


When these particular solutions are transformed back to the X', Y', ¢ 


coordinate system, X' and Y' of these solutions are 
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S 
a = (sin 26R, 4) (9.20) 
B Ro 


The projection of the particular solutions of equations (9.20) 
onto the transverse plane X'-Y' is plotted in figure (9.2). This figure 
indicates that the particle trajectories in the circular twisted structure 
are in fact essentially the same as the trajectories in the linear structure, 


except that the mean orbit of particles no longer lies on the structure axis. 
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Instead, the mean orbit is now a helical orbit whose axis is displaced 


452 
radially outward in the X-direction by an amount = . Ihe radius of 
é S gk 
this helical orbit-is equal’ to >= . . 
BR 
0 


direction of 
twist 


Figure (9.2) The Mean Orbit of Particles Along the 
Circular Twisted Structure. 


Because of the displaced mean orbit the acceptance conditions for the 
circular structure are difficult to compute and no expression for the 
confinement criterion is obtained. However, it is of interest to study at 
least one specific example to illustrate what order of magnitude the radius 


of curvature Ro must have relative to the aperture radius a. 
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Example: Consider a particle which is injected with initial conditions 


5 eae , Vos-—32 (9.21) 


As mentioned above the complete solutions of equations (9.19) are the sum 
of the homogeneous solutions, analogous to those given by equations (2.31) 
and the particular solutions. It can be shown that the trajectory of the 
particle whose injection conditions are given by equations (9.21) is just 
that of the particular solutions of equations (9.20) and the homogeneous 
solutions for these particles are zero. 

If the maximum radial displacement of this particle is now 


chosen as a, then 


A (Lt 4s) (9.25) 
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it then follows from (9.25) that Ro = 180a and that the circular twisted 
quadrupole must execute a total of 57 twists in going once around the 


circular structure. 
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CONCLUSION 


During the past few years special attention has been focused 
on the charging and acceleration of heavy ions as well as heavy 
micron-sized solid and liquid particles. The capability of the 
uniformly twisted electrostatic quadrupole structure to guide beams 
of heavy microparticles has been investigated in detail. A first order 
analytic solution for particle trajectories has been developed. This 
analytic solution shows that stable trajectories exist for particles 
with arbitrary charge-to-mass ratios provided that the structure 
parameters are properly chosen. It was shown that it is possible to 
confine particles stably provided their injection conditions are within 
the acceptance limits set by an analytically derived confinement 
criterion. This criterion was used to study the acceptance limits of 
the twisted structure for some special cases of injection. Comparison 
of the analytically computed trajectories with those computed through 
numerical integration of the exact equations of motion showed that the 
analytic solution is accurate over the entire quadrupole aperture. 
Furthermore the comparison showed that the confinement criterion 
successfully predicts the acceptance limits of the twisted structure. 
The effect of gravitational forces on the motion of heavy particles 
along the twisted structure also has been studied. It was found that 


the influence of gravity is negligible provided the charge-to-mass 
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ratio of the particle is higher than a value that depends on the 
structure parameters. 

The guiding properties of the twisted structure were compared 
with those of the classical electrostatic quadrupole channel. It was 
found that the twisted structure has a stronger focusing action and 
hence is especially suitable for guiding beams emerging from point 
sources. 

The investigation of the analytic solution showed that for 
Suitably chosen parameters a finite length of the twisted structure 
can be operated as a novel type of lens. By using this imaging lens, 
the initial injection coordinates and momenta of a system of particles 
were found to be identically reproduced in a coordinate plane located 
at the exit of the lens. The imaging takes place with a small amount 
of aberration. Numerical investigations showed that this aberration 
can be eliminated almost entirely by constructing the lens of two equal 
but oppositely twisted sections. 

The identical imaging property of the twisted quadrupole has been 
experimentally tested. The lens structures have been fabricated by 
moulding. A microparticle contact charger was designed to provide the 
high density microparticle beam that was used in the experiments. The 
results of the experiments were in good agreement with the theoretical 
analysis inspite of the presence of some aberration. It is expected 
that this aberration can be reduced if more accurate methods of 
alignment and machining are used. 

Also, particle motion has been studied along a twisted structure 


whose axis is bent in a circle. The results of this study are important 
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to the design of beam bending devices and to the design of storage 
rings for heavy microparticles. The analysis showed that for proper 
choice of the structure parameters stable trajectories exist and that 
the mean particle orbits are displaced from the structure axis. 

The twisted electrostatic quadrupole structures can be used as 
beam transport elements in conjunction with very heavy ion accelerators, 
or with accelerators for simulation of micrometeoroids. A further use 
could be the guiding of electrostatically sprayed liquids, such as the 
exhaust streams of colloidal propulsion devices. Another application 
could be as a confinement device for clouds of positively and negatively 


charged microparticles for the study of microparticle plasmas. 
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APPENDIX A 


POTENTIAL DISTRIBUTION IN THE UNIFORMLY TWISTED ELECTROSTATIC QUADRUPOLE. 


The potential at any point in a charge free region should satisfy 


Laplace's equation, which in polar coordinates is 


2 Z 2 
a V 1 ov now \. a V 

+ — — + — ——~ + ——~ = 0 
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For an infinitely long, uniformly twisted structure,the potential at a 
point (r, 6, Z) is exactly that at the points (r, 6 + g2, Z + 2%) where 
2 can be any distance. Therefore the potential should be in the form 


V(r, 6 - BZ) = R.(r)#. (6 ---BZ) 
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Substitution of this potential in Laplace's equation, and separation of 


the variables gives 
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where I (ner) is the modified Bessel function of the first kind and 


order n. 


For a symmetrical twisted quadrupole structure the following 


conditions must be satisfied (see figure (A.1)) 


V(r, @ - BZ) = V(r, -6 + @Z) 


ul 


V(r, 6 - BZ) = V(r, 6 - BZ + 7) 


Vie. e0= 67) 


-V(r, 6 - (Les) ; 
These conditions imply that V(r, 6 - 8Z) should be 


V(r, @ - BZ) = j I,,,(2mBr )cos 2m(e@ - BZ) (A.1) 


in 
m=1,3,5 


Equation (A.1) gives the potential produced by a twisted structure whose 
electrodes can be of any shape and twisted at any rate. The coefficients, 
On? must be numerically evaluated by matching the potential at points on 
the electrode surfaces with V(r, 6 - 8Z) of equation (A.1). This matching 
produces a set of independent algebraic equations whose number is equal 
to that of the points considered. 

For the special case of slowly twisted structures the Bessel 


functions of equation (A.1) can be approximated as follows: 
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For guiding and focusing structures the only values of r which are 
of interest are those within the aperture, r sa. (hus, the! potential 
at points within the aperture of a slowly twisted quadrupole structure, 
with any electrode shape, is 


V ‘2 kr? cos 2(6 - BZ) + kor° cos 6(6 - BZ) 
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provided only that (ga)? << 1. The coefficients Ky can be evaluated by 
potential matching at points on the electrode surfaces. For the special 
case of hyperbolic shaped electrodes, the surfaces of the electrodes are 


described by 


(a) Positive electrodes, V = Vo 
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ry cos 2(6 =8Z) =a 


(b) Negative electrodes, V = -V 
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Therefore, for the slowly twisted structure with hyperbolic electrodes, 


the coefficients of equation (A.2) are 


V 


K = ae k = 0 foram = 1. 
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The potential at points within the aperture when (ga)¢ << ] is thus given 


approximately by 


V 
V=Kr* cos 2(6 - BZ) where k = + (A.3) 
a 
Since 
er 
and e = tan” : 


equation (A.3) may be rewritten in terms of the fixed coordinates (X,Y,Z) 


as: 
V = kL(x? - ¥*)cos 282 + (2XY)sin 262] (A.4) 
Alternatively 
ne 
and 6 = e7 = ne + 


and in terms of the rotating coordinates (x,y,z) 


Pe kOe 5) (A.5) 


Equation (A.5) shows that the potential of the slowly twisted 
electrostatic quadrupole, when viewed in the rotating coordinate system, 


is negligibly different from the potential distribution in a classical 


quadrupole. 
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APPENDIX B 
THE AXIAL VELOCITY Z 


The objective of this appendix is to obtain an expression for the 
axial velocity of the particle, in terms of the transit time, the initial 
axial velocity, and the values of the transverse displacements and 
momenta at injection. It will also be shown that the axial velocity of a 
particle in the structure varies only slightly from its initial value, 


Ze 
The axial velocity Z is given by equations (2.32) and (2.33) as, 
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Examination of equation (B.2) shows that the following upper 


bound may be placed on the magnitude of the quantity g(x Py ay, PL): 
| ap Ce BC, CoM 
| 8 (x Py ~ ¥ P,) | >. 2 7 Co 2f | 
» 65 : ae a 
2 C.Y2qkM 
= |62p eh ales 
ig ney or aoe 
2 
etgaiepee ta cece a 
Ze vera Walecr 


Therefore 


[C., | 
[B(x Py ¥ Px) Imax < 8 Py Le | one eeu T 15! dinax 


(B.4) 


It is shown in Chapter 3 that for a particle to be confined within the 


aperture, 
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Substituting this expression in equation (B.8) gives 


ees E 
[ P Ace Akl 
ey) 2 
o[]l-—(X P, -Y. P, )] 
peooe Yo 0) Xo 


Using the inequalities (B.6) and (B.7) the axial velocity is found to be 
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APPENDIX C 


DERIVATION OF THE MATRIX FORM OF THE ANALYTIC SOLUTION 
Equation (2.46) gives the x coordinate as 
aie gree 
Ay ieee aed 


Substituting for A, Ao A. and A, from equations (2.48) gives 
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2qkM {cos fat SC Saar t}Py (C.1) 
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Similarly, 
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Equations (C.1), (C.2), (C.3) and (C.4) can be combined in a matrix form 
giving equations (2.49) and (2.50). 

Further consideration of equations (C.1) to (C.4) shows that they 
can be rearranged and the trigonometric expressions can be combined to 
obtain the form of the trajectory solution given by equations (2.31). 


Equations (C.1) to (C.4) can be rearranged to take the form 
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Comparing equations (C.6) and (C.7) with equations (2.31) and (2.37) 
it is clear that equations (C.6) and equations (2.31) are the same and 
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APPENDIX D 


COMPUTATION OF ACCEPTANCE FOR POINT SOURCE INJECTION 


Equation (3.13) sets the limits for the transverse momenta at 


injection, from a point source at Z = 0, as 
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This equation is more easily interpreted when expressed in terms of the 


polar coordinates: 
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This equation was used to plot the transverse momenta acceptance diagrams 
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shown in figure (3.5). These curves intersect the normalized <2 and 
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Solving equation (D.3) numerically gives S, = 1.3295. Substitution of 
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this value of Sy into the expression for ls—| gives 
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APPENDIX E 


COMPUTER PROGRAMS LISTINGS 


The first program in this Appendix was used for plotting the 
two sets of trajectories that were used for the evaluation of the 
analytic solution developed in Chapter 2. To obtain the first set 
of trajectories, the analytic trajectories, equations (2.31), (2.34) 
and (2.35) are used to compute the trajectory points in the rotating 
(x,y,z) coordinate system. Then equations (2.7) and (2.13) are used 
to refer these trajectories to the fixed (X,Y,Z) coordinate system. 
The trajectories are then plotted in this fixed (X,Y,Z) coordinate 
system. The second set of trajectories are computed directly in terms 
of the fixed (X,¥,Z) coordinate system by numerically integrating 
equation (2.6) using the Runge-Kutta method for numerical integration. 

The Runge-Kutta method is a self starting method in which the 
calculation of the variables at the end of each interval depends on the 
values of the variables and their derivatives at the beginning of this 
interval. For all intervals the values of the variables are computed 
in the same manner by using for the initial values those at the 
beginning of each interval. The method may be best understood by 
studying Table E.1 which illustrates the steps for solving a second 
(36) 


order differential equation in onevariable x This method was 


easily applied to the system of differential equations, (2.6). 
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The second program in this Appendix computes the trajectory 
points of particles along a twisted electrostatic quadrupole lens 
Structure. The structure may be a Single twisted section or it may 
consist of two equal but oppositely twisted sections. The axial 
distance traversed by the particle was computed using the exact 
expression of equation (2.35). The value of s was determined for 


each particle individually in terms of the injection conditions. 
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FIRST PROGRAM 258 


PROGKAM TO COMPUTE AND PLOT PHE APPROXIMATE ANALYTIC TRAJECTORIES 
AND THE ACTUAL TKAJECTORIES (RUNGE- KUTTA TRAJECTORIES) OF 

A PARTICLE ALONG THE TWISTED ELECTROSTATIC QUADRUPOLE STRUCTURE 

THE ACTUAL TRAJECTORIES ARE COMPUTED THROUGH NUMERICAL INTEGRATION 
OF THE SECOND ORDER COUPLED DIFFERENTIAL EQUATIONS OF MOTION 
USING THE RUNGE - KUTTA METHOD OF NUMERICA INTEGRATION. 

DIMENSION AX (8002) ,AY (8002) ,AZ(8002) , DVAZ (8002) 

CALL PLOTS 

CREAT ORIGIN mds THE PLOTS 


CG) 


16 


a 


‘CHL IMELOT(0.0,150,-3) 


JIM= NUMBER. A PARTICLES WHOSE TRAJECTORIES ARE TO BE COMPUTED 
AND ELCTTED. ONE AT A TIME. 

READ (5,16) JIM 

FORMAT (12) 


= NUSBER OF POINTS COMPUTED ON EACH OF THE TRAJECTORIES. 


CHG) 


(@ 


VO= ACCELERATING POTENTIAL OF THE PARTICLES. 

AR= RADIUS OF THE APERTURE IN METERS. 

QM= CHARGE TO MASS RATIO (COULOMB PER KILOGRAM). 

SOO00= STABILITY FACTOR WHEN ( XO * PYO -YO * PXO )=0.0. 
REAL (5,1) J, VO,AR,QM,SO000 

FORMAT (I4, WE12. Vk 


PT=357015927 
ddd G=d 
AL=2.¥*PI*¥SQKT (10.0) *AR 
NI=250, 
REG, 
V=(2.*PI*AK/AL) * (2. *PI*AR/AL) * (VO/SOOOO) | 


VZO=SQRT (2. *OM*VO) 
AK=V/(AR¥AR) 
HO=AL/ (AJ*VZO) 

B= Sep /AL 

DO 555 MIN=1,J 

XO , YO , ZO THE INITIAL COORDINATES OF THE PARTICLE. 


DXO,DYO THE INITIAL NORMALIZED TRANSVERSE MOMENTA OF THE PARTICLE. 
DVZO THE NORMALIZED AXIAL INJECTION VELOCITY. 

READ {SP 2) X07 10,20 

READ (5; 2) 0x0, DYO,;DVZO 

FORMAT (3E12.6) 

XOO=X0 


YOOSYO ~~ 
Poo-20 

DXOO=DxXO 
DYOO=DYO 
S THE STABILITY FACTOR OF THE PARTICLE COMPUTED ACCORDING TO ITS 


INJECTION CONDITIONS. 


‘DELS=1. +B¥ (XO*¥DYO-YO*DXO) 


LEDS=1./DELS 
S=S0000*DELS* DELS 
FX=SQRT (2. ¥OM*AK*(S+1.)) 
FY=SQRT (2. ¥OM*AK* (S-1.)) 
ZOMNELO 


i> O'ER Go) 


XX11=X0 
si tl Ras PNG 
Trs0 210 
PVZ11=DVZO 
Dye DyO! =! 


ancrosns au ¢ a 
r Sas roeta ae , a8 


ate 
oe 


no 

aya DAITANNA ATE agaawe 

WHT SDYOEMUE TO MOWTEN APRA = SOUL f ver 

a siaiiaaiaile (5008) na OAL ae ee 
2T014 a 

_ Bross 7 7 fe sTTh , <5) 7 


= = =e Ge == = =< = 


at e050) 
cLrnsMOy Bu OF BHA CALHOTISEADT aie en sgu3hes vo ‘eanue. . bo 
A TA SHO .CSTTONY GHA OD 


lab ih phages, ” 


MSIGOTORGANT BH? TO tea BO ITE gps 10 aan = : 
HRIDL VA SUF AG TAT TH PrawZIsI5A =01 

-<eeUt@A 22 WavTAAGs GET TO meinak To 

e(NAMOOLIN 45% SMODNCD) OLTAR SAN OT BOMAND =8Q 

0.0e | ORD © OF- Od © OF ) HERW FOPAAI VTi11eaTe ae oe 
7 J os = 4 eee Efse.a3 yh Tr 
*“Tebetit kara 

Le0GLL 

dae (0.01) YAgeerte -SaIs 

«VeS=ta 


~ DALE 
(OOOCE\GV)* [FAV N497 94.5 (JA\GAOTE# AS) ——_ 
sage &) 1022059 
i (AAPRADNV =] 1A 
(OSV #04) \IA=08 
JAN IG* S28 - 
MIC, Perm eB? OG © 
-S292TH24 2NP to eSyanicwouy SAITIAL Bir OF » OY2, OX I ® 
foreak’ awry Va WeeA NDE FeLeTvER ALY GIST IAAROE JATPINT WAT OfG oro 
YRISGISY MOTTISGHT JAINA GESTIANGON Ant “OSV, © 
os sOl« #9 ere) aka 
oxida, oY, oxg (<, 
(0 -SFHE) 


7 - a oe - ay 


0 
G 
err Of Swiediaos ayry§ noe SJOLTHAY WET YO GOTUAD TYTITeATE 


7 dno> ROL TDS He 
rou ke is 204 


SVP Ea) Gala Silo ola: ‘ei — 


g1zZ 
314 


313 


DXXX1=Dxo cei 2 a me Ca = 
MMM=1 

CCMEUTE THE CONSTANTS OF THE MOTION UI , VI, ALP1 AND ALP2. i 
AP 1=~-SQRT(S*(S+1.)) *DXO*LEDS/ABS (B) 
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ALP1=ATAN (CP1) +PI 


~ 320 AP2=SORT (S* (S-T.)) ¥DYO*LEDS/ABS(B) . oe 


rn 


BP2=YO-S*DXO*LEDS/B 
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1-UL*UL*COS (THX2) /(S+1.)) -UL*VI* (SIN (THP12) * (1. /SQRT (S*5-S) 
241. /SORT (S¥S+S) ) SIN (THN12) * (1. /SQRT (S*S-S) -1. /SQRT(S*S#5) ))) +1. 


THZ=B*ZZ11 * la? roll , . cr 
COMEUIE THE ROVATING COORDINATES AND THEIR ADJOLNT NORMALIZED 


for sey vey hne* 


ie 


_—_ a THY USE SEI 


a\e Rena 1 s0XuS%e 
tre,0re,e6 teak 
ibaa 59 a 


L4 (yen | \ECRL RG aeCRs* UI (. reepeg) cee (@\2GRL80YA*e-OR}) PAVE 
( (see) VEU se 26s 19 OTd O10 * (fe) 8) shee So ee a 
_—_ il Sib 

s2it—1 auA=20 
Vee ol 
ralae peXt= 
SSIA+TEST = 


— _ = ~~ 


sa 


——— a om 


rare «i wk a 
rir ene ist a 
Yiit=2H} 4 - 


7 


gaeyakeww TAT ee OWN Ts woeNTeTO Taree arn SF — 
suewenpyy ©: A Tiyasdes.0i+ (x47 Hoy -OYAA SOR 
‘Sys WEB) rae sey at =2)\e) F002) *BFeS Se (EE nivy mid 
Bree gobs) 914) 0108209 tf 92 fo Ry up! ae cE ide 
¢ wt) aie Lied ; a 
Je (SUNT) 29I~ (HOUT) 2OD= see bin 


yeti at-a) \aver¥—{ 142) \aerly 


1yriv 
GEN. 1)* faut) UPS) *L Vel 9 (fee 
oy 2\ st eeeney woe» 180m Z 


— 
GIShT UMA BATAWI UO) DUTT 


eae 


44 


 TTSTT+HAO 


46 


MOMENTA XXX,YYY,DXXX,DYYY, 
YYY=VI*COS (THY) -UL*SORT (S/(S#1.)) ¥SIN (THX) ¥ABS (B) /B 

XXX=UL*COS (THX) +VI*SQRYT (S/(S-1.) ) *SIN (THY) ¥ABS(B) /B 260 
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XX11=XXX*COS (THZ) -YYY*SIN (TH) ent: eae ae as 
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DVAZ (MMM) =PVZ11 

MMM=MMM+ 


LF (MMM-J) 45,45,46 
MANNN=O 

PLOT FRAMES , AXES AND WRITE SCALES ON AXES. FOR PLOTTING AX , AY 
AND CVAZ VERSUS AZ AND AX VERSUS AY ONCE FOR THE ANALYTIC 
TRAJECTORIES COMPUTED ABOVE AND FOR THE RUNGE - KUTTA 


TRAJECTORIES COMPUTED BELOW. 


CALE PLOT (0.0; 9.25, 2) 

PRUE, (PLOU(5+s0. 9525, 2) 

CRUEL PELOT 16. 0, 0.0; 2) 

SALE PLOT. 0 70.0, 2) 

CALUeo EBOLd 0. 6,0. 75,0. 12, ANALYTICAL’ 0.0, 70) ase 
Chpe ATS (055 0. 5, '2 (RY, 9,55 0,0. 0, ise te Oe eo 0) 

Cpr, SVEBOL {0.43705 35, 0.08, 0.02 0.053) 
CALLS YNSOL(ls tae 06 35,0208, 81.5" 0.0633 
CAUd7 SIMBOi 24 43,0.325,0.080, 13.0%, 00,3) 
CAD) SY" 80134 43,02 35,0.08,°4.5",0.0,3) 
CALL SI ReeL (4. 43701.35,0.08, 16.0" , 05043) 
Ci Ped OL 1 o. B 3G ed pO ey oo plays) 
ChLis ARIS10S 5,0. 75, * KACO) oy S42 0F 90805 J 0E10s 1 20520203 
CALL SITRBOLI0. 19 ,0..75,0.08, "=32 524500078) 

CALL SYMBOL (9.119, 1% J5,0.08, *= 05 81602 074) 

Cat. SYeR0L(0..25,2. 75,0. 08, 70.84 ,0,0,3) 
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Cat Maser SS OF Stouts TMT joa Re Pe 
CALL SYMBOL (3.1,8.9,0.12,"Y (CM) ',0.0,5) 
CPL HOLS Wey 3 0,2 8 10.0.0, 1, 0810. 1.0, 20.0) 
CALL SYMBOL (5.25,7.00,0.12,*X (CM) ',0.0,5) 

___ CALL SYMBOL (0.6,5.0,0.12, *RUNGE-KUTTA',0.0,11) 
CALL “CPRETE (4, 87517 7.0,0.0,000., 10879, 10079,0.5) 
CALL PLOT (-10.5,0.75,-3) 


GOTO 23236 
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MJM=1 

D1Z11=SQRT (2. *QM*VO) 
D1X11=DX00*D1Z11 

D1Y11=DY00*D1Z11 

AQM=-2. *QM*AK 
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~~ )1Z712=D1Z114+0. 5*H*D2Z211 © 


BT12=2.*B*¥Z12 
CiB=605 (BT-12) 
S$12=SIN (BT12) 
D2X 12=AQM* (X12*C 12+¥12*512) 
D2¥12=AQM* (X12*S 12-Y12*C12) 


HOU D= NOU eB 20* KX T2*Y 12*C 12 (2 128 122292 Fk 12) FS 72) 


113=111+0.5*H 
RVS=MVIF0SS*H*D 1X12 

VIS =Vala es os SHADY 12 

@ 03=2 UT40. Oe HAD 12 12 
Drei seD tel ihOr, OF HD 2 12 


DIVII=DIVI1+0.5*H*D2Y12-°°°¢ & #23 
D1Z13=D1Z11+0.5*H*D2Z212 

BT13=2. *B*213 

C13=COS (BT 13) 

$13=SIN (BT13) 

D2X 13=AQM* (X12¥*C 13+Y 13*513) 


DZY 13=AQ0M* (X13*S T3-¥ T3¥*C13) 

D2Z13= Fo pode oer ioent Ri S*213 e593) 
T14=111+H 

¥14=K 174H*D1X 13 

Yieav Tithe D1Y13 


240=2 147tH*#D 1213 
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DIX14=D1X114+H*D2xX13 
DIY 14=D1¥114#H*D2Y13 
D1IZ14=D1Z211+H*D2213 263 
BT 14=2. *B*Z 14 

C14=COS (BT14) 

-S14=SIN (BT14) 

D2xT4=AQn* (XTU*CTU+Y1G%STQ) a = 
D2Y 14=AQM* (X14*S 14-Y 14*C14) 

D2Z 14=AQM*B* (2.*X14*Y 14*C 144 (YIG*Y 14-X14*X 14) *514) 

XA VA=7 14 

YUASA 1 

PATR=E FA 


“hhi=e ities. aia ; ; 3 
K11=X114 (H/6.) *(D1X114#2. *D1X1242.*D1XK 13+D1X 14) 
Y 11=Y¥ 11+ (H/6.0) * (D1Y¥11+2. *D1¥124#2.¥*D1Y 134D1Y14) 
Z11=211+ (H/6.0) * (D1Z114+2. *D1Z124+2.*D1Z134+D1Z14) 
611 DIX11=D1X114(H/6.0) * (D2X114+2. *D2K 1242. *D2X 134+D2X 14) 
DIY 11=D1¥114+(H/6.0) *(D2Y11+2. *D2Y12+2. *D2Y134#D2Y14) 


D1TZ11=D 1211+ (1/6. 0) ¥(D2Z 71142, ¥D2Z 1242. ¥D2Z 13+D 2215) 
BT11=2.*B*Z11 

Ci i-cosar 1) 

S$11=SIN (BT11) 

D2X11=AQM* (X114C11+¥11*S11) 
D2Y11=AQM* (X11*S 11-Y11*C11) 


D2Z TT=AQM*B* (2. *X TT¥Y 1 TEC T 14 (YT T¥Y TI-X1 1¥X 17) ¥ST7) 
Be Wiz Tio: 
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MIM=MIM+1 

AX (MIM) =K11*100. 


Miho ey e100. SCO : - 

AZ (MIM) =211 

DVAZ (MJ M) =D1Z11/¥Z0 

7772 LF (MdM-J) 3,3,33 
33 MANNN=MANNN+1 

GO TC 3333 
eos ool -ddo0t 1 

Jao =300 189 

AE od 17 =050 

WAT AEE Palesiie o 

AX (000 1) 2-16 

Aveda Jdoy=0. 0 


a wid =-1.6 
AYdodoZyH0. 6 
Helv 10.0 
HVAT Jado =0.3 
IF (MANNN) 811,811,812 

O1c ALT ELOT (0. 0,4, 25,-3) 


Bott CAG riINE (AZ PAX; OUUUylsU,0) 
CRUE MPHOTCO. Ope. y= 3) 
Cult 2INE (AZ, AY, I9100,51.070) 
CALI ELOT(1120,0.0,-2) 
CALL LINE(AZ,DVAZ,JJJ0,1,0,0) 
CRUL PLOT (0.5,—12.0,—3) 


ee ert LINE (AL AY;Jd00,1,0,0) 
IF (MANNN) 821, 821,822 
Bo SP CAdh (ELON (=1145,0<0,—3) 
GO TO 567 
822 MANNN=0 
GAUL UPLOT (11.0725. 0,-3) 
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555 CONTINUE 


223 CALL FLOT(0.0,0.0,999) 
WRITE (6,6666) 
6666 FORMAT (17H**PLOT FINISHED**) 


264 


STOP 
__ END ; ; 
LS 


SECOND PROGRAM 265 


PROGRAM FOR TRACING PARTICLES ALONG THE TWISTED LENS STRUCTURES. 
PLOTS ARE PRODUCEO FOR THE x=-z. x 


PLOTS ARE PROOUCED FOR THE X-Zy Y-Zs PX-Zs PY-Z» PX-PY AND X-Y 


PROGECTICNS OF THE MOTION. 

NUM = NOs OF (XOsYG) CCORDINATES CONSIDERED. 

KLJM = NOs OF (FXOsPYC) MCMENTA CONSIDERED FOR EACH INJECTION 
POINT (XOsYO). 

QM = CHARGE-TO-MASS RATIC (CCULOMB/KGe)e 

AR = RADIUS OF APERTURE. 


M sN = DESIGN PARAMETERS SPECIFYING LENGTH OF THE LENS. BOTH M oN 
GDD OR BOTH EVEN INTEGERS AND MDND>O. 
NNN = PARAMETER SPECIFYING INITIAL CONDITION 
NNN = PARAMETER SPECIFYING DIRECTION OF TwIST OF THE ELECTRODES. 
NNN=0 FOR COUNTER-CLOCKWISE TwIST(8B POSITIVE)s 
NNN=1 FOR CLOCKWISE TwIST (B NEGATIVE). 


MAAA = PARAMETER SPECIFYING THE STRUCTURE CONSIDERED. 
MAAA= —1 FOR SINGLE TWISTED LENS» 
MAAA= 0 FOR THE STRUCTURE CONSISTED OF TWO OPPOSITELY 
TWISTEO LENS@ SECTIONS. 
MTRA = 020 OR NEGATIVE THRE X-Y AND PX-PY CURVES ARE PLOTTED AND 
TRE EXIT CGORDINATES AND MOMENTA ARE MARKED. 


See ee Silom Bee, an lah (aver iar (a) tei vay fejltel ieiial tepiel layed <a) 10) 


MTRA = POSITIVE ONLY THE EXIT COORDINATES AND MOMENTA MARKED. 

3 = PARAMETER CEFINING THE ACCURACY IN COMPUTING THE EXIT 
COORDINATES AND MOMENTA AT THE EXIT PLANEeIT IS THE NUMBER 
OF DIVISICNS OF THE TIME INTERVAL BETWEEN THE TRAJECTGRY 
POINTS AS THE PARTICLE APPROACHES THE LENS EXIT. 

vo = ACCELERATING POTENTIAL OCF THE PARTICLES. 


AM NUMBER OF TRAJECTCRY POINTS COMPUTED PER PERIODIC LENGTH 
OF TWISTe 

DIMENSION AX(2090)sAY(2000)s4Z(2000) »PX( 2000 )sPY( 2000) 

DIMENSION RRPX(30)sRRPY(30)»RRPZ( 30) 

CALL PLOTS 

CALL PLOT(020120+-3) 


REAC(Ss2001 )KLIM 

FORMAT( 12) 
READ(S5s,101) JS eMeNsAMoVOeARs GM 

FORMAT (212e4E12¢€&) 

COMPUTE IMAGING VfLUE GF S UGING MsNe 
BAM=M 


BAN=N— 
EAP=EAM*PAN+BAN¥BAN 
BAND=FBAM*BAM—-BAN*¥BAN 
SO=EBAP/ BAC 

PI =321415927 

AL =2e*PI*AR¥*SQART(10¢0) 


“COMPUTE LENGTH OF LENSe © 


ALLENS=SQRT(BAP/9 20) ¥AL 

B=2e*PI/AL 

AK =8*B*(VO/SO0) 

V =AK*AR*AR 

RFAO(5s102)NUMsNNN the al - 


“READS. 102) MAAA o MIRA 


102 


FORMATC 212) 
ANN =NNN 


AB =COS(PT*ANN) = 
SL Gi. ERAMES Go AXES! AND, PRTNT SCALE? ONO THE AXES AND DASHED CIRCCE 


CN THE X-Y PLOT TC REPRESENT THE GUADRUPOLE APERTURE e 
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PLOT (0200902552) 


PLOT (620 .9.2532) 
PLOT (6203902082) 
PHOTO. 0505 052) 


AXTS(1e0s2o7Sot “ol e4oOe0eOe le OE 0s 1le0s 2000) 


SYMBOL (5 0259207590129 "X(CM) "0 0005) _ 


CTRELE(45.27532.75%0 Co 36000108750108752005) 
AXETS( 30090075 0% %515400290000100E10010032040) 
SYMBOL (361 94 e630 e125" V(CM)%300025) 
AXTS€032095005* %515450990005 100E10010002020) 
SYMBOL (3 0159 eSsCel2s *PY/SPZ° 900085) 

AXIS (1209720! Seo l140020 20051 e0E10 6120220040) 


SYMBOL (502557005 Ce 12s *PX/PZ* 602005) 
PLOT (020312202-2) 

PLOT (0 200922502) 

PLOT (620:39225£.2) 

PLOT (620202032) 

PLOT (02090 e032) 


AXIS 50065. "Zi CM) s-555.000e0s 1 SUEID s 16052050) 
SYMBOL (0.4 3s 0235200 08s "D060" 50003) 

SYMBOL (12590203540 008s 98" 900001) 

SYMBOL ( 2 0520 035902089416" 502002) 

MMi os Ge 45906085 "24%,050s2) 

SYMBOULN 4 0550635. 0ls CBs 9327505052) 


SYMBOL (5 0520035202085 "'40% 00082) 
AXITS(0e5sDe75Ss*X(CM) % 980400 eID e Os Le 0E10910092000) 
SYMBOL (0 e619 eOe75s0el&s *—-1e06*% 202054) 

SYMBOL (0 0199107520 0eChs*—008* 800054) 

SYMBOL (0 e269 20752 00l&s*000% 500023) 

SYMBOL (0 6263327590 068s 40 0B* 300093) 


SYMBOL (0 0260407590 008s *106% 200093). 

AXIS (0 055500 oF VY(CM)% 54 54 ce Oe GO Ole AELOs 10022020) 
SYMBOL (0219352030208 e'—1 06? 302004) 

SYMBOL (0 0199600906089 '-048% 302094) 

SYMBOL (0019 e27eO090e 08s 9020" 902083) 

SYMBOL (0215 28 e090eC8s "00 8%2000 +3) 


SYMBOL (021996 Ge0e0os “leom yO. Gea) 
PEON COs S ee 6 1553) 

PLOT (5252207592) 

PLOT (52527203) 

PLOT (025070022) 

PLOT(102e0s02e0s-3) 


PLOT (0 «0539 «25s2) 

PLOT (6200902592) 

PLOT (60020 0052) 

PLOT (06010 00.2) 

AXLS (00530 050" Z(CM)*® 9-595 e000 091 e0E10s 16022000) 
SYMBOL (0 04320 03550008s*000' 200093) 


SYMBOL (10500 03590008s "8B" 500001) 
SYMBOL (2 052 023500008 416% 800002) 
SYMBOL (3259023590008 s*24" 100082) 
SYMBOL (4 050023590008 s' 32% 500002) 
SYMBOL (5 0520 035900C8s 240% 200092) 
AXIS (00500075 eo PX/SPZ# 959400290 0s be OE10s 10042060) 


SYMBOL (00190007590 0C8s'—002® 100094) 

SYMBOL (00199107590 00ls*—-Oel® 900094) 

SYMBOL (00199207520 0C8s 2000% 900093) 

SYMBOL (001993075 90eCEs*00l* 200083) 

SYMBOL (0 0199407590 e08s 002% 300023) 

AX1S (005 95e009*PY/PZ® 95940099000 Le 0E109 10092000) 
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CALL SYMBOL (0619550000 e08s*—Oe2® o0c0e4) aces 
CALL SYMBOL (0 .19s6.0s0.08s *-0.1%s000 04) 
CALL OSYMBOL (0 226 »7 0020008 s "020% s0e0e3) oy 
CALE SYMBOL (0226980020008 s"00l%s0 0023) 
CALL SYMBOL (0 226 05 20200C8Bs "002% 00093) 

WEAS HPUIOT (0S), 2.75.3) 


CAINS PLOT (54552075 .2) 7 = 
SAAT (56567003) 

EC ACAS FPO (0S 997 610s 2 ) 

GA PLOT (—10 0 .—-12.0.—3) 

DO 5566 LTIL=1+KLUM 


9566 READS, 103) RRPX( TIL) »RREYCTIL)» KRPZ{(TIL) 


DO 666 IEX=1+NUM ; : Yi'eea 7 
READ (5,1023)X0sYC,20 
DO 666 L[ABC=1.KLYUM 
RPXO=RRPX(IABC) 
RPYO=RRPY (IABC) 
RPZG=RRPZ(TABC) 

~103 FORMAT (3612.6) 
PZODY 360 

PZ0 =RPZO 

ZOO=0%0 

THOOOD=0 0 

r= 


CALCULATE S FOR THE PARTICLE ACCORDING TO ITS INJECTION CONDITIONS. 
1033 PXO=RFX0*PZO 

PYO =RFYO*PZO 

PZ=PZO*(1.64+8*( XO*RPYC-YO*RPXO) ) 

DXO*=FxX07E2Z 

DYO =PYD/PZ 


S =SO*PZ*PZ 
GMK =2¢*GM*AK 

FX =SGRT(QMK*(S41e)) 

FY =SCRTCQMK*(S5-1.2)) 

VZU=SORT(2e*QM*VO ) 

START CALCULATICN OF TRAJECTORY PCINTS USING ANALYTIC SOLUTION. 
“PELE LAMEV EO) | PYCAT Tee mies oe ea Bs gi a 
UI=SQRT( (XG—SGRT (VO¥S/AK ) #FYO )*¥*24+(VO*(St1 2) *PXOXPXO/AK) ) 
VI=SORT( (YO-SQRT(VO*S/SAK ) ¥PXO )**24+(VO¥(CS—1.) *¥PYO*PYO/SAK) ) 
AP1=—SQRT (VO*(S41e)/AK)¥*PXO 

BP1=(XO-SORT(VO*S/AK)*PYO) 

IF (8P1) 30942104311 


~~ 310 IFCAP1)212,232135314 


312 ALPI=—-PI/2. 
GOP Los 29 

314 ALP1=PI1/2- 
GO phe es 20 

313 ALP1=0-0 
Gory (Oes2 @ - 

311 CPL=API1/6P1 
ALP1=ATAN(CP1) 
GO ti 20 

309 CP1=AP1/6P1 
ALPI=ATAN(CP1)+PI1 

JO APS=SORTCUVI*(S—Te )AZAKY¥PYCT 
BP2=YO-SQRT (VO*S/ AK) ¥PXD 
IF (BP2)4099410+411 

410 TF(AP2)4122413+414 

412 AL P2=-P1/2e6 
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ALP2=PI/2. 

GO TO 420 
ALP2=0.0 

GO TO 420 
CP2=AP2/BP2 
ALP2=ATAN(CP2) 


GQ TO 420 
CP2=AP2/8P2 
ALP2=ATAN(CP2) 4P1 
HO=H 

KK=0 

eo. 


ZFI1=6 36 
ATHX=FX*TT 
ATHY =FY*TT 
CX=COS(ATHX) 
SX=SIN(ATHX) 
CY=CGS(ATHY) 
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SY=SIN(ATHY) 
FACTG=SQRT(VO/AK) 
ATI=A€x 
A1L2=AB*ESY*SQRT(S7( 


Sz ier 


A13=(SX*SOQRT(S4#1.)-S*SY/SGRT(S—1.-)) ¥FACTO 


Al 4=( AB*SORT(S)¥*( CY-CX) ) ¥FACTO 


~ A21=-AB¥EEX#ESOART(S/(S#1.2)) 


A22=CY 
A23=-A14 


A24=((SX*SSSQRT(S+1 2) I-SY*¥SCGRT(S—1le) )¥ FACTO 
A31=-SX/SOQRT(S#+1.) 


A32=9.90 
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AS 3Z=CK*EACTO | 


A34=-Ace1*FACTO 
A41=020 
A42=SY/SQRT(S-1e) 
A43=AL2*FACTO 
A44=CY*X FACTO 


~AXLI=ALIXXOFA1L2*YC+A1 34PXCHALG*PYO | 


AYL1I=A21 *®XOF+A22*Y OAL I #PXTHAL4*PYO 
PXLLI=CABSLI*XXOFAZ2*YOCAIZ¥PXOFAZ4*PYO) SFACTO 
PYLI=(CAG1L¥XO+A424*YO4A43 ¥P XO4+A444 *%PYOISFATCTO 


THX=ATHX+ALP1 
THY =ATHY +ALP2 


THX 2=2 6 ® THX 

THY 2=2-%*THY 
THP12=THX4#THY 
THN 12=THX-THY 
CAGQARZ 11 

THX 20=22¥*ALP1 
THY 20=24¥*ALP2 
THPUG=ALPI+ALP2 


THNOG=THPOG-2 e* ALP] 


Z1L1=VZ04TT¥*¥(1-+B*( 


XO*XRPYO-YO*¥RPXO) +0 eS5S¥*B*SQRT(AK*S/VO) *(CUT*UTSIS 


Lele) )—CVIEVIS(S— 10) ))940025*B*(SQRT(S/(S— Le) X* 3) FVIRV IFC SING THY 2) 


2-SIN(THY20) )-SQRT ( 


- 3BKUTKVI*® (COS (THP12) —COS(THFOO) -COS (THNI 2) +COS(THNOO) 07 


4SQGRT(S*S—-1.6) 
CHECK THE AXIAL DI 


S/( Stilo )*#*3) UI #UL*(SINCTHX2)-SIN(THX20) ))+0e5* 


STANCE TRAVERSEC BY TRE PARTICLE.’ 


IF (Z11-ALLENS )601 56014602 


TT=TT-HO 
DECREASE TRE TIME 
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PARTICLE BE FINALLY EGUAL TC THE LENS LENGTH WITHIN THE DEGREE OF 
Cc ACCURACY SPECIFIED BY Je 


HO=HO/2.0 
KK=KK+1 
TT=TT+HO 

Oe PACU—KK US 30.830 ,920 

601 T=1+1 mii) ih) 
TT=TT+4+HO 
THO 1=B*Z114THOOO 
CCTH=COS{( THO!) 
SCTH=SIN(THOTI) 
AXI=AX11 
AYI=AY11 a seer es zy 
Se a a ie | 
PXI=PX11*PZO 
PYI=PY11*PZO0 
PZI1=P2-B*( AXI*PYI-AYI*PXI ) 

c STORE TRAJECTGRY POINTS IN TERMS GF THE FIXED COORDINATE SYSTEMe 

—_ AX (TY =AX11*®CCTH-AYI1%*SCTH ~ i 
AY(C1)I=AX11*SCTHEAY11¥*CCTH 
PX(L)=(PX11*CCTH-PY11¥*SCTH) *PZO/PZI 
PY( IT )=(PX11*SCTHEPYLIL*CCTH) ¥PZO/PZI 

AZ(1)=Z11+4Z00 
6006 FOURMAT(SXe5(E1023 »3X)0I4) 
- Sito s20 Gebreiestop- ie °° UOC ae 
530 MAAA=MAAA+1 
ITF (MAAA )5305% £3052 5306 
FGR THE STRUCTURE CCNSTRUCTED CF TwO OPPGSITELY TWISTED SECTIONS 
CONSIDER THETRANSEVERSE COORDINATES AND MCMENTA AT THE EXIT PLANE 
OF THE FIRST SECTION IN THE ROTATING COORDINATE SYSTEM AS THE 
INJECTION CONDITICNS CF THE PARTICLE FOR THE SECOND TWISTED SECTION. 
5305 XG=AxI 
YO=RAYI 
ZOO=AZ1 
PZ0=PZI 
VO=VO¥#PZI*PZI 
RPXO=PXI/PZ0 
RPYO=PYI/PZO 
THOCO=8* Z0G 
B=-B 
HO=H 
GO 1o 1093 
Biot Lisie: | 
f=. 42 
f CREAT SCALES FOR PLOTTING THE TRAJECTORIES 
AX(1I1)=-0.016 
AX(112)=0.008 
AY OC TI=AX UT 1) 
AYCTSP=AxXCI2) 
PX(I11)=-O402 
Px tL ei=O. 1 
PY(11)=-042 
PNAC LO) =0\6 1 
AZCTAS=050 ee 2 Seer nccerccei = _ 
— AgdtI2)=0-08 ; we 
GAC PEG .000 el 50-22 
ABAX1=(AX(1)-AXCT1) )ZAXCT2) 
ABAX2=C(AX(IT)-AXC 11) )ZAX (12) 
ABAYL=C(AY (1) -AYC TI) IZAYCT2) 
ABAY2=(AYC(L)-AYC IL) SAY CTA) 
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MARK. ENTRY. ANDO EXT -“TRANSEVERSE COORDINATES ON THE X-Y PLANE. 
CALL SYMBOL (1 ABAX1 sABAY1.0. 045003 20e0s-1) 
GALA SYMBOL (AB AX2 sABAY2s0204.011+0005,-1) 
IF (MTRA) 10551053106 
DOS MGARE EINE CAX s AY sf sls 0) s.0)) 
LUG CALL PLOT(0 00 94.6 tA ie} 


ABP X1= (PXC1)-PX(11))/PX( 12) 
ABPX2=(PXC1I)I-PX(11))/FX(12) 
ABPYIL=(PYCLI-PYCT1))/PY( 12) 
ABPY e=(PY(T)I-PY(C1I1))/PY( 12) 
MARK ENTRY AND EXIT TRANSEVERSE MOMENTA ON THE PX-PY PLANE. 
CALL SYMBOL (AEPX1 »sABPY150. 044003%020s-1) 


CALL SYMBOL (ABPX2 sABPY250204401120.02-1) 
START PLOTTING THE TRAJECTORY PROJECTIONS. 
TF (MTRA) 11591155116 
115 CALL LINE (PX sPY 51912050) 
CALL PLOT (-0.557275s-3) 
CALL LINE (AZsAXo1013090) 


CALL PLOT(0-0+74425.4-3) 
CALL LINE (AZ 5AY 91912050) 
CALL PLOT(10.-05;020,-3) 
CALi LINETAZ PY 6101200) 
CALL PLOT (020.-4- 255-3} 
CALL LINECAZ sPXoisls Qs 0) 


CALL PLOT (-10255-12075s-3) 
GO Ne ees 
ite CALL PLOT (—1205=5 20 s—3) 
666 CONTINUE 
CALL PLOT(02059202999S) 
WRITE(6s7788) 


7788 FORMAT{(10X.9***xx PLOT FINISHED *#* xx *) 
3333 FORMAT(10X%»s "HHRHHHHH ® ) 

Sieh 
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APPENDIX F 
SPATIAL MODES OF THE MOTION ALONG THE TWISTED STRUCTURE 


In order to express the analytic solution of the particle motion 
along the twisted structure in terms of the fixed (X,Y,Z) coordinates, 
equations (2.31) are substituted into equations (2.7). The X co- 


ordinate of the moving particle will be, 
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[C, cos F in eearcal C, sin F,]cos BZ 


ie Same n . : 
- [C., cos Fh mor ni Ale! C, sin F,Jsin BZ (F.1) 


while the Y coordinate will be 
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Xx Sin BZ Py COS ieZ 


[C, cos F, pe = i ee Cy sin Fo ]sin Bz 


S n 
+ [C,, cos F., oY eras | (-1) C, sin F, Jcos Bz (E.2) 


where , F, = ft + a, and F., = fat + do. Equation (F.1), after the multi- 


plication of the trigonometric terms becomes , 


a - : ie 04 


. : ia : 
gauTouat2 GiTelwi WT Ay aorton at ea eH | ea 
rr | — _ Cae 


i — 


re is 5 patee ‘ at a 
avhtee efotaveq ong Yo sotruloe ofayfons oat aoe oe webt6 a 
 estootivocs (5. %yX) Seat? att 70, enrved nt ures, bere twt ait oat 
~0o % ott (5.0) enoPetpe oda) baj et treduz fab ({£. $) “not supe = 


pad fttw sist? 18g ror ont + sie vo ; 
. 7 Ye 


Te eee a 
santzy>S8@o7 8K 
is ee 


7 ' 


ane ane 
<4 2ont§ nla d “(I-) +9 209.42) © 7 
a a * i = [ "ae Gi - - 


(1,3) sa ntel 1 ote po M- \13\, 3 2921 


 $h 200 y Se Fe x =¥ 


. 


7 
ah nal.) ate od "J +\;" Pe pl 


on an 
eat) 
a vee 
(8.7) sq 2oa[ 7 nte ,9 "(t-) oi 3 s mo 1+ | 7 
- a as _ 
7 : rs 0 
tetum one save «(1 4) motieupd .ge + Jet = cP link iy penis 


7 
_ ae att 7 


_<seneaed tes 3 ote : 


Pa 


ra fs 


X = C, cos(f,t + a, )cos Bz 


S 
ae oy Gana Cy sin(f,t + a, )sin BZ 


- Cy cos(fpt + a, )sin Bz 


S 
Ne = a C, sin(f,t + 5) COS Bz 
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Equation (F.2), after the multiplication of the trigonometric terms 


becomes 


Y=C cos(f,t + a, )sin BZ 


i: ape 


1 sin(fjt + a, )cos BZ 
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+ C, cos(fpt + Gy )COS BZ 


+ FSS (-1)" 6, sin(fot + a)sin az 


] ‘ { 
7 Ci [tsin(f,t + a4 + Bz) - sin(f,t + a - Bz) } 


- he Sey init +a, + Bz) + sin(f,t ta, - Bz)}] 


] 
7 Coltcos(fat + Oy + Bz) + cos(-f,t - a, + Bz) } 


+ (-1)"\ {<> {cos(-fpt - G + BZ) - cos(fot + a, + Bz)}] 


x Cy{1 + (-1)"\faa? sin(f)t - 82 + 04) 
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